A/m/9 3 ^^ 


17 


6 00166 4649 


DOE/NASA/0029-1 
NASA CR-1 65380 




Wind Flow Characteristics in the Wakes oi 
Large Wind Turbines 


Volume I— Analytical Model Development 


W. R. Eberle 

Lockheed Missiles and Space Company, Inc. 


September 1981 


Prepared for 

National Aeronautics and Space Administration 
Lewis Research Center 
Under Contract DEN 3-29 


for 

U.S. DEPARTMENT OF ENERGY 
Conservation and Renewabie Energy 
Division of Wind Energy Systems 


LIBRARY COPY 

NOV 4 1981 

LS.NGL2Y RGSGASCH CENTER 
LIBRARY, NASA 
HAM.nON, VIRGINIA 


NF01568 



NOTICE 


This report was prepared to document work sponsored by 
the United States Government. Neither the United States 
nor its agent, the United States Department of Energy, 
nor any Federal employees, nor any of their contractors, 
subcontractors or their employees, makes any warranty, 
express or implied, or assumes any legal liability or 
responsibility for the accuracy, completeness, or useful- 
ness of any information, apparatus, product or process 
disclosed, or represents that its use would not infringe 
privately owned rights. 



DOE/NASA/0029-1 
NASA CR-1 65380 


Wind Flow Characteristics in the Wakes of 
Large Wind Turbines 

Volume I— Analytical Model Development 


W. R. Eberle 

Lockheed Missiles and Space Company, Inc. 
Huntsville, Alabama 


September 1981 


Prepared for 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 
Linder Contract DEN 3-29 


for 

U.S. DEPARTMENT OF ENERGY 
Conservation and Renewable Energy 
Division of Wind Energy Systems 
Washington, D.C. 20585 
Under Interagency Agreement EX-76-1-01 -1028 






FOREWORD 


The material in this report documents the efforts of the 
first phase of an overall program conducted by Lockheed Missiles 
and Space Company, Inc. , of Huntsville, Alabama, under NASA 
contract DEN3-29. The first phase of the program called for a 
computer model of the wind turbine wake, while the second phase 
dealt with wake measurements on the Mod-OA wind turbine at Clay- 
ton, New Mexico. 

The computer program of this report is considered an initial 
calculation procedure for the wake development of wind turbines 
because additional expansions and refinements are possible to 
make the program more comprehensive. Areas for further consider- 
ation include additional methods for calculating the effects of 
atmospheric turbulence, non-equal vertical and lateral wake growth 
rates, and power recovery factor based on real turbine power char- 
acteristics . 


iii 




V 


TABLE OF CONTENTS 

Page 

1.0 SUMMARY 1 

2.0 INTRODUCTION 3 

3.0 DESCRIPTION OF WAKE MODEL 5 

3. 1 Approach 5 

Turbulence relationships. 5 

Wake geometry 6 

3 . 2 Data Input 8 

3.3 Calculation of Program Constants 11 

Recalculation of input geometric parameters 11 

Calculation of the initial wake radius. . . 12 

Calculation of wake growth due to ambient 

turbulence 12 

Calculation of wake characteristics at 

the end of Region 1 17 

Downwind extent and radius at end of 

Region II 21 

3.4 Wake calculations in Region III 22 

Wake growth in Region III 22 

Turbine power factor 25 

3 . 5 Wake plots" 29 

3.6 Calculation of Wind Speed Profiles 29 

Coordinates 29 

Region 1 31 

Region III 31 

Region II 32 

3.7 Calculation with Ground Effect 32 

3.8 Calculation of Wind Speed 33 


V 



4.0 COMPARISON WITH PREVIOUS MODEL 51 

4.1 Wake Growth Rate Due to Ambient 

Turbulence 51 

4.2 Support Tower 52 

4.3 Calculations for Region I 52 

4.4 Wake Growth in Region III 53 

5.0 SAMPLE PLOTS 55 

5.1 Wake Plots 55 

Abramovich solutions 55 

Wakes in ground effect 57 

5.2 Effect of Principal Wake Parameters 58 

Effect of ambient turbulence 58 

Effect of initial velocity ratio 59 

Effect of turbine rotor hub height 60 

Effect of the power law exponent 60 

6.0 CONCLUSIONS 143 

7.0 REFERENCES 144 

APPENDICES 

A - LIST OF SYMBOLS 145 

B - PROGRAM LISTING 152 

C - WAKE CALCULATIONS ON THE TI-59 CALCULATOR . . .177 

D - ALTERNATIVE ALGORITHMS FOR CALCULATION OF 

THE DOWNWIND EXTENT OF REGION 1 196 


vi 



1 . 0 SUMMARY 


One of the proposed approaches to reducing the consumption 
of fossil fuels Cin particular, petroleum and natural gas) in 
the United States is the use of large wind turbines. It is 
expected that such turbines will be constructed in arrays in 
large fields of wind turbines. To minimize cost, it is 
desired to space the turbines as close to each other as 
possible. However, the turbines must not be spaced so closely 
that a turbine lies within the wake of another turbine and 
thereby the output power of the downwind turbine is reduced. 

In order to determine the appropriate spacing between 
turbines, a research program was initiated to characterize 
the recovery of the wake behind a large wind turbine. The 
research program had two aspects. The first was the develop- 
ment of an anlytical model of wake recovery downwind of the 
turbine. The analytical model was developed to calculate the 
wake properties as functions of the downwind coordinate and 
to generate wind speed profiles (i.e., wind speed as a function 
of the lateral coordinate for selected altitudes and wind 
speed as a function of the vertical coordinate) at selected 
distances downwind of the turbine. The inputs to the model 
included wind speed, ambient turbulence, and turbine geometric 
parameters. The second aspect of the program was an 
experimental program to measure the wake behind the Mod-OA 
wind turbine at Clayton, New Mexico. The Lockheed laser Doppler 
velocimeter was used to make the measurements. This volume 
contains a description of the analytical model development. 

The analytical model development is based upon the 
results presented by G. N. Abramovich in "The Theory of 
Turbulent Jets". The term "jet" is used to describe a flow 
surrounded by another flow of a different velocity, regardless 
if the velocity of the inner flow is less than or greater than 
that of the surrounding flow. Although the Abramovich results 
were taken as the basis for the analytical model for turbine 
wakes, several modifications of the Abramovich results were 
necessary. The Abramovich model assumes that the initial 
jet emanates from an opening. Therefore, it was necessary 
to add calculations to relate the turbine to the initial 
wake as presented by Abramovich. Second, the results 
of Abramovich were developed for flow with no ambient 
turbulence. Therefore, the effects of ambient turbulence 
were added to the Abramovich model for the turbine wake model. 
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The wake growth rate due to ambient turbulence was taken from 
atmospheric diffusion theory as developed by F. Pasquill. The 
method of combining the wake growth rate from the Abramovich model 
with the wake growth rate due to ambient turbulence was taken 
from the work of P. B. S. Lissaman from a prior model for the wake 
recovery behind large wind turbines > 

Third, the concept of a turbine power factor, which is the 
ratio of the power generated by a turbine centered in the wake of 
another turbine to that generated by a turbine in the free stream 
air, was added to the model. 

The mathematics of the wake model was written in a FORTRAN 
computer program. The FORTRAN computer program was implemented on ' 
a PDP-10 computer using the DISPLA graphics package for graphical 
presentation of results. In addition to the computer program, a 
version of the wake program was formulated for the TI-59 program- 
mable calculator. The calculator version of the wake program is 
presented in an appendix. 

Several runs of the computer program were made to illustrate 
the effect of wake variables on wake recovery. The effect of 
ambient turbulence, initial wind speed deficit in the wake, the 
height of the turbine rotor hub above the ground, and the power 
law coefficient for the free stream wind were investigated. 

The results of these effects are presented in figures in the 
report . 

The principal indicator of wake recovery is the turbine 
power factor. Much of turbine operation occurs under conditions 
of neutral atmospheric stability, represented by Pasquill stability 
class D. For this condition and for an initial velocity ratio of 
1.5 [typical for large wind turbines), the turbine power factor 
is approximately 0.5, 0.75, 0.9, and 0.96 at 10, 20, 35, and 50 
rotor radii downwind, respectively. 
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2.0 INTRODUCTION 

One o£ the proposed approaches to reducing the consumption 
of fossil fuels (in particular, petroleum and natural gas) in 
the United States is the use of large wind turbines. It is 
expected that such turbines will be constructed in arrays in 
large fields of wind turbines. To minimize cost, it is 
desired to space the turbines as close to each other as 
possible. However, the turbines must not be spaced so closely 
that a turbine lies within the wake of another turbine and 
thereby the output power of the downwind turbine is reduced. 

In order to determine the wake characteristics of tur- 
bines, a research program was initiated to characterize the 
recovery of the wake behind a large wind turbine. The research 
program had two aspects. The first was the development of an 
analytical model of wake recovery downwind of the turbine. The 
second aspect of the program was an experimental effort to 
measure the wake behind the Mod-OA wind turbine at Clayton, 

New Mexico. The Lockheed laser Doppler velocimeter, (LDV) was 
used to make the measurements. This report presents the re- 
sults of the first phase of the investigation, which involved 
the development of an analytical wake model to calculate the 
wind speed profile (i.e., wind speed as a function of the 
lateral coordinate for selected altitudes) at selected dis- 
tances downwind of the turbine. The principal inputs to the 
model include wind speed, ambient turbulence, and turbine ge- 
ometry and operating parameters. 

The initial wake behind a wind turbine can be represented 
as a volume of flow of lower momentum than the surrounding flow. 
This mass gradually acquires the speed of the surrounding flow 
as it progresses downwind. Some features are similar to those 
of a wake behind a bluff body as described in reference 1, 
while other features are like those of coflowing jets as de- 
scribed in reference 2. There are few published reports on 
the detailed wake behind wind turbines. Templin (ref. 3) and 
Crafoord (ref. 4) have made theoretical estimates of large 
scale effects, that is, without considering detailed flow near 
each turbine. Some experimental wind tunnel work has been done 
by Ljungstrom (ref. 5), while more detailed wind tunnel work 
has been conducted by Builtjes (ref. 6). A comprehensive re- 
view of this work is given by Lissaman (ref. 7) . 

The analytical model presented in this report was de- 
veloped through an evolutionary process. At the beginning of 
the contract covering the work reported herein, a subcontract 
was issued to AeroVironment , Inc. (Lissaman and Walker) for 
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the development of an analytical model and computer program 
for the wake flow. The model developed by AeroVironment was 
based on previous work done for the Swedish National Board for 
Energy Source Development. The report describing the details 
of the AeroVironment model and calculation procedure is refer- 
ence 8. The general principles of the model are also described 
in reference 9, which has a wider distribution than reference 8. 

Upon close examination of the AeroVironment model, it was 
determined that there were certain aspects of the model which 
could be modified to improve the accuracy of the modeling of 
the fluid mechanics of the wake. There were also some re- 
finements in the model which were originally requested, but 
which, in retrospect, should have little influence on the 
overall results. The effects of these refinements were 
deemed to be less than the uncertainties of some of the major 
assumptions of the model. The revised version of the model 
is contained in this report. 

Both versions of the analytical model are based on the 
theory of coflowing jets as presented by Abramovich in ref- 
erence 2. The term "jet" is used to describe a flow emanating 
from an opening into a larger surrounding coflowing flow. The 
velocity of the jet may be greater or less than that of the 
surrounding flow. In the application of the coflowing jet to 
the analysis of wind turbine wakes, ambient turbulence must 
be considered. Thus, for both the AeroVironment model and 
the revised model as reported herein, the main task was the 
addition of the effect of ambient turbulence and the ground 
plane to the Abramovich model, computerizing the results, and 
providing automated graphical output for the wake character- 
istics . 

This report contains a description of the analytical model 
and a comparison with the earlier AeroVironment version. 

Sample calculations using the analytical model are presented, 
and the effects of the principal input parameters on wake 
characteristics of a large wind turbine are explored. Also 
included in the report are appendixes containing a listing 
of the computer model, a version of the model for use on the 
TI-59 programmable calculator, and a discussion of the alter- 
nate algorithms considered in the development of the model. 
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3.0 DESCRIPTION OF WAKE MODEL 
3.1 Approach 

The mathematical model used for calculating the wind speed 
profiles in the wakes of large wind turbines is described in 
this section. The basic concept of the mathematical model is 
that the wind speed profile, in the wake is of different math- 
ematical form in different regions of the wake (as defined 
later) , but this form is uniquely defined in each region at 
each position downwind of the turbine. For each profile region, 
knowing the wake radius at the beginning of the region and the 
wake growth rate in the region, the wake radius can be determined 
at any desired downv/ind location. Then, for a given mathematical 
form of the wind speed profile, the centerline wind speed is 
uniquely determined by the condition that the momentum deficit 
is conserved from the initial wake to the downwind location. 

Turbulence relationships .. - The total effective growth 
rate of the wake is given by the Pythagorean sum of the growth 
rate due to the mechanical turbulence (i.e., turbulence generated 
because of velocity gradients in the flow) and that due to the 
ambient turbulence. The growth rate due to mechanical turbulence 
is obtained by the experimental results of Abramovich (Ref. 2), 
and the growth rate due to ambient turbulence is assumed constant 
along the wake for a given value of the atmospheric turbulence. 
The Pythagorean sum (i.e., square root of the sum of the squares) 
was chosen because the total kinetic energy of the combined 
turbulence equals the sum of the kinetic energy of the mechanical 
turbulence and the kinetic energy of the ambient turbulence. 

Abramovich (ref. 2) has done extensive work on coflowing 
jets in fluids. Much of his work has been based on experimental 
results with no ambient turbulence in the jet or in the 
surrounding fluid. There are two mechanisms of momentum 
transfer in the coflowing jets for which Abramovich has presented 
results. The first mechanism is the result of the viscosity 
of the air acting on the velocity gradients in the flow. This 
mechanism of momentum transfer could be represented by the 
Navier-Stokes equations for laminar flow (although Abramovich 
does not use the Navier-Stokes equations). The second mechanism 
of momentum transfer is turbulence generated by velocity 
gradients in the flow. Because of energy extraction by the 
turbine rotor, the wind speed of the air flowing through the 
rotor is reduced, thus creating a stream tube of flow with a 
wind speed less than that of the surrounding free stream. The 
velocity gradient across the flow creates turbulence. This 
creates an eddy viscosity (ref. 1C) , which increases the momentum 
transfer in the flow. This mechanism of momentum transfer is 
much more significant than momentum transfer in laminar flow. 
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In the following discussion, the parameters calculated 
from Abramovich are termed as being due to "mechanical 
turbulence" because it is the gradients in the flow which cause 
the turbulence which causes the momentum transfer in the cases 
which Abramovich has studied. This term is used to distinguish 
the effects which Abramovich has presented in ref. 2 from the 
effects of ambient turbulence which were not included in the 
effects which Abramovich included in his analysis. 

The wake growth rate due to ambient turbulence has been 
taken from the theory of diffusion of pollutants by turbulence 
in the atmosphere. Because of recent regulations related 
to air quality, the theory of pollutant dispersion has been 
developed extensively over the past few years. The relationship 
between the dispersion of pollutants and the transfer of momentum 
in a wake is given by Abramovich. 

Wake Geometry . - The wake of a single turbine is shown in 
Figure 3-1, as idealized for the unbounded flow Ci*e., no effect 
of the ground). The wake is divided into three regions for 
the calculations. The wake radius in the first two regions 
increases linearly with downstream distance at a rate set by 
the effective turbulence, which is a combination of the mechanical 
turbulence and the ambient turbulence. The influence of the 
mechanical turbulence diminishes and asymptotically approaches 
zero as the wake moves downwind. Eventually, the effect of the 
mechanical turbulence becomes negligible, and the wake growth is 
esseptially determined by ambient turbulence alone.. In order to 
accurately model wake growth with a diminishing influence of 
mechanical turbulence, the wake boundary of Region III was 
calculated by a numerical integration. The details of the 
geometry and flow in these regions will be discussed in the 
following sections. 

In the initial region (Region I) , the potential core is 
that portion of the wake which has been unaffected by the shear 
between the wake and the outer (ambient) flow. Region I extends 
to X = Xy, the point at which the shear due to the outer flow 
has completely eroded the potential core of uniform flow down- 
stream of the extraction disk. Wind speed profiles across the 
wake in this region are not self-similar at various downstream 
locations due to the change in relative size of the core flow 
and the turbulent mixing zone. At the end of Region I, a 
continuous shear layer-like wind speed profile has completely 
developed but is represented by a slightly different functional 
form from that used later in the far wake regime. The transition 
region. Region II, allows for the smooth transition of the 
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completely developed wind speed profile of Region I to that 
used in the far wake, which is self-similar for all subsequent 
downstream locations. 

Region I also includes the expansion of the wake from the 
diameter of the physical extraction disk, R^, to the expanded 
slipstream value, Rq , that is, the slipstream expansion due to 
potential effects. The computer model assumes that this 
expansion occurs at the station of the disk itself so that the 
wake develops from R = Rq at X = 0. Reasons for this assumption 
are discussed in a following section. 

The wake growth rate is identical for Regions I and II and 
is given by a combination of ambient and mechanical turbulence 
as discussed in the following sections. The end of Region II 
occurs at = nX„, where n is derived from the form of the 
wind speed profile in Region I and the form of the wind speed 
profile in Region III as described by Abramovich. The wake 
radius is calculated at the end of Region I and at the beginning 
of Region III based upon a momentum balance with the initial 
wake and the form of the wind speed profiles at each of the 
respective locations. The downwind extent of Region II is then 
calculated from the ratio of these radii and the wake growth 
rate in Region I. 

Region III is the far wake. In the initial part of Region 
III, wake growth is controlled by a combination of mechanical 
and ambient turbulence. Further downwind, the effect of 
mechanical turbulence asymptotically approaches zero, and the 
wake growth is essentially controlled by ambient turbulence 
alone. In the numerical integration in Region III, wake growth 
due to mechanical turbulence is never mathematically eliminated. 
However, its influence asymptotically approaches zero. 

Figure 3-1 also shows the notation to be used for the 
geometrical parameters in the following discussion. Upper case 
R denotes the wake radius in physical units. Notation for lower 
case r is similar to that shown in Figure 3-1, but denotes wake 
radius normalized by the turbine rotor radius, R , . The normalized 
radius of the potential core is r. . The normalized outer radius 
of the wake at the end of Region K is ^ 2 ]^. The notation, r 2 , 
is used for the normalized outer radius'^at any general location 
in the wake. The normalized downwind distance at which Regions 
I and II end are x^j and Xj^, respectively. In the following 
discussion, upper case U denotes wind speed in physical units. 
Lower case u denotes wind speed which has been normalized by 
the free stream wind speed, U^. 
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The following sections describe the computer model. The 
description follows the calculation procedure of the computer 
program sequentially, giving derivations of equations used as 
appropriate. The equations follow the development of Abramovich 
(ref. 2) closely, and no derivation of equations taken directly 
from Abramovich is given here. The reader is referred to 
ref. 2 for derivations of these equations. A flow diagram for 
the computer program is shown in Figure 3-2. A list of symbols 
is given in Appendix A, and a program listing is given in 
Appendix B. The equation number in the following description 
is underlined when the equation is used directly in the computer 
program. In the program listing in Appendix B, the equation 
numbers for equations presented in this section are given. 

In addition to the computer program, a version of the 
wake model was developed for the TI-59 programmable calculator. 
The TI-59 version of the model is presented in Appendix C. 

The TI-59 version of the program uses the same equations as 
those developed in this section for the computer model. 


3.2 Data Input 

The first step of the program is data input. Table 3-1 
shows the input data and card formats for the input data. The 
definitions of input parameters are given in Table 3-2. 

The output of the program consists of a set of plots. The 
plots consist of the wake radius, the wind speed deficit at the 
center of the wake, the wind speed at the center of the wake, 
and the turbine power factor (ratio of the power which would 
be generated by an identical turbine in the wake of the first 
turbine to that power which would be generated if the downwind 
turbine were in the free stream wind) plotted as functions of 
the downv;ind coordinate. At the user's option, the geometric 
parameters may be plotted in physical units or as normalized 
by the rotor radius. Plots of the wind speed in the wake as a 
function of the lateral coordinate for eight altitudes for each 
of the selected downwind locations are also made. At the option 
of the user, the plots are normalized by the free stream wind 
speed at the hub altitude or by the free stream wind speed at 
the altitude at which the plots are made. In addition, a plot 
showing the vertical wind speed profile at each of the downwind 
locations of the lateral wind speed profiles is made. 

The parameter, Og, is an indication of the atmospheric 
turbulence. It may be entered in one of two ways. First, it 
may be entered as the standard deviation of wind direction as 
measured by tower-mounted anemometers. Alternatively, the 
atmospheric turbulence may be represented by the Pasquill 
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atmospheric stability class. The stability class varies from 
class A (highly unstable atmosphere) to class F (highly 
stable atmosphere) . Table 3-3 (taken from ref. 11) provides 
a key to stability classes. 

"Strong" incoming solar radiation corresponds to a solar 
altitude of greater than 60° with clear skies; "slight" inso- 
lation corresponds to a solar altitude from 15° to 35° with 
clear skies. Cloudiness will decrease incoming solar radiation 
and should be considered along with solar altitude in deter- 
mining solar radiation. Incoming radiation that would be strong 
with clear skies can be expected to be reduced to moderate with 
broken (5/8 to 7/8 cloud cover) middle clouds and to slight 
with broken low clouds. These methods will give representative 
indications of stability over open country or rural areas, but 
are less reliable for urban areas. This difference is due 
primarily to the influence of the city’s larger surface 
roughness and heat island effects upon the stability regime over 
urban areas . 

In general, inter turbine spaceing for wind turbines should 
be based on Pasquill stability class D, which corresponds to 
neutral atmospheric stability. From Table 3-3, stability classes 
other than class C and class D occur only when wind speeds are 
too low for turbine operation near rated power. Since wake 
recovery will be faster for class C and for class D (because of 
greater ambient turbulence for class C) , class D is the 
appropriate design condition. 

The value of the initial veloctiy ratio, m, to be used 
as input may be determined from one-dimensional momentum theory 
in the following manner taken from Ref. 12. Let be the 
free stream wind speed, be the wind speed through the turbine 
disk, and Uq be the initial wind speed of the wake (i.e., after 
expansion due to potential effects as discussed in the 
following discussion) . For the one-dimensional momentum 
analysis, power is extracted from the disk uniformly over the 
disk area, A. The axial thrust on the disk is 

T = Momentum flux in - Momentum flux out 


or 


T = m(U^-UQ) = pAU^(U^-Uo) (3-1) 

where m is the mass flow rate of air passing through the 
turbine disk, and p is the mass density of the air. Also, 
from pressure considerations, the thrust can be expressed as 
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T = ACp'^-p') (3-2) 

where p^ and p are the static pressures on the upwind and down- 
wind sides of the turbine, respectively. Applying the momentum 
equation (Bernoulli's equation) on the upwind side of the 
turbine gives 


= p'^+^pU^^ (3.33 

where p^ is atmospheric static pressure. On the downwind side 
of the turbine 


V~+hp\ij^ = p^+%pUp^ (3-4) 

Subtracting these equations to get (p^-p ) and using equation 
(3-2) gives 


T = A(p^-p') = J2PA(U^2 -Uq2) (3-5) 

Equating this with equation (3-1) gives 


U.J, = %(U^+Uq) (3-6) 

From equation (3-6), it is seen that the wind speed through the 
turbine is the average of the free stream wind speed ahead of 
the turbine and the wind speed in the expanded wake of the 
turbine . 

The axial induction factor, a, is defined by 

E U^(l-a) (3-7) 

Therefore, from equation (3-6), 

Uq = U^(l-2a) (3-8) 

and the initial velocity ratio, m, is 


m 


l-2a 


(3-9) 
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The initial velocity ratio can also be expressed in terms of 
the turbine output power. Because power is given by the mass 
flow rate times the change in kinetic energy, AKE, of the wind 
flowing through the turbine, the power, P^, is 


= mCAKE) 



(3-10) 


or, substituting for U.j, from equation (3-7) and for Uq from 
equation (3-8) and simplifying, 

P^ = 2pAU^'a(l-a)2 (3-11) 

The power, P , has a maximum at a = 1/3 and a minimum at a = 1. 
Therefore, i? the turbine is operating at its maximum power for 
the given free stream wind speed, m = 3. If the turbine is not 
extracting the maximum power available at its free stream wind 
speed, the appropriate value of "a” must be calculated from 
equation (3-11). An iterative solution is necessary because 
equation (3-11) is a cubic equation in "a". The power to be 
used in equation (3-11) is the power extracted from the wind, 
not the electric power output of the generator. The aero- 
dynamic power is the power extracted from the air by the tur- 
bine blades. For an operating tubine, it may be obtained from 
the output power and the generator and shafting efficiencies. 

After the values of the input parameters have been read, 
the program writes the values of the input parameters . 


3.3 Calculation of Program Constants 

Geometric constants for the program are calculated first. 

The constants include the wake growth rate due to ambient 
turbulence, the initial wake radius, the downwind extent of 
Region I, the wake radius at the end of Region I, the downwind 
extent of Region II, and the wake radius at the end of Region II. 

Recalculation of input geometric parameters . - The program 
allows input of geometric parameters in physical units or in 
rotor radii. If they have been input in physical units, they 
are corrected to units of rotor radii. 


(3-12) 
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^0 


(3-13) 

Az 

= iZ/Rj 

(3-14) 

Ay 


(3-15) 


Calculation of the initial wake radius . - From the 
momentum analysis presented above, the wake expands from the 
wind speed through the physical extraction disk, U-j,, to the 
initial wake wind speed, Uq . Let be the radius ^o£ the 
turbine disk, and let Rq be the radius of the wake after 
expansion to speed, Uq . The mass flow rate of air through 
the disk is 


m = pttR^^U^ = pttRq^Uq 
F rom equations (.3-6) and (3-9) 


' Uq ■ 2Ug ■ 2 


(3-16) 


(3-17) 


Therefore, the initial wake radius is given in terms of the 
disk radius as 


r 


0 


R. 


R 


/(m+l)/2 


(3-18) 


For the calculation of wake parameters, as indicated 
previously, the results of Abramovich were derived for 
mechanically- generated turbulence only (no ambient turbulence). 

In the following sections, the results of Abramovich are modified 
to include the effects of ambient turbulence. 

Calculation of wake growth due to ambient turbulence . - 
The program next calculates wake growth due to ambient turbu- 
lence. The theory for the turbulent dispersion of plumes 
(e.g., pollutants) in the atmosphere is used as the basis 
for the calculation of wake growth due to ambient turbulence. 
There are three steps in the following discussion. The first 
step is the determination of the plume growth as measured by 
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the concentration o£ atmospheric pollutants. The second step 
relates the concentration o£ atmospheric pollutants to the wind 
speed de£icit in the wake o£ turbines. The third step relates 
the pro£ile parameters £or pro£iles used £or pollution concen- 
tration work to the pro£ile parameters £or pro£iles used £or 
wakes . 

In atmospheric dispersion work, the Gaussian distribution 
is used to describe the concentration o£ pollutants. This 
distribution is given in axisyminetric £orm as 


Xc'Xco ^Xc 


(3-19) 


where x 



r 

a 


concentration o£ pollutants 

concentration o£ pollutants at plume center 

free stream concentration of pollutants (usually zero) 

radial coordinate 

pollution dispersion coefficient. 


As mentioned previously, the atmospheric turbulence may be 
input as the Pasquill stability class or as the standard devia- 
tion of wind direction. Input as the Pasquill stability class 
is considered first. Figure 3-2 of Ref. 11 gives a as a function 
of X, the distance downwind of the pollution source, for the six 
Pasquill stability classes. That figure is reproduced here and 
shown as Figure 3-3. From this figure, the following values 
for da/dx v/ere calculated. 


DERIVATIVE OF POLLUTION DISPERSION COEFFICIENT 
FOR PASQUILL STABILITY CLASSES 


Pasquill 

stability class do/dx 

A 0.212 

B 0.156 

C 0.104 

D 0.069 

E 0.050 

F 0.034' 
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Reference 13 relates the Pasquill stability class to Og, 
the standard deviation of- wind direction of the atmosphere. 
This relation is given as follows. 


STANDARD DEVIATION OF WIND DIRECTION FOR 
PASQUILL STABILITY CLASSES 


Pasquill 
stability class 

A 

B 

C 

D 

E 

F 


Oq 

Cdegrees) 

25.0 

20.0 

15.0 

10.0 
5.0 
2.5 


A least squares curve fit of the data given in the two pre- 
ceding tables gives 


n Rrr 

— = O.O31e'^®'^0 C3-20) 

dx 

Since a = 0 at x = 0 , 

aU) = O.OSlxe'^^'^O C3-21) 

If the atmospheric turbulence is input as a Pasquill 
stability class, table pg. 13 is used to generate a value for 
do/dx. If the atmospheric turbulence is input as a value 
of Og, equation C3-20) is used to generate a value for do/dx. 

As the second step, it is desired to relate the distribution 
of concentration of pollution to the distribution of wind speed 
deficits. In the discussion proceeding equation C5.25), 
Abramovich states that according to Prandtl's old and new 
assumptions of free turbulence, the dimensionless profiles of 
temperature and wind speed are the same, but according to 
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Taylor's theory of free turbulence, they differ. In the 
discussion proceeding equation CS-27), Abramovich states that 
the mechanism of lateral transfer of heat and of admixture 
is the same; consequently, the profiles of concentration 
difference must be similar to the profiles of temperature 
difference. This is stated in equation C5.27) of Abramovich. 
Combining this result with equation (5.25) of Abramovich 
gives 


Ax 


AXc 



u-u 


-u 

c 


CO 


00 


(3-22) 


where U = wind speed in the wake 

= wind speed at the center of the wake 
= free stream wind speed. 

Combining equation (3-22), the result of the Taylor theory 
mentioned above, with equation (3-19) for a Gaussian wind speed 
profile 


Au 


Au, 


-2rV2o2 

= e 


(3-23) 


As a third step, it is desired to relate the a of the 
Gaussian wind speed profile given by equation (3-23) to the 
wake radius, r 2 > for the wind speed profile used in the model 
presented herein. In the far wake, the wind speed profile is 
given by equation (5.23) of Abramovich as 


Au 


Au 


c 



(3-24) 


The wind speed profiles of equation (3-23) and (3-24) are 
very similar and are compared in Figure 3-4. In determining 
the relationship between o and Au is the same for both 
profiles, and the mass deficit is the^same for both profiles. 
The equality of mass deficits is 

CO ^ 


27tAu 

c 



0 


2ry2a2^ 


dr 


2irAu^ 



2 

r dr 


(3-25) 
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Dividing by 27rAu and integrating between the stated limits 
gives 


a" 9 

— = — r-2 (3-26) 

2 70 ^ 


or 


= 1.97a (3-27) 

Therefore, the wake growth rate due to ambient turbulence is 

a = Cdr 2 /dx)g^ = 1.97(da/dx) ( 3-28 ) 

where (da/dx) is obtained from the table on p. 13 or from equation 
(3-20) according to the method of input of atmospheric turbulence. 

The above derivation is based on the Taylor theory of free 
turbulence. This is the approach used by Abramovich, and 
Abramovich presents a curve (Figure 5.10 in Abramovich) which 
shows excellent agreement between experimental data and the 
Taylor theory. Therefore, the Taylor theory was accepted for 
use in this model. The Prandtl assumption is stated as 


Ay Au 


AXc Au^ 


(3-29) 


which is analogous to equation (3-22) for the Taylor theory. 
Under the Prandtl assumption, equation (3-23) would be 


Au 


Au, 




(3-30) 


Tracing through the derivation above gives the result 


x^ = 2.79a (3-31) 

Therefore, the Prandtl assumption can be used in the model by 
replacing the 1.97 in equation (3-28) with 2.79. 
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Calculation of wake characteristics at the end of Region I . 
Let be the free stream wind speed, U be the initial wind 
speed in the wake (c.f. equation (3-9)), and U be the local 
wind speed in the wake. Abramovich assumes the wind speed 
profile in Region I to be 


Uq-u 


= £(n) 


d-n'*')' 


where, in the notation of Figure 3-1, 


(3-32) 


n 


R^-R t2-r 


^ 2'^1 ^ 2'^1 


(3-33) 


Equation (3-32) is based on experimental data. 

By definition, the end of Region I is that point at 
which the potential core vanishes. Thus, the wind speed at 
the center of the wake at the end of Region I is U^. A 
momentum balance between the initial wake and the end of 
Region I gives 


^21 ^ ^ ^ ^0 

R^ Rq R^ /O. 214+0. 144m 


(3-34) 


where ^ 21/^0 given by equation (5.19) of Abramovich, and rg 
is given'^by^equation (3-18). 

From the equation for boundary layer growth given by 
his equation (5.1), Abramovich gives the length of the initial 
region of the wake as (equation (5.20)) 


rgCl^m) 

Rg R^ 0.27(m-l)/0. 214+0. 144m 


(3-35) 


The m subscript denotes that the quantity is associated with 
mechanically-generated turbulence (and is given by the Abram- 
ovich model) . 
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Equation (3-34), which gives the wake radius when the 
potential core has been completely eroded, is derived from 
conservation of the mom.entum deficit from the initial wake, 
the fact that U = Uq at the center of the wake at the end 
of Region I, and the assumption of the wind speed profile 
given by equation (3-32). It is therefore valid regardless 
of the presence or absence of ambient turbulence. Therefore, 
the presence or absence of ambient turbulence only affects the 
distance, Xj^, at which the end of Region I occurs. 

In the development of the model, seven approaches for the 
downwind extent of Region I were considered. The approaches 
are described and compared in Appendix D. Two of the approaches 
are slightly different in concept, but yield identical results. 
They are physically more justifiable than the other approaches 
and give numerical results that lie near the middle of the 
numerical results of all of the other approaches. Therefore, 
these two approaches have been selected for the model. The 
reader is referred to Appendix D for a description of the other 
five approaches . 

Figure 3-5(a) shows Region I for the Abramovich solution. 

The three areas shown are the free stream flow, the potential 
core, and the boundary layer between the free stream flow and 
the potential core. Also shown is the line which passes 
through the initial wake boundary and the midpoint of the wake 
radius at the end of Region I. The effect of ambient turbulence 
is shown in figure 3- 5(b). 

The first approach is based upon the growth of the boundary 
layer, b. Under this approach, the downwind extent of Region I 
is defined as that point at which the width of the boundary 
layer is T 2 \’ Since the wake radius ia also X 2 i at the end 
of Region I and 

b = r2 - rj^ (3-36) 

the radius of the potential core, r^ , must be zero at this point. 
For mechanical turbulence, the growth rate of the boundary layer 
is 



(3-37) 


Furthermore, for this formulation, the ambient turbulence exists 
in both the free stream flow and in the potential core. Since 
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the ambient turbulence affects both sides of the boundary layer, 
and a is the wake growth rate due to ambient turbulence on one 
side of the boundary layer (since it was developed as the rate 
of growth of the wake radius) , the total growth rate of the 
boundary layer is 


db 

dx 



(3-38) 


where the Pythagorean sum of the mechanical turbulence and the 
ambient turbulence has been used. Since the wake radius is r-. 
at the end of Region I, the downwind extent of Region I is 


X 


H 


db/ dx 



(3-39) 


An alternate approach to the downwind extent of Region I 
results from the assumption that the boundary layer develops 
about its own center. For this assumption, the erosion of 
the inner core, or the growth of the outer radius, is measured 
relative to the line passing from the initial wake radius to 
half the wake radius at the end of Region I, as shown in fig- 
ure 3-5(b). For the erosion of the inner core, let b^ be the 
distance from the edge of the potential core to the midpoint 
of the boundary layer. Hence, for mechanical turbulence. 


d ^\ ^ 0 • ^^21 

dx /m 


(3-40) 


Adding the ambient turbulence from the inner core by the 
square root of the sum of the squares gives 



Since b^ = 0.5r22 the end of Region I , 


(3-41) 
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/ 0-5t'2i V 

/ 


(3-42) 


+a‘ 


It is noted that multiplying both the numerator and the 
denominator o£ equation (3-42) by 2 gives equation (3-39) . 


The same result is obtained when the growth of the outer 
radius of the wake is considered. Let h 2 be the distance 
from the mid point of the boundary layer'^to the outer radius of 
the wake. Then at the end of Region I 



m 


0 . 5r 


21 




Adding ambient turbulence gives 



Since b^ 


O.Sr^^ the end of Region I, 


(3-43) 


(3-44) 



(3-45) 


which is identical to equation (3-42) . 

The parameter for the wind speed deficit at the center of 
the wake is 


Au_ = 


AU U -U 
c _ °° c 

U -U. 
Oc 00 0 


(3-46) 
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which is the wind speed dificit at the center o£ the wake 
divided by the initial wind speed deficit. The c subscript 
denotes the wind speed in the center of the wake. Because 
the wind speed in in the core of the initial region, 

Au^ = 1 in Region I; 

Downwind extent and radius at end of Region II . - Region 
II is a transition region from the wind speed profile form 
of the initial wake to the wind speed profile form of the 
far wake. The form of the wind speed profile in Region I 
is given by equation C3-32] . Using this form for the wind 
speed profile and the condition that the wind speed at the 
center of the wake is Uq , the wake radius at the end of Region 

I (denoted by is given by equation (3-34) and is derived 

from the fact that the momentum deficit at the end of Region I 
must equal the momentum deficit of the initial wake. As shown 
later in the discussion of Region III, in Region III, the 
wind speed profile is of the form given by equation (3-53) . 
Using the same conditions (i.e., wind speed at the center of 
the wake is Uq and the momentum deficit of the wake must equal 
the initial momentum deficit of the wake), the wake radius 
calculated is greater than t^ 21 ‘ this wake radius be the 

wake radius at the end of Region II and be denoted by ’^ 22 " 

The downwind extent of Region II is the distance required to 
allow wake growth from r 2 \ to X 22 ^t the wake growth rate 
of Region I. The downwind extent of Region II is given by 

Xj^ = nxj^ = nXj^/R^ (3^47) 

where n is taken from equation (5.124) in Abramovich as 

»/0. 214 + 0. 144m 1-/0.134 + 0 . l24m 

n = 1 3-48 1 

1-/0. 214+0. 144m /O . 134+0 . 124m 

The wake growth rate in Region II is identical with that 
in Region I. Therefore, the wake radius at the end of Region 

II is given by (cf . , Figure 3-1) 


(^71 


"21 '^ 0 - 


^ 22'^0 


( 3 - 49 ) 
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or, dividing by gives 


^22 “ 


The relationship of 

from, the matheTnatical forms of the wind speed profiles in 
Region I and Region III and the assumption that the 
wake growth rate in Region II equals that in Region I. 

If these same assumptions are made in the presence of 
ambient turbulence the presence of ambient turbulence 

does not change the form of the wind speed profile in 
Region I or in Region III, and the wake growth rate 
in Region II equals that in Region I), then the relation- 
ship of equation (3-48) is valid in the presence of ambient 
turbulence . 


^ = Tg+nCr^^-rd] 

equation (3-48) is derived solely 


3.4 Wake Calculations in Region III 


Wake growth in Region HI - - Region III of the wake is 
the region in which the mecHanically-generated turbulence decays. 
Thus, at the beginning of the region, wake growth is governed 
by both mechanically- generated turbulence and ambient turbulence. 
The wake growth due to mechanically-generated turbulence asymp- 
totically approaches zero as the downwind coordinate, x, 
increases. In Region III, the wind speed profiles are self- 
similar, that is, they have the same mathematical form at all 
downwind locations. If the applicable expressions from 
Abramovich are used, the wake growth must be calculated by 
numerical integration. 

Let it be assumed that the wake radius, r , is known 
at distance, x, downwind of the turbine, where T 2 and x are 
values of the wake radius and downwind distance which have 
been normalized by the rotor radius, R^. Let T 2 ' = r 2 +Ar 
be a wake radius which is slightly larger than r 2 . It is desired 
to find the downwind distance at which the wake radius is ^ 2 ' • 


For the main region of the jet, the unnumbered equation 
proceeding equation (5.97) of Abramovich gives the wake radius 
as a function of the wind speed deficit at the center of the 
wake as 


^2 = ^0 


no -mni 

2u lu 


.^l-m)'(A2A^c'^A^(l?iVu^)J 


(3-51) 
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where Au^ is defined by equation (3-46) and 


"lu = "2u " 1 


(3-52) 


in the case of uniform fields of velocity and density at the 
initial cross section of the jet (equation following equation 
(5.39) of Abramovich). 

In the main region of the jet, Abramovich assumes (based 
on experimental data) that the form of the wind speed profile 
in the far wake is 


AU U -U 

CO 

AU U -U 

C 00 c 



(3-53) 


where U is the wind speed at the center of the wake. From 
this assumed form of the wind speed profile, equation (5.86) 
of Abramovich gives 


A^ = 0.258 


(3-54) 


and 


= 0.134 

Therefore, equation (3-51) is 


^2 = 


/0.258m\ 

0 .134Au2-t-( jAu^ 

\ m-1 / 


(3-55) 


(3-56) 


Equation (3-56) is a momentum equation which equates the momentum 
deficit in the main region of the wake to the initial momentum 
deficit of the wake. Rearranging equation (3-56) to solve for 
Au gives 


0 . 134Au^ 


'0 . 258m\ 


-lAu + 


m 


-1 / (m-1) 


0 


(3-57) 


or 


Au 


c 


3.73 


(0.258m 

yo. 258mY 0.536rQ^' 


( m-1 

_\ m-1 / (m-1) . 

) 


(3-58) 
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where the negative sign for the quadratic equation was chosen 
so that Au^ goes to zero as T 2 becomes large. 

From his equation (5.28), Abramovich assumes that the wake 
growth rate is directly proportional to the difference between 
the wind speed of the free stream flow and the wind speed at 
the wake center and is inversely proportional to the mean wind 
speed of the wake. Using the first part of equation (5.31) of 
Abramovich gives 


+ c- 


dx 

dr. 


2U 


2U 


= 1 +- 


= 1 - 


U^-U 

C 00 


Au (U -U„) 

C oo 0 


(3-59) 


where the definition of Au given by equation (3-46) has been 
used in the last term of the equation. Dividing by U^, using 
the definition of m given by equation (3-9) and rearranging gives 
the wake growth rate due to mechanical turbulence (using 
c = -0.27 for m > 1) as 


'dr. 


i.dx /m 


.27 - 


2m 


(m-1) Au 



(3-60) 


The calculation procedure for wake growth in Region III 
is as follows. Let it be assumed that the wake radius, r 2 > is 
known at distance, x, downwind of the turbine. Let r 2 * ='^r 2 +Ar 
be a wake radius which is slightly larger than r 2 * Let 
(dr„/dx) be the wake growth rate due to mechanically-generated 
turoul e^:e calculated from equation (3-60) for wake radius, r 2 - 
The value of Au^ used in equation (3-60) is calculated from 
equation (3-58). Let (dr„/dx) * be the wake growth rate due to 
mechanically-generated turbuli^ce fPr wake radius, 1 ^ 2 '’ 
before, the wake growth rate due to ambient turbulence is 
taken as a. Therefore, the total wake growth in the interval 
Ar between r 2 and r 2 * is given by 



(drz/dx).^ 



( 3 - 61 ) 


The downwind distance, x' , at which the wake radius is r 2 * is 
given by 


Ar 

x' = x+ 


(dr2/dx)g 


(3-62) 
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Equations C3~60) through (3-62) are solved iteratively, and 
r-, is thereby generated as a function of x. In the calculation 
procedure, the calculations of equations (3-58) and (3-60) are 
performed in Subroutine R3 since the calculations must be 
repeated many times for different values of V 2 and ' 

If one of the downwind locations at which wind speed profiles 
are to be calculated is encountered during the numerical inte- 
gration of Region III, the wake radius at that point is retained. 
The wake radius at the desired value of x is obtained by linear 
interpolation between the values of x which bracket the desired 
value of X. The numerical integration is continued downwind 
until X > 50 or r 2 > 7, whichever occurs first. 

Turbine power factor . - As mentioned previously. Subroutine 
R3 is used to perform the calculations of equations (3-58) and 
(3-60). The subroutine also calculates the turbine power factor 
which is defined as 


P = P/P 

' CO 


(3-63) 


where P = power which is generated by a turbine which 

is geometrically identical to another turbine 
and centered in the wake of the other turbine 

P^ = power generated by the turbine in the free 
stream wind. 

The geometry related to a wind turbine in the wake of 
another wind turbine is shown in Figure 3-6 . This geometry 
represents a worst case condition because the downwind turbine 
is centered in the wake of the upwind turbine. For actual 
turbines, the wake oscillates because of changes in the wind 
direction. However, in the present analysis, it is asssumed 
that the downwind turbine is fixed in the center of the wake 
of the upwind turbine. 

Immediately upwind of the turbine, the stream tube of the 
flow passing through the turbine expands from a radius, R_^, and 
wind speed U^, to the disk radius, R,, and the wind speed, U,p, 
through the turbine disk. It is assumed that the distance, 
over which the potential effects occur (cf. equations (3-1) P’ 
through (3-6)) is small compared with the distance between the 
turbines. Let A be the cross-sectional area of the stream 

00 

tube in the free stream. Let the turbine power for a turbine 
in the free stream air be 



3 


(3-64) 
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where, from comparison with equation (3-11), the constant, 
K, is 


K = 2pAa (l-a) 


(3-65) 


For a turbine in the wake of another turbine, equation (3-64) 
must be modified to show the variation of the "free stream" 
wind speed across the stream tube for the downwind turbine. 
For a turbine in the wake of another turbine, the power is 


P = 



(3-66) 


or, normalizing R by and normalizing U by gives 


P = 27tKU ^R/ 


/ [u“r)] 
0 


'r dr 


(3-67) 


Combining the definition of the turbine power factor from 
equation (3-63) with equations (3-64) and (3-66) gives 

r 


27tKU ^R,^ 
00 


/ CO 

[u(r)] 


r dr 


P = — = 
P 


ttKR ^ 

CO oo 


(3-68) 


or 


2 r 00 

P = — /[u(r)]^r dr 

oo 0 


(3-69) 


From the conservation of mass in the stream tube upwind of the 
turbine, the mass flow rate of air through the turbine is 

m = irpR^^U^j, = TTpRoo^U^ 


Using equation (3-6) for and rearranging gives 


(3-70) 
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U +U. 

R 2 = _^1_0 R 2 

” 2U 


(3-71) 


Dividing by Rj^U and using the definition of m given in 
equation C3-9j gives 



C3-72) 


The preceeding discussion has been related to the downwind 
turbine (i.e., the turbine which lies within the wake of the 
upwind turbine) . The wake of the upwind turbine forms the 
free stream wind flow for the downwind turbine. In the far wake 
(Region III), the form of the wind speed profile is given by 
equation (3-53). Combining this equation with the definition 
of Au^ given in equation (3-46) and the definition of m given 
in equation (3-9) and dividing by gives 


u = 



(3-73) 


If the center of the downwind turbine coincides with the 
center of the wake with the wind speed profile given by equation 
(3-73), the turbine power factor is 


P = 


2 /r / i\r /r 

/ h-Au^(l— ) l-( — 

"n t \ "'A \r?/ -I 


2 3 


r dr (3-74) 


For convenience in the derivation, let 


B = Au^ (1 — 


1 

m 


(3-75) 


Then 


P = 


2 /-“t r r A r 

/ <1-3B l-( — +3 bMi-( — -bMi-( — > 

I L \r2/ -I L Vr^/ L \r2/ -I ) 


r dr 
(3-76) 
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P = 


oo 0 


/> -3B+3B^-B^] + [6B-12B^+6B^]^ — ^ 


T 

+ [-3B+18B2-15B^]( — ) + [-12B2+20B^] 

^ 2 ' \^2 

' lV I 


'r 

r. 


. 4.5 


r x6 .7.5 


+ [3B^-15B^]( — +6bY — 


B^f — j [ r dr 
^2/ ^ 


(3-77) 


Integrating between the indicated limits and simplifying gives 

_ (1-3B+3B^-B^ 6B-12B^+6BVr - 3B+18B^ - 15B Vr„ 

P = 2 + — + 

(2 3.5 Xr^/ 

12B2 + 20BVr„\‘*-= 3B^ - 15B Vr.. 6B ‘ | 


— (3-78) 


6.5 \r. 


9.5\r. 


ll\rJ f 


The power factor F can then be expressed explicitly 
as a function of m and ri by evaluating B from equation. (3-58) 

J:? equation (3-18), and r from equation 

(3-72). The expression for B is 


B = 


3.731 


|o. 258- 1^0. 


0.268 


06656- 



Figure 3-7 shows calculated values of P as a function of t 2 for 
three values of m as obtained from equation (3-78) and the above 
expression for B. This plot is not generated during the computer 
run . 


This analysis was developed for a wake out of ground effect. 
That is, there is no effect of the ground plane and no wind 
shear because of the ground. Two other assumptions are involved. 
The first is that no significant wake growth occurs over the 
distance of the expansion of the stream tube from radius to 
the disk radius R, for the downwind turbine. The second assump- 
tion is that the axial induction factor, a, is the same for both 
turbines. Because the turbine power factor is based upon the 
form of the wind speed profile that exists in Region III, the 
turbine power factor is undefined in Regions I and II. 
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3.5 Wake Plots 

As mentioned previously, the numerical integration in 
Region III is terminated when x > 50 or r 2 > 7, whichever occurs 
first. Four plots are then made. The first plot shows the wake 
radius, the wake width, and the height of the top of the wake 
above the ground as a function of the downwind coordinate, x. 

The wake width is 


w = 2r2 C 5-79 ) 

The height of the top of the wake above. the ground is 


^ " ^2'^^a r 3-80 1 

The second plot shows Au^, as a function of x. The third plot 
shows the normalized wind speed at the center of the wake, u^. 
This is obtained from rearranging equation (3-46) to give 

Uc = (3-81) 


and then dividing by to give 


U, 

— = u^ = 1-Au^ ( 1 — 

U ^ \ my 

00 ' ' 


(3-82) 


The fourth plot shows the turbine power factor as a function 
of the downwind coordinate, x, in Region III. 

3.6 Calculation of Wind Speed Profiles 

After the plots described above have been completed, plots 
of wind speed profiles are made. The first set of wind speed 
profiles shows the wind speed as a function of the lateral 
coordinate for eight altitudes specified by the input values 
of Zq and Az. Plots are made for each of the input desired 
downwind locations. A plot of the vertical wind speed profile 
at the wake centerline is then made. A single plot shows the 
vertical wind speed profile at all of the downwind locations 
at which lateral wind speed profiles were made. 

Coordinates ♦ - If the downwind locations at which wind 
speed profiles are desired have been input in physical units, 
they are converted to units of rotor radii by 
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Xdj = X^j/Rd (3-83) 

where the notation Xj- denotes the downwind coordinate o£ the 
j th location at whicn^wind speed profiles are desired. 

During the calculation of the wake radius as previously 
described, the wake radius corresponding to each x , . is 
calculated and retained. Let r . be the wake radiui at x = x, • . 
For the x^^ which are less than^Xj^ 



ro-^(r22-ro)Xd^/x^ 


(3-84) 


For the x^- which are larger than Xj., the values of r - are 
retained as the values of the x,. aYe determined duriAg the 
numerical integration of Region^lII. The wake radius at x , - 
is obtained by linear interpolation between the values of x^^ 
which were calculated in the numerical integration and which 
bracket x, . . 


The generation of the lateral wind speed profiles begins 
at z = zg and y = 0 (axis of turbine) . Values of the wind 
speed in the wake are calculated at increments. Ay, beginning 
at y = 0 . Once the value of y is in the free stream (indicated 
by u = 1) , the free stream value of the wind speed is extended 
across the plot to y = 6 rotor radii. The altitude is then 
incremented by Az, and the wind speed profile at the next alti- 
tude is generated. The process is repeated until wind speed 
profiles have been generated for eight altitudes. Since the 
lateral profiles are symmetrical about y = 0, the wind speed 
for only positive values of y are determined. 


After the lateral wind speed profiles have been generated 
for a given downwind location, , the vertical wind speed 
profile is generated. The vertical wind speed profile is cal- 
culated at y = 0 only. The profile begins at z = 0 with in- 
crements of the input value of Ay used for successive calcula- 
tions (i.e., Az = Ay for this calculation). 


The altitude of the point on the wake profile (y,z) rela- 
tive to the axis of the turbine, denoted by Zy, is also calcu- 
lated. It is 


z 


V 


z-ha 


where z is the altitude as measured from the ground. 


(3-85) 
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Subroutine CALCU is used to calculate the wind speed in 
the wake for given values of y and z . The subroutine first 
calculates the radius of the point (y,z^) as 

r = /y^+z^^ ( 3-86 ) 

Region I . - If x 
of the potential core 




< X 


H’ ^dj 


is in Region I. The radius 


r^ = rgCl-x^j/x^) ( 3-87 ) 

Consider n as defined by equation (3-33) where T 2 = • . If 

n > 1, r is inside the potential core and 

u = 1/m (3-88) 


If n < 0, r is outside the wake boundary, and 

u = 1 (3-89) 


If 0 < n < 1 , r is in the boundary layer. Rearranging equation 
(3-32), dividing by U^, and using the definition of m from equation 
(3-9) gives ” 


u = 




(3-90) 


After the value of u is calculated, control is returned to the 
main program. 


Region III . 


If X 


dj 


> X 


N> 


X 


dj 


is in Region III. 


If r > r 2 . 


u = 1 


(3-91) 


If r<To, Subroutine R3 is called to calculate Au^. from equation 
(3-58).^ Then, with (U^-U) from equation (3-53) and (U^-U^,) from 
equation (3-46) , 


U -U = Au (U -U.) 
OO c o° 0 


1-1 


1 . 5 - 


2 


(3-92) 
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Dividing by and using the definition of m in equation (3-9) 
gives 


For this equation 


which agrees with 


u = 1-Au 1 — 

‘'V m/ 


H- 


1 . 5 


it is seen that if r = r^, u = 1. If r 


u = 1-Au^(l-- 

^ V m/ 


the definition of Au^ given in equation 


(5-95) 


= 0 , 

(3-94) 

(3-46) . 


Region II . - If . is in the transition 

region. The wind speed^is linearly interpolated between the 
profile at the end of Region I and the profile at the beginning 
of Region III. The actual radius of the wake in Region II is 
used. Let u-, be the value of u calculated from equations (3-33) 
and (3-90) using the value of t / t 2 in Region II. Let be the 

value of u calculated from equation (3-93) using the same value 
of r/r«. At the beginning of Region III, Au = 1. Then, for 
Region'^ II , 


u 


-X n 

H' ' N / 


(5-95) 


After calculation of u, control is returned to the main 
program. 


3.7 Calculation with Ground Effect 

The above equations describe the wind speed in the wake 
for an isolated turbine. The effect of the ground is to 
shield the lower part of the wake from the effect of the ambient 
wind which would otherwise act to accelerate the flow in the 
wake. The effect of the ground thus retards the acceleration 
of the wake flow by the surrounding free stream. Thus, in the 
presence of the ground effect, the wind speed in the wake is 
less than it would be if the ground were not present. 
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The presence of the ground is modeled by placing an image 
turbine at distance, h^, below the ground. The imaging technique 
is shown in Figure 3-8. Let denote the altitude of the 

point Cy,z), relative to the axis of the image turbine. Then 

z^* = z+h^ ( 3-96 ) 

Using the value of z^* instead of z , Subroutine CALCU is called 
to calculate u* , the^normalized wind speed at point (y,z) in the 
wake of the image turbine. The total wind speed deficit in the 
wake in ground effect is the sum of the wind speed deficits in 
the wakes of the real turbine and of the image turbine. If u 
is the normalized wind speed for a wake in ground effect, then^the 
sum of the wind speed deficits is 

1-Ug = (l-u)+(l-u*) (3-97) 

where the 1 is the normalized wind speed of the free stream. 
Rearranging equation (3-98) gives 


u = u+u*-l ( 3-98) 

o 

By definition, the approach conserves the total mass deficit. 

The mass deficit in the real wake is (1-u) . The mass deficit 
in the wake of the image turbine is (1-u*), and the mass deficit 
of the wake in ground effect is (1-u ). Equation (3-97) shows 
that the mass deficit for the wake iS ground effect is the sum of 
the mass deficits for the wakes of the real and image turbines. 

It is noted that the presence of the ground does not affect 
the shape of the wake boundary above the ground. This is shown in 
Figure 3-9. The only portion of the wake which is affected by 
the ground effect is that portion of the wake which lies in the 
intersection of the wake of the real turbine and the wake of 
the image turbine. 


3.8 Calculation of Wind Speed 

For the horizontal wind speed profiles, the input parameter, 
NP, specifies whether the wind speed is to be normalized by the 
free stream wind speed at the altitude of the profile or is to 
be normalized by the free stream wind speed at the hub altitude. 
If NP = 0, the wind speed in the wake is normalized by the free 
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stream wind speed at the altitude of the wind speed profile. The 
parameter, u , given by equation (3-98) has been developed as the 
normalized w?nd speed in the wake in a uniform free stream wind 
speed. Therefore, u = 1 in the free stream. For a nonuniform 
free stream wind speed profile, if the wind speed in the wake is 
normalized by the free stream wind speed at the altitude of the 
wind speed profile, the wind speed in the free stream is 1. 
Therefore, u is the wind speed in the wake normalized by the free 
stream wind speed at the altitude of the wind speed profile. 

In NP = 1, the wind speed in the wake is normalized by the free 
stream wind speed at the hub altitude. Then, the wind speed at 
the altitude of the wind speed profile is 

U = U u (3-99) 

where U (z/h^)^ is the free stream wind speed at the altitude of 
the wind speed profile, and (as indicated above) u is the wind 
speed in the wake normalized by the free stream wind speed at the 
altitude of the profile. 

If has been input in physical units, the graphical output 
of the program will be in physical units of wind speed instead 
of normalized by the free stream wind speed. In this case, the 
use of NP =1 will give output in physical units of wind speed 
for each altitude. The output for NP = 0 has no desirable 
interpretation if ^ 1. 

For the vertical wind speed profile, equation (3-99) is used, 
regardless of the value of NP or whether the output is normalized 
by the free stream wind speed or given in physical units. Thus, 
the wind speed profile of the free stream is always evident on the 
vertical wind speed profile is y > 0. 



Table 3-1 


FORIvIAT 

OF INPUT 

PARAMETERS 

FOR TURBINE WAKE 

COMPUTER 

PROGRAl^i 


Columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

Format 

Card 1 

ST 

CPM 

AH 

VHO 

WEXP 

RRR 

6F10.4 

Card 2 

J, 10, NP 
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Card 3 

XNPT(l) 

XNPT(2) 

XNPT(3) 

XNPTC4) 

XNPT(5) 

XNPT(6) 

6F10.4 

Card 4 

ZO 

DZZ 

DYY 




3F10 .4 




Table 3-2 


TABLE OF INPUT PARAMETERS FOR TURBINE WAKE COMPUTER PROGRAM 


Symbol 


Computer 

Symbol 


Definition 


a 


0 


ST 


Atmospheric turbulence parameter. If input as a 
positive number, it is the standard deviation of 
the wind direction. If input as a negative number, 
it is the Pasquill stability class as follows: 

-1 for Pasquill stability class A 
-2 for Pasquill stability class B 
-3 for Pasquill stability class C 
-4 for Pasquill stability class D 
-5 for Pasquill stability class E 
-6 for Pasquill stability class F 
If the input value is 0 . , the wake growth due to 
ambient turbulence is zero, and the Abramovich 
wake solution results. 


m 


CPM 


\ »a 


AH 


Ratio of free stream wind speed, U^, to the initial 
wind speed in the wake, Uq . 

Hub height of the turbine in rotor radii or in 
physical units as specified by 10. 


U 


00 


Y 


VHO 


WEXP 


Ambient wind speed at the hub altitude. If out- 
put is desired in physical units, input 
value should be in physical units. If out- 
put normalized by free stream value is de- 
sired, input should be 1.0. 

Coefficient of the power law profile for the 
free stream wind speed. 


ON 




Table 3-2 


TABLE OF INPUT PARAMETERS FOR TURBINE WAKE COMPUTER PROGRAM (Continued) 


Computer 

Symbol 


Definition 

Rotor radius. An input value of 1.0 will give 
all output of units of length in rotor radii. 

An input value of 0.5 will give all output of 
units of length in rotor diameters. An input 
value greater than 2.0 will give all output of 

units of length in the physical units used for 

the rotor radius. 

Number of downwind locations at which wind speed 
profiles are to be calculated. 

Input option for parameters with physical dimensions 
of length (AH, XNPT, ZO , DZZ, DYY) 

10 = 1 for input in rotor radii 

10 = 2 for input in physical units (must be same 

physical units as used for rotor radius) , 


NP NP Specifies how velocity is to be normalized 

for plots of wind speed profiles in the lateral 
direction . 

NP = 0 normalizes wind speed by the free stream 
wind speed at that altitude 
NP = 1 normalizes wind speed by the free stream 
wind speed at the hub altitude. 




Xd or 


XNPT(j) 


Downwind locations at which wind speed profiles 
are to be calculated (rotor radii or physical units 
as specified by 10). 



Table 3-2 


TABLE OF INPUT PARAMTERS FOR TURBINE WAKE COMPUTER PROGRAM (Concluded) 


Symbol 

Computer 

Symbol 

Definition 

Zq or Zq 

ZO 

Minimum altitude at which wind speed profiles 
are to be calculated (rotor radii or physical 
units as specified by 10) . 

Az or aZ 

DZZ 

Increment in altitudes at which wind speed 
profiles are to be calculated (rotor radii 
or physical units as specified by 10) . 

Ay or AY 

DYY 

Increment in y by which calculations are to 
by made in generating the wind speed profiles 
(rotor radii or physical units as specified 
by 10) . 




Surface wind 
speed at 10m 
(m/sec) 


Day 


Night 


Incoming 

Strong 

solar radiation 
Moderate Slight 

Thinly overcast 
or 

>4/8 Low cloud 

<3/8 Clouc 

<2 

A 

A-B 

B 



2-3 

A-B 

B 

C 

E 

F 

3-5 

B 

B-C 

C 

D 

E 

5-6 

C 

C-D 

D 

D 

D 

>6 

C 

D 

D 

D 

D 







Turbine 

disk 



Rj = Radius o£ turbine disk 
d 

Rq = Initial wake radius 

R^ = Radius o£ potential core 

R 2 = Outer radius o£ wake 

R-i= Outer radius o£ wake at end 
o£ Region I 


R-~= Outer radius o£ wake at end 
o£ Region II 

Xu = Downwind extent o£ Region I 
n 

Xj^ = Downwind extent o£ Region II 


Figure 3-1. - Wake geometry £or the wake computer model. 
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Figure 3-2. Flow diagram for wake model computer program. 
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Figure 3-2. Flow diagram for wake model computer program 
(continued) . 














3 


i 


r N. 

^ Close plot files / 



Figure 3-2. 


Flow diagram for wake model computer program 
(concluded) . 




Pollution dispersion coefficient. 
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Figure 3-3. Pollution dispersion coefficient as a function 
of distance downwind of the source. 




Wake radius parameter, r/r 



Figure 3-4. Comparison of the Gaussian wind speed profile 
used in plume dispersion analysis with the 
wind speed profile used in the far wake of 
the wake model. 
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Distance downwind of turbine, x 


Ca) Region I for Abramovich solution 
(no ambient turbulence) 


!h 


V) 

•H 

cd 
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c3 


(b) Region I with ambient turbulence 


Figure 3-5. 


Geometry of Region I of the wake for 
calculating the downwind extent of 
Region I. 




Figure 3-6. - Geometry o£ one turbine centered in the 
wake of another identical turbine. 



Turbine power factor. 



Wake of wind turbine 
Wake of image turbine 



Figure 3-8. - Geometry for calculation of ground plane effect by 
image technique . 
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Figure 3-9. Illustration o£ the portion of the wake 
affected by the ground plane. 
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4.0 COMPARISON WITH PREVIOUS MODEL 


There are some differences between the model presented 
in the previous section and the AeroVironment model presented^ 
in reference 8. For the purpose of completeness, those dif- 
ferences are documented in this section. 


4.1 Wake Growth Rate Due to Ambient Turbulence 

For the calculation of wake growth due to ambient tur- 
bulence, the AeroVironment model uses the expression 

dR, 

dx - ITTsT ^ 

where a is a direct data input for the calculation. Methods 
for determining the value of a were suggested, but not de- 
tailed. These included evaluating a from the effective 
turbulence intensity of the atmosphere at the site, or from 
the growth rates of smoke or polutant plumes as given in 
atmospheric dispersion theory. The present model uses the 
growth rate due to ambient turbulence as 


dR2 
dX " 


(4-2) 


where a is calculated internally from input considerations 
of Pasquill stability class designation for the atmosphere as 
given in dispersion theory. This is done in conjunction with 
the Taylor assumption for the relationship between velocity 
deficit profile and admixture profile (corresponds to the 
1/0.51 in the AeroVironment formulation). 


The AeroVironment model has provision for different values 
for vertical and lateral (normal to the wind direction) com- 
ponents of wake growth due to ambient turbulence. Under these 
conditions, the wake (out of ground effect) assumes an ellip- 
tical cross section instead of a circular cross section. This 
condition was not included in the present version for sim- 
plicity and because of the uncertainty about the appropriate 
procedure for relating unequal wake growth rates to an axisym- 
metric flow solution. In the present model, the wake growth 
rate due to ambient turbulence was taken as the lateral growth 
rate. 
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4.2 Support Tower 

The effect of the wake of the turbine support tower was 
not included in the present version because it was believed 
to be an unwarranted complexity. Tower wake relationships de- 
pend on the type of construction (i.e., lattice or tubular) 
and the interaction between the tower wake and the turbine 
wake. Furtherfore, the effects of the tower wake should be 
diminished in the far wake of the turbine, which is the prin- 
cipal region of interest. Thus, the tower wake effect was 
deemed negligible compared to the other uncertainties in the 
model formulation. 


4.3 Calculations for Region I 

The calculation of the wake radius at the end of 

Region I was done differently in the revised model than in 
the AeroVironment model. The AeroVironment model used equation 
(5.21'0 of Abramovich. At the end of Region I, the boundary 
layer width, b, is i^21’ Abramovich's notation, y is the 
distance from the initial wake radius to the boundary of the 
potential core. Therefore, at the end of Region I, y^ = rQ, 
and equation (5.21') of Abramovich can be written as 


= 0. 416 + 0. 134m+0.021^ici+0.8m-0.45m2) 
^21 ^0 


(4-3) 


This is a quadratic equation in ^2l/^0 which is solved for 


to give 


-f+/fi+4g 


2g 


(4-4) 


where 


f = 0.416+0. 134m 


g = 0. 021(1+0. 8m-0.45m2) 


(4-5) 

(4-6) 


The positive sign was chosen for the solution to the quadratic 
equation so that i' 21 /^O positive. 
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Although this approach is correct, it was deemed to be 
unnecessarily cumbersome. The parameter, given 

directly by Abramovich, equation (5.19) as 


= /O. 214 + 0. 144m 

^21 


(4-7) 


This equation (same as eq. (3-34) was deemed to be simpler than 
that presented in the AeroVironment model. Numerically, 
the difference between these two methods is less than 0.41 for 
values of m between 1 and 3. 


In the AeroVironment model, only one approach for 
calculating the downwind extent of Region I was presented. 

This was the first approach (r^ approach) described in Appendix 
D. After consideration of this approach, it was realized that 
the second approach (r- approach) described in Appendix D was 
equally plausible from^a physical phenomenological point of view, 
but the to approach gave much different results. This paradox 
prompted the investigation of other approaches, and the seven 
approaches outlined in Appendix D resulted. The approach which 
was eventually chosen was not the approach originally pre- 
sented in the AeroVironment model. The reader is referred 
to Appendix D for the description of the AeroVironment ap- 
proach and the six other approaches considered. 


4.4 Wake Growth in Region III 

The AeroVironment model attempted to simplify the 
calculation of wake growth in Region III in order to circum- 
vent the necessity of the numerical integration in Region III. 

The basic idea by which the calculation procedure was to be 
simplified was to perform the numerical integration externally 
to the wake program for different sets of values of a and m. 

The numerical integration was performed over an interval of x 
of lOrg beginning with the wake radius ^ function 

of m only. Once the wake radius was known at the beginning 
and end of Region III, the effective ambient turbulence was 
defined as that' ambient turbulence which alone would give the same 
wake growth in Region III as the combination of mechanical 
turbulence and ambient turbulence had given in Region III by 
the numerical integration. For example, for m = 3 and a = 0.1, 

Tq = 1.414 (4-8) 

T 22 = 1.988 (4-9) 


5.488 


(4-10) 
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The downwind extent of Region III is lOrp or at 

X = Xj^+lOrQ = 19.63 (4-11) 

From the numerical integration of equations (3-57) through 
(3-62), ^2 - 3.43 at x = 18.42, and = 3.57 at x = 19.83. 

By linear interpolation, X 2 - 3.55 at'^x = 19.63. Therefore, 
in Region III 


/dr-\ 3.55-1.99 

a = ( — -j = = 0.1088 (4-12) 

\dx /avg 19.83-5.49 


Therefore, for m = 3 and a = 0.1, the effective growth rate 
in Region III is = 0.1088. 

In the AeroVironment model, the numerical ingetration 
was conducted external to the main wake computer program. 

Plots of Ug as a function of a were generated for several 
selected values of m, and was an input paramter for the 
wake computer program. In Region III, the wake growth rate was 
assumed to be constant at a value of a . In Region IV (downwind 
of Region III), the wake growth rate was assumed to be constant 
at a value of a. That is, downwind of x = Xj^+lOrQ, the wake 
growth was assumed to be due to ambient turbulence alone with 
no contribution to wake growth due to mechanical turbulence. 

In this discussion, the definition of a used in this report is 
the growth rate of the wake radius due to ambient turbulence. 

In the AeroVironment report, a was defined slightly differently, 
and the growth rate of the wake radius due to ambient turbulence 
is shown in the AeroVironment report as a/0.51. 

The method used in the AeroVironment report is reasonably 
accurate, acsept for low values of ambient turbulence. At low 
values of ambient turbulence (a<0.1), there is still a sig- 
nificant value of wake growth due to mechanical turbulence down- 
wind of X = Xj^+lOrQ. Also, it was deemed appropriate that the 
model should assume the exact form of the Abramovich solution 
for a = 0. The AeroVironment model does not do that since the 
wake growth rate in Region IV goes to zero for a = 0. For these 
reasons, in the revised model, the numerical integration of 
equations (3-57) through (3-62) was made an integral part of 
the model. 
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5.0 SAMPLE PLOTS 

This section contains sample plots made with the turbine 
wake computer program. The plots were made from seventeen 
runs of the program. Illustrated in this section are the in- 
put of parameters for the computer program and the output ob- 
tained. The first part of this section presents complete sets 
of plots from the computer program. These figures include 
plots of the Abramovich solutions (i.e., no ambient turbulence). 
The second part of this section shows the effect of changes 
in wake variables on wake recovery. Results are presented 
for changes in ambient turbulence factor; power law exponent, 
height of the axis of the turbine rotor, and initial veloc- 
ity ratio. 


5.1 Wake Plots 

Abramovich solutions . - Initially, it was desired to ob- 
tain plots for the Abramovich solution for the coflowing jet 
(a = 0). The input parameters , are shown below. 


VALUES OF INPUT PARAMETERS FOR TURBINE WAKE PLOTS 
FOR ABRAMOVICH SOLUTION 


ST = 0.0 

J = 5 

CPM = 1.5 

10 = 1 

AH = 50. 

NP = 0 

VHO = 1.0 

XNPT = 5., 

WEXP = 0. 

ZO = 48.5 

RRR = 1.0 

DZZ = 0.5 
DYY =0.05 


The input value, ST = 0.0, makes the ambient turbulence 
zero. The value of CPM = 1.5 was chosen as a typical value, 
for large wind turbines. The value AH = 50. was used to 
make the wake isolated from the ground effect. VHO = 1. was used 
to make all wind speed plots normalized by the free stream 
wind speed. The power law exponent, WEXP = 0., gives no variation 
in free stream wind speed with altitude. The rotor radius, 

RRR = 1.0 gives all geometric output in rotor radii. J = 5 
indicates five downwind locations at which wind speed profiles 
are to be calculated. 10 = 1 indicates that geomteric input 
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is in rotor radii. NP = 0 indicates that the wind speed 
profiles are to be normalized by the free stream wind speed 
at the altitude of the profile. The values of downwind 
location at which wind speed profiles are to be generated 
are defined by XNPT(j) = 5., 10., 20., 35., and 50.. For the 
wind speed profiles, ZO = 48.5, and DZZ = 0.5. The altitudes 
at which the wind speed profiles are to be calculated are 

Zi = Zo + iAz i = 0, . . ., 7 C5-1) 

Therefore, the altitudes of the lateral wind speed profiles are 
48.5, 49.0, 49.5, 50.0, 50.5, 51.0, 51.5, and 52.0. Since 
h^ = 50. for the input values of Zq and Az, the fourth wind 
speed profile is at the hub altitude. The second wind speed 
profile is at the lower edge of the turbine disk, and the 
sixth wind speed profile is at the upper edge of the turbine 
disk. The value, DYY = 0.05, indicates that calculations are 
made at intervals of 0.05 rotor radii. 

Figure 5-1 shows the set of computer-generated plots for 
the set of input data shown in the preceding table. The first 
plot shows the wake radius and the wake width as a function 
of the downwind coordinate, x. Since the wake is not near 
the ground, the wake height above the ground is not shown. 

The first plot also lists the important wake parameters for 
the set of plots. The second plot shows Auq as a function of 
X, and the third plot shows Uc/^(s ^ function of x. The 
fourth plot shows the variation of turbine power factor. Since 
the turbine power factor is undefined in Regions I and II, the 
plot is shown downwind of Xj.^ only. The fifth through the 
ninth plots show the lateral wind speed profiles for the down- 
wind locations shown in the table. The altitude of each pro- 
file is shown in rotor radii relative to the turbine hub. 

Since the wake is far from the ground, the vertical wind speed 
profile at y = 0 is symmetrical to the lateral profile at 
Zv = R^. It was, therefore, not shown. 

Figure 5-2 shows a set of similar plots for m = 2.0, and 
figure 5-3 shows the values for m = 3.0. The first plots of 
each of figures 5-1, 5-2, and 5-3 is similar to the plots given 
by Abramovich in figure 5.17 of reference 2. The second plot 
in each of these figures is similar to figure 5.20 in refer- 
ence 2. In the Abramovich plots, the geometric parameters 
are normalized by the initial wake radius, Rq. In figures 
5-1, 5-2, and 5-3, the geometric parameters are normalized by 
the turbine dish radius, R^j. The relationship between Rq 
and R(j is given by equation (3-18). 
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Wakes in ground effect . - The input parameters for the 
wake in ground effect are listed below. The hub of the tur- 
bine is 1.35 rotor radii above the ground. This is typical of 
very large wind turbines. The first lateral wind speed pro- 
file (i = 0; ZQ = 0.35) is at the lower edge of the turbine 
disk, and the third lateral wind speed profile (i = 2) is at 
the rotor hub. 


VALUES OF INPUT PARAMETERS FOR TURBINE WAKE PLOTS 
FOR WAKE IN GROUND EFFECT 


ST = 0.0 

J = 5 

CPM = 1.5 

10 = 1 

AH = 1.35 

NP = 0 

VHO = 1.0 

XNPT(j) = 5., 

WEXP = 0.15 

ZO = 0.35 

RRR = 1.0 

DZZ = 0.5 
DYY =0.05 


Figure 5-4 shows the set of wake plots for the input data 
shown above. Since the wake is near the ground, the height 
of the top of the wake above the ground is also shown in the 
first plot, and a vertical wind speed profile plot is included. 
In the plot of the vertical wind speed profiles, the solid 
line going across the plot is the wake center. The short 
horizontal lines at the left of the plot are the bottom and 
top of the turbine disk. 

In the plot of the normalized wake boundaries (Fig. 5-4(1)), 
the point at which the wake radius equals h , the height of the 
rotor hub above the ground, is the point at which the wake 
intersects the ground plane. Initially, the effect of the ground 
plane affects only the lower part of the wake, but it can spread 
to affect the wake center, although the magnitude of the effect 
at the wake center is negligibly small because of the large wake 
radius (and correspondingly small wid speed deficit) necessary 
for the effect of the ground plane to reach the wake center. 

It is recalled from the development of the turbine power factor 
that the turbine power factor does not include the effect ot the 
ground plane. 
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The effect of the ground on the wind speed profiles is 
seen by comparing the wind speed profiles for x = 50 from 
Figures 5-1 and 5-4. In Figure 5-1, the profiles for altitudes 
of -1.0 and 1.0 rotor radii from the rotor hub are identical 
because of the symmetry of the wake. However, because of the 
ground effect, the wind speed profiles at altitudes of -1.0 and 
1.0 rotor radii relative to the hub are not identical in Figure 
5-4. The lower altitude exhibits a lower wind speed in the wake 
because of the ground effect. At x = 50, it is only the 
lowest altitude which has been affected by the ground. As the 
wake expands further, the effect of the ground affects the wind 
speed profiles of higher altitudes. This is seen in later plots 
where ambient turbulence causes faster expansion of the wake. 

For all of the previous figures, there is no ambient 
turbulence. Figure 5-5 shows a plot set for ambient turbulence 
represented by Pasquill stability class E. The data input 
is identical with that shown in the table, except that ST = -5. 
for Pasquill stability class E. From Table 3-4 arid equation 
(3-28), the wake growth rate due to ambient turbulence is 
a = 0.099. 

Figure 5-6 shows a plot set for Pasquill stability class 
D (ST = -4.). The wake growth rate due to ambient turbulence 
is a = 0.136. Figure 5-7 shows a plot set for Pasquill sta- 
bility class C (ST = -3.). The wake growth rate due to ambient 
turbulence is a = 0.205. 


5.2 Effect of Principal Wake Parameters 

In the following discussion, the effect of the principal 
wake parameters upon the wake is investigated. The parameters 
which are varied are the ambient turbulence as represented by 
the Pasquill stability class; the initial velocity ratio, m; 
the height of the rotor hub, hg.; and the power law coefficient 
of the free stream wind speed profile, y. 

Effect of ambient turbulence . - Figure 5-8 shows the 
effect of ambient turbulence on the wake radius. Figure 5-8 
is a composite of the wake radius as shown in the first plot 
of Figures 5-4 through 5-7. Figure 5-9 shows the effect of 
ambient turbulence of AUr. Figure 5-10 shows the effect of 
ambient turbulence on Uc/U„,. Figure 5-11 shows the effect of 
ambient turbulence on the turbine power factor. It is seen 
that ambient turbulence is a very important factor in the 
recovery of wakes of large wind turbines. From Table 3-3, 
it is recalled that most turbine operation occurs with Pasquill 
stability class D, and Pasquill class D represents the ambient 
turbulence level to which inter- turbine spacing should be 
designed. 
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The turbine power factor is the most significant indicator 
of wake recovery. From Figure 5-11, for Pasquill stability 
class D, the turbine power factor is approximately 0.9 at a 
distance of 35 rotor radii downwind of the turbine. While the 
ambient turbulence has a significant impact on the growth in 
the wake radius, its impact on the turbine power factor is not 
great. For comparison with the 0.90 power factor for Pasquill 
stability class D, the turbine power factor is 0.85 for Pasquill 
stability class E and 0.95 for Pasquill stability class C at 
35 rotor radii downwind of the turbine. 

It is clear that a point of diminishing returns is reached 
in the re-energization of the turbine wake. From Figure 3-7, 
for m = 1.5, a turbine power factor of 0.9 is reached at T 2 = 5.4 
Figures 5-8_and 5-11 confirm this number since Figure 5-11^ 
shows that P = 0.9 at approximately x = 34 for Pasquill stability 
class D, and Figure 5-8 shows that r-=5.4at x = 34 for Pasquill 
stability class D. From Figure 3-7, “^the turbine power factor 
reaches a value of 0.,95at wake radius, r 2 = 8. From extrap- 
olation of the wake radius for Pasquill stability class D in 
Figure 5-8, a wake radius of 8 occurs at x = 56. Thus, a 
distance of 35 wake radii is required for recovery from the 
initial wake to a turbine power factor of 0.9, and an additional 
distance of 21 rotor radii is required for recovery from 0.90 
to 0.95 for the turbine power factor. Certainly, power recovery 
is very much slower after the turbine power factor reaches 
0.9 than it is before it reaches 0.9. 

For the far wake, Figure 3-7 is essentially a plot of 
turbine power factor as a function of downwind distance since 
the wake radius is a straight line function of the downwind 
coordinate, x. In the far wake, the wake growth rate is 
almost entirely due to ambient turbulence with very little 
contribution from mechanical turbulence. For Pasquill 
stability class D, a scaling factor of 1/0.136 should be applied 
to the wake radius axis of Figure 3-7 to convert it to distance 
downvand of the turbine;, s ince the wake growth rate due to 
ambient turbulence is 0.136 for Pasquill stability class D. 
Interpreting Figure 3-7 in this way shows that a very large 
downwind distance is required for recovery of the last few 
percentage points for the turbine power factor. In practice, 
recovery to very high values of turbine power factor may be im- 
practical because of the large additional distances required. 

Effect of initial velocity ratio . - Figure 5-12 shows 
the effect of the initial velocity ratio, m, upon the wake 
radius. Figure 5-13 shows the effect of m upon Au^. Figure 
5-14 shows the effect of m upon Uc/Uoo, and Figure 5-15 shows 
the effect of m upon the turbine power factor. 
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The initial velocity ratio has almost no effect upon 
the wake radius. The downwind extent of Region II, Xj^, is 
almost totally unaffected by m. This is shown in Figures 
5-13, 5-14, and 5-15. Figure 5-13 shows the relative wind 
speed deficit (i.e., the wind speed deficit at the wake center 
normalized by the wind speed deficit at the center of the 
initial wake) . This relative wind speed deficit is not a 
strong function of m. However, because the wind speed 
deficit of the initial wake is much larger for a large value 
of m (by definition of m) , the normalized wind speed at 
the wake center (as shown in Figure 5-14) strongly depends 
on m. Correspondingely , the turbine power factor shown in 
Figure 5-15 is a strong function of m. For the larger values 
of m, more power is extracted by each wind turbine (cf. 
equations (3-9) and 3-11)), but the inter- turbine recovery 
of the wdke requires more space. The lines shown in Figure 
5-15 approximate the cube of the lines shown in Figure 5-14. 

Effect of turbine rotor height . - The wind speed profile 
plots are the only plots which show the effect of the ground 
on the wake. The ground effect can be seen best from the 
vertical wind speed profiles with Y = 0 (i.e., no variation 
in the free stream wind speed with altitude). Figures 5-16, 
5-17, and 5-18 show the vertical wind speed profiles for 
values of h of 1.2, 1.35, and 1.65, respectively. 

a 

In the near wake (x = 5 and x = 10) , the wake is 
approximately symmetric for h^ = 1.35 and h = 1.65. For 
the far wake (x = 20, x = 35,^and x = 50) the wake is clearly 
not symmetric for all values of h . The effect of the ground 
upon the wind speed profiles can be clearly seen by comparing 
the profiles for x = 20 from the three plots. The retarding 
effect of the ground is clearly evident. For comparison. 
Figures 5-19, 5-20, and 5-21 show the same data as Figures 
5-16, 5-17, and 5-18, but with a power law coefficient, y, 
of 0.15. 

Effect of the power law coefficient . - The power law 
exponent, y, only affects the vertical wind speed profile. 

It can be seen on the vertical wind speed profile and on the 
lateral wind speed profile (if NP = 1 so that the wind speed 
in the wake is normalized by the free stream wind speed at 
the hub altitude). The effect of y on the vertical wind 
speed profile may be seen by comparing Figures 5-13, 5-20, and 
5-22. The value, y = 0.25 in Figure 5-22 represents a 
very turbulent lower atmosphere. 
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(1). NORMALIZED WAKE BOUNDARIES 



Figure 5-1. Wake plots for Abramovich solution for m = 1,5. 
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(2). NORMALIZED WIND SPEED DEFICIT 
AT THE WAKE CENTER 
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Figure 5-1. Wake plots for Abramovich solution for m = 1.5 
(continued) . 
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(3). NORMALIZED WIND SPEED 
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Figure 5-1. Wake plots for Abramovich solution for m = 1.5 
(continued) . 
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(4). TURBINE POWER FACTOR 
FOR AN UNBOUNDED WAKE 
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Figure 5 - 1 . Wake plots for Abramovich solution for m = 1.5 
(continued) . 
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( 5). LATERAL WIND SPEED PROFILE 
AT X/R^ = 5.00 





Figure 5-1. 


Wake plots for Abramovich solution for m - 1.5 
(continued) . 




WIND SPEED. U/U 


66 


( 6). LATERAL WIND SPEED PROFILE 
AT X/R^ = 10.00 





Figure 5-1. Wake plots for Abramovich solution for m = 1.5 
(continued) . 
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( 7). LATERAL WIND SPEED PROFILE 


AT X/Rd = 20.00 



Figure 5-1. Wake plots for Abramovich solution for m = 1.5 
(continued) . 
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( 8). LATERAL WIND SPEED PROFILE 
AT X/R^ = 35.00 
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Figure 5-1, Wake plots for Abramovich solution for m = 1.5 
(continued) . 
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( 9). LATERAL WIND SPEED PROFILE 
AT X/R^ = 50.00 
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Figure 5-1. Wake plots for Abramovich solution for m = 1.5 
(concluded) . 
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(1). NORMALIZED WAKE BOUNDARIES 
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Figure 5-2. Wake plots for Abramovich solution for m = 2.0. 
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Figure 5-2. Wake plots for Abramovich solution for m = 2.0 
(continued) . 
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Figure 5-2. Wake plots for Abramovich solution for m = 2.0 
Ccbntinued) . 
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(4). TURBINE POWER FACTOR 
FOR AN UNBOUNDED WAKE 



Figure 5-2. Wake plots for Abramovich solution for m = 2.0 
(continued) . 
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( 5). LATERAL WIND SPEED PROFILE 
AT X/Rd = 5.00 
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Figure 5-2, Wake plots for Abramovich solution for m = 2.0 
Ccontinued) . 
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( 6). LATERAL WIND SPEED PROFILE 
AT X/R^ = 10.00 
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Figure 5-2. Wake plots for Abramovich solution for m = 2.0 
(continued) . 
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Figure 


( 7). LATERAL WIND SPEED PROFILE 
AT X/R(j = 20.00 



5-2. Wake plots for Abramovich solution for m = 2 0 
(continued). 
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( 8). LATERAL WIND SPEED PROFILE 
AT X/R^ = 35.00 
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Figure 5-2. Wake plots for Abramovich solution for m = 2.0 
(continued) . 
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{ 9). LATERAL WIND SPEED PROFILE 
AT X/R^ = 50.00 
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Figure 5-2. Wake plots for Abramovich solution for m = 2.0 
(concluded) . 
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(1). NORMALIZED WAKE BOUNDARIES 



Figure 5-3. Wake plots for Abramovich solution for m = 2.5. 
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Figure 5-3. Wake plots for Abramovich solution for m = 2.5 
(continued) . 
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(3). NORMAUZED WIND SPEED 
AT THE WAKE CENTER 
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Figure 5-3. Wake plots for Abramovich solution for m = 2.5. 
(continued] . 
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(4). TURBINE POWER FACTOR 
FOR AN UNBOUNDED WAKE 
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Figure 5-3. Wake plots for Abramovich solution for m = 2.5 
Ccontinued) . 
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( 5). LATERAL WIND SPEED PROFILE 


AT X/R<j = 5.00 
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Figure 5-3. Wake plots for Abramovich solution for m = 2.5 
(continued) . 
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( 6). LATERAL WIND SPEED PROFILE 
AT X/R^ = 10.00 



Figure 5-3. Wake plots for Abramovich solution for m = 2.5 
Ccontinued) . 
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( 7). LATERAL WIND SPEED PROFILE 
AT X/Rd = 20.00 
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Figure 5-3. Wake plots for Abramovich solution for m = 2.5 
(continued] . 
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( 8). LATERAL WIND SPEED PROFILE 
AT X/R^ = 35.00 



Figure 5-3. Wake plots for Abramovich solution for m = 2.5 
(continued) . 
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( 9). LATERAL WIND SPEED PROFILE 
AT X/R^ = 50.00 



Figure 5-3. Wake plots for Abramovich solution for m = 2.5 
(concluded) . 
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(1). NORMALIZED WAKE BOUNDARIES 
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Figure 5-4. Wake plots for wake in ground effect with no 
ambient turbulence. 
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Figure 5-4. 


Wake plots for wake in ground effect with no 
ambient turbulence (continued) . 
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(3). NORMAUZED WIND SPEED 



Figure 5-4. 


Wake plots for wake in ground effect with no 
ambient turbulence (continued) . 


POWER rACTOR, P/P, 


91 


(4). TURBINE POWER FACTOR 
FOR AN UNBOUNDED WAKE 


3-1 




3-^ 


§. 


oo 


"Jo" 


m 


ISjO 


20j0 




3ao 


26.0 4ao 


~T 

4SA 


eao 


DISTANCE DOWNWIND FROM TURBINE . X/R^ 


Figure 5-4. 


Wake plots for wake in ground effect with no 
ambient turbulence (continued) . 
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( 5). LATERAL WIND SPEED PROFILE 
AT X/R^ = 5.00 
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Figure 5-4. Wake plots for wake in ground effect with no 
ambient turbulence Ccontinued) . 
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( 6). LATERAL WIND SPEED PROFILE 
AT X/R^ = 10.00 
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Figure 5-4. Wake plots for wake in ground effect with no 
ambient turbulence (continued) . 
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( 7). LATERAL WIND SPEED PROFILE 
AT X/R^ = 20.00 



Figure 5-4. Wake plots for wake in ground effect with no 
ambient turbulence (continued) . 
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( 8). LATERAL WIND SPEED PROFILE 
AT X/R<j = 35.00 



Figure 5-4. Wake plots for wake in ground effect with no 
ambient turbulence (continued) . 
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( 9). LATERAL WIND SPEED PROFILE 
AT X/R^ = 50.00 



Figure 5-4. Wake plots for wake in ground effect with no 
ambient turbulence (continued) . 
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(10). VERTICAL WIND SPEED PROFILES 
AT WAKE CENTER 



Figure 5-4. Wake plots for wake in ground effect with no 
ambient turbulence (concluded) . 



WAKE BOUNDARY PARAMETER 


98 



Figure 5-5. 


Wake plots for wake in ground effect for Pasquill 
stability class E. 
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Figure 5-5. 


Wake plots for wake in ground effect for Pasquill 
stability class E (continued) . 
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Figure 5-5. 


Wake plots for wake in ground effect for Pasquill 
stability class E (continued). 
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Figure 5-5. Wake plots for wake in ground effect for Pasquill 
stability class E (continued). 
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( 5). LATERAL WIND SPEED PROFILE 
AT X/Rd = 5.00 
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Figure 5-5. Wake plots for wake in ground effect for Pasquill 
stability class E Ccontinued) . 
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( 6). LATERAL WIND SPEED PROFILE 
AT X/Rj = 10.00 
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Figure 5-5 


Wake plots for wake in ground effect for Pasquill 
stability class E Ccontinued) . 
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( 7). LATERAL WIND SPEED PROFILE 
AT X/Rd = 20.00 
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Figure 5-5. Wake plots for wake in ground effect for Pasquill 
stability class E (continued). 
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( 8). LATERAL WIND SPEED PROFILE 


AT X/Rd = 35.00 
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Figure 5-5. 


Wake plots for wake in ground effect for Pasquill 
stability class E Ccontinued) . 
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( 9). LATERAL WIND SPEED PROFILE 
AT X/Rd = 50,00 
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(10). VERTICAL WIND SPEED PROFILES 
AT WAKE CENTER 
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Figure 5-5. 


Wake plots for wake in ground effect for Pasquill 
stability class E (concluded) . 
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Figure 5-6. Wake plots for wake in ground effect for Pasquill 
stability class D. 


WIND SPEED DEFICIT (Uoo-UcVCUoo-Uq) 


109 


(2). NORMAUZED WIND SPEED DEFICIT 
AT THE WAKE CENTER 
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Figure 5-6. Wake plots for wake in ground effect for Pasquill 
stability class D (continued) . 
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(3). NORMALIZED WIND SPEED 
AT THE WAKE CENTER 
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Figure 5-6. Wake plots for wake in ground effect for Pasquill 
stability class D (continued). 


Q 



Figure 5-6. 


Wake plots for wake in ground effect for Pasquill 
stability class D (continued). 
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( 5). LATERAL WIND SPEED PROFILE 
AT X/R^ = 5.00 

Altitude of wind speed profile in 
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Figure 5-6. 


Wake plots for wake in ground effect for Pasquill 
stability class D (continued). 
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( 6). LATERAL WIND SPEED PROFILE 
AT X/Rji = 10.00 
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Figure 5-6. 


Wake plots for wake in ground effect for Pasquill 
stability class D Ccontinued) . 
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( 7). LATERAL WIND SPEED PROFILE 
AT X/Rd = 20.00 



Figure 5-6. 


Wake plots for wake in ground effect for Pasquill 
stability class D Ccontinued) . 
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( 8). LATERAL WIND SPEED PROFILE 
AT X/Rd = 35.00 
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Figure 5-6. Wake plots for wake in ground effect for Pasquill 
stability class D (continued) . 



WIND SPEED. U/U 


116 


( 9). LATERAL WIND SPEED PROFILE 
AT X/Rd = 50.00 
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Figure 5-6. 


Wake plots for wake in ground effect for Pasquill 
stability class D Ccontinued) . 
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(10). VERTICAL WIND SPEED PROFILES 
AT WAKE CENTER 



Figure 5-6. 


Wake plots for wake in ground effect for Pasquill 
stability class D Cconcluded) . 
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Figure 5-7. Wake plots for wake in ground effect for Pasquill 
stability class C. 
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(2). NORMAUZED WIND SPEED DEFICIT 
AT THE WAKE CENTER 



Figure 5-7. Wake plots for vrake in ground effect for Pasquill 
stability class C (continued). 
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Figure 5-7. 


Wake plots for wake in ground effect for Pasquill 
stability class C Ccontinued) . 
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Figure 5-7. Wake plots for wake in ground effect for Pasquill 
stability class C Ccontinued) . 
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( 5). LATERAL WIND SPEED PROFILE 
AT X/Rd = 5.00 
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Figure 5-7. 


Wake plots for vrake in ground effect for Pasquill 
stability class C Ccontinued}. 
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( 6). LATERAL WIND SPEED PROFILE 
AT X/Rd = 10.00 



Figure 5-7. 


Wake plots for wake in ground effect for Pasquill 
stability class C Ccontinuedl. 
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( 7). LATERAL WIND SPEED PROFILE 


AT X/Rj = 20.00 
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Figure 5-7. Wake plots for wake in ground effect for Pasquill 
stability class C Ccontinued) . 
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( 6). LATERAL WIND SPEED PROFILE 
AT X/Rd = 35.00 
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Figure 5-7. Wake plots for wake in ground effect for Pasquill 
stability class C Ccontinued) . 
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( 9). LATERAL WIND SPEED PROFILE 
AT X/Rd = 50.00 
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Figure 5-7. Wake plots for wake in ground effect for Pasquill 
stability class C (continued). 
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(10). VERTICAL WIND SPEED PROFILES 
AT WAKE CENTER 
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Figure 5-7. Wake plots for wake in ground effect for Pasquill 
stability class C (concluded) . 



m = 1 . 5 



Figure 5-8. 


Effect of ambient turbulence on wake radius. 
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m = 1.5 



Figure 5-9. 


Effect of ambient turbulence on wind speed 
deficit parameter. 
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Figure 5 -10. 


Effect of ambient turbulence on wind speed 
at the wake center. 
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-11. Effect of ambient turbulence on turbine 
power factor. 
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Pasquill stability- 
class D 



Figure 5-12. 


Effect of initial velocity ratio in the wake 
on wake radius . 
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Figure 5-13. 


Effect of initial velocity ratio in the wake 
on wind speed deficit parameter. 
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Figure 5-14. 


Effect of initial velocity ratio in the wake 
on wind speed at the wake center. 
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Figure 5-15. 


Effect of initial velocity ratio in the wake 
on turbine power factor. 
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(10). VERTICAL WIND SPEED PROFILES 
AT WAKE CENTER 
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Figure 5-16. Effect of the gound plane on the vertical 

wind speed profile for a rotor hub altitude 
of 1.2 rotor radii. 
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(10). VERTICAL WIND SPEED PROFILES 
AT WAKE CENTER 
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Figure 5-17. Effect of the ground plane on the vertical 

wind speed profile for a rotor hub altitude 
of 1.35 rotor radii. 



ALTITUDE, Z/R^ 


138 


(10). VERTICAL WIND SPEED PROFILES 
AT WAKE CENTER 



1.0 0.2 0.4 0.8 0.8 1.0 i£ 


WIND SPEED, U/Uoo 


Figure 5-18, 


Effect of the ground plane on the vertical 
wind speed profile for a rotor hub altitude 
of 1.65 rotor radii. 
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(10). VERTICAL WIND SPEED PROFILES 
AT WAKE CENTER 



Figure 5-19. Effect of the ground plane on the vertical 

wind speed profile for a positive wind speed 
profile exponent and a rotor hub altitude of 
1.2 rotor radii. 
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(10). VERTICAL WIND SPEED PROFILES 
AT WAKE CENTER 
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Figure 5-20. Effect of the ground plane on the vertical 

wind speed profile for a positive wind speed 
profile exponent and a rotor hub altitude 
of 1.35 rotor radii. 
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(10). VERTICAL WIND SPEED PROFILES 
AT WAKE CENTER 
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Figure 5-21. Effect of the ground plane on the vertical 

wind speed profile for a positive wind speed 
profile exponent and a rotor hub altitude 
of 1.65 rotor radii. 
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(10). VERTICAL WIND SPEED PROFILES 

AT WAKE CENTER 



Figure 5-22, Effect of the ground plane on the vertical 
wind speed profile for a power law exponent 
for the free stream flow of 0.25. 
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6.0 CONCLUDING REMARKS 

An analytic model for the calculation of the recovery of 
wakes of large wind turbines has been developed. The model is 
based upon the theory of coflowing turbulent jets as developed 
by G. N. Abramovich. The Abramovich model has been modified 
to relate turbine parameters to the wake parameters used by 
Abramovich, to add the effects of ambient turbulence to the model, 
and to calculate the turbine power factor, which is defined as 
the ratio of the power which a turbine in the wake of another 
turbine would generate the power which the turbine would generate 
if it were in the free stream air flow. The theory for the dis- 
persion of pollutants in a turbulent atmosphere as developed 
by F. Pasquill has been adapted to describe wake growth due to 
ambient turbulence, and the methods used by P. B. S. Lissaman 
have been used to combine the effects of wake growth as de- 
scribed by the Abramovich model with wake growth due to ambient 
turbulence . 

This approach is currently considered to be the most ap- 
propriate approach for an analytical model of turbine wake re- 
covery. However, several intrinsic uncertainties about the model 
remain. A numberical approach (e.g., finite difference inte- 
gration of the three-dimensional Navier-Stokes equations) could 
be used. However, it is not considered appropriate, because 
comparable gross uncertainties would exist for this approach 
(e.g., magnitude of the turbulent fluctuations) and large costs 
would be incurred to generate solutions. Any further work on 
the model might best be directed toward expansions or refine- 
ments of factors within the existing formulation. 

The literature has shown only very limited experimental 
evaluation of the model. Comparisons of computer results with 
experimental data are therefore, very much desired to establish 
the adequacy of the model and to resolve the questions involved. 
However, it is recognized that extensive data from full scale 
turbines are very difficult to obtain. An attempt to obtain 
such data is described in the companion report on the second 
(experimental) phase of the research program conducted by Lock- 
heed. Wind turbine model tests in wind tunnels may also be of 
use in validating some aspects of the model formulation. Al- 
though it would be almost impossible to duplicate the Reynolds 
number and geometry of an actual turbine in the wind tunnel, 
meaningful applicable results could be obtained. In another 
approach, the computer program may be used to determine the 
sensitivity of the final results to the factors in question. 



144 


REFERENCES 


1. Schlichting, H. , Boundary Layer Theory , Fourth Edition, 
McGraw-Hill, New York NY, 1962 . 

2. Abramovich, G. N. , The Theory of Turbulent Jets , MIT Press, 
Cambridge MA, 1963. 

3. Templin, R. J. , "An Estimate of the Interaction of Windmills 
in Wind Spread Arrays", Canada NRC, LTR-LA-171, Dec 1974. 

4. Crafoord, C., "An Estimate of the Interaction of a Limited 
Array of Windmills", University of Stockholm Sweden, 

Nov. 1975. 

5. Ljungstrom. 0., "Aerodynamic Interference Between Adjacent 
Wind Generators of the Conventional Horizontal Axis Type", 
Appendix 3-3, "Wind Power in Sweden", 1974. 

6. Builtjes, P., Private communication with P. B. S. Lissaman, 
July 1977. 

7. Lissaman, P. B. S. , "Research Program on Energy Effectivess 
of Arrays of Wind Energy Collection Systems", Report 

AV SP 637, AeroVironment , Inc., Pasadena CA, 1976. 

8. Anon., "Wind Flow Characteristics in the Wake of Large 
Wind Turbines", Report AV FR 8147, AeroVironment, Inc., 
Pasadena CA, Sept 1978. 

9. Lissaman, P. B. S. , "Energy Effectiveness of Arrays of 
Wind Turbines", Journal of Energy, v. 3, n. 6, Nov-Dec 
1979 , 

10. Bird, R. B., W. E. Stewart, and E. L. Lightfoot, Transport 
Phenomena , John Wiley § Sons, Inc., New York NY, 1965. 

11. Turner, D. B. , "Workbook of Atmospheric Dispersion Estimates", 
Environmental Protection Agency, Office of Air Programs, 
Research Triangle Park NC, 1970. 

12. Wilson, R. E., P. B. S. Lissaman, and S. N. Walker, 
"Aerodynamic Performance of Wind Turbines", Oregon State 
University, Corvallis OR, June 1976. 

13. Slade, D. H. (editor), "Meteorology and Atomic Energy," 

U. S. Atomic Energy Commission TID-24190, 1968. 



145 


Appendix A 
LIST OF SYMBOLS 


This appendix contains a list o£ symbols used in the 
development of the analytic model. Where appropriate, 
variable names used in the FORTRAN computer program are 
also listed. The equation number, figure number, or table 
number following the definition is the equation, figure, 
or table where the parameter is defined or first used in 
the description of the analytic model. 

A few parameters used in the deriviation have been 
omitted from this list. They are omitted if they are of 
minor importance and are used for only two or three successive 
equations in the derivation so that the location of their 
definition in the text is never in doubt. 


Symb o 1 
A 


A 


00 


a 


B 


b 

bi,b2 


c 


H 


a 


h 


a 


Computer 

Symbol 


B 


AH 

AH. 


Definition 


Cross-sectional area of the turbine 
disk. Equation (3-1). 

Free stream cross-sectional area 
of the stream tube which, passses 
through the turbine disk. Equation 
(3-64). 

Axial induction factor for the 
turbine. Equation (3-7). 

Parameter used for the derivation 
of the turbine power factor. 
Equation (3-75) . 

Full width of the boundary layer 
in region I. Equation (3-36). 

Partial widths of the boundary 
layer in region I. Figure 3-5. 

Growth rate of the boundary layer 
in region I. Equation (3-59). 

Hub altitude of the turbine in 
physical units. Table 3-2. 

Hub altitude of the turbine nor- 
malized by the rotor radius. 

Table 3-2. 



Symbo l 


Computer 

Symbol 

AH 
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Definition 


h 

I 


i 

J 


K 


m 


lii 


n 


P 


P 


CO 


p 


10 


J 


CPM 


AN 


PR 


Height of the top of the wake above 
the ground. Equation C3-80). 

Input parameter which specifies 
whether geometric input parameters 
are in physical units or normalized 
by the rotor radius. Table 3-2. 

Index to define altitude of lateral 
wind speed profile. Equation (5-1). 

Number of downwind locations at 
which wind speed profiles are to be 
calculated. Table 3-2. 

Constant used in derivation of the 
turbine power factor. Equation 
(3-64) . 

Ratio of the free stream wind speed 
to the initial wind speed in the 
wake (after expansion by potential 
effects to wake radius, Rn) . 

Equation (3-9) . 

Mass flow rate of air passing 
through the turbine disk. Equation 
(3-1). 

Ratio of the downwind extent of 
Region II to the downwind extent 
of Region I. Equation (3-47). 

Power extracted from the air by 
the turbine. Equation (3-63). 

Power extracted from the air by 
a turbine in the free stream. 
Equation (3- 10) . 

Turbine power factor, which is the 
ratio of power extracted by a tur- 
bine in the wake of another turbine 
to the power which an identical 
turbine would extract in the free 
stream flow. Equation (3-63). 
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Computer 

Symbol Symbol 


P 


P 


R 

R. 

c 

R. 


R 


21 


R 


22 


R 


RRR 


r 


Definition 


Static pressure on the upwind 
side of the turbine. Equation 
C3-2). 

Static pressure on the downwind 
side of the turbine. Equation 
(3-2). 

Ambient static pressure. Equation 
(3-3). 

Radial coordinate in physical units. 

Radius of the turbine rotor disk. 
Figure 3-1. 

Initial wake radius in physical 
units. Figure 3-1. 

Radius of the potential core in 
Region I. Figure 3-1. 

Wake radius in physical units. 

Figure 3-1. 

Wake radius at the end of Region 
I in physical units. Figure 3-1. 

Wake radius at the end of Region II 
in physical units. Figure 3-1. 

Free stream radius of the stream 
tube which passes through the 
turbine disk. Figure 3-6. 

Radial coordinate normalized by the 
turbine disk radius . 


r 


sj 


r 


0 


R2SAV Wake radius at the j th downwind 

location at which wind speed profiles 
are to be plotted. Equation (3-84). 

RO Initial wake radius normalized by 

the turbine disk radius. Equation 
(3-18) . 
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Symbol 


r 


1 


r 


2 


r 


2 


I 


r 


21 


r 


22 


r 


CO 


Ar 


T 

U 


U 


0 


Computer 

Symbol Definition 


R1 

R2 

R2P 


R21 


R22 


RI 


DR 


Radius of the potential core in 
Region I normalized by the turbine 
disk radius. Equation (3-36). 

Wake radius normalized by the 
turbine disk radius. Equation 
(3-24) . 

Incremented value of R2 used in 
numerical inetgration for Region 
III. Discussion proceeding 
Equation (3-51). 

Wake radius at the end of Region I 
normalized by the turbine disk 
radius. Equation (3-34). 

Wake radius at the end of Region II 
normalized by the turbine disk 
radius. Equation (3-50). 

Free stream radius (normalized by 
the turbine disk radius) of the 
stream tube which passes through 
the turbine disk. Equation (3-67). 

Increment of r used for numerical 
integration in Region III. 

Equation (3-62) . 

Axial thrust on the turbine disk. 
Equation (3-1) . 

Local wind speed in the wake. 

Equation (3-32) . 

Wind speed at the center of the 
wake. Equation (3-81). 

Wind speed of the air passing 
through the turbine disk. Equation 
(3-1). 

Initial wind speed of the wake (after 
expansion to wake radius, Rq) . 
Equation (3-1) . 
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Computer 

Symbol 


Definition 

Free stream wind speed. Equation 
(3-1). 


XNPT 


Wind speed deficit at the center 
of the wake in physical units. 
Equation (3-53). 

Wind speed deficit in the wake 
normalized by the initial wind 
speed deficit. Equation (3-22). 

Wind speed at the center of the 
wake normalized by the free stream 
wind speed. Equation (3-82). 

Wind speed deficit at the center 
of the wake normalized by the initial 
wind speed deficit at the center 
of the wake. Equation (3-22). 

Wind speed in the wake normalized 
by the free stream wind speed. 
Equation (3-73) . 

Wind speed in the wake normalized 
by the free stream wind speed and 
calculated for the image turbine. 
Equation (3-96) . 

Wake width normalized by the rotor 
radius. Equation (3-79). 

Coordinate distance measured down- 
wind from the turbine in physical 
units . 

Downwind location at which wind 
speed profiles are to be calculated. 
Table 3-2. 



Symbol 


Computer 

Symbol 


Definition 


X 


H 


X 


N 


X 


x’ 






H^m 



Y 


AY 


Ay 


Z 


X 


XP 

XH 

XHM 


XNPT 


DYY 

DYY 


Downwind extent of Region I in 
physical units. Figure 3-1. 

Downwind extent of Region II in 
physical units. Figure 3-1. 

Coordinate distance measured 
downwind from the turbine, 
normalized by the turbine disk 
radius. Equation C3-20). 

Incremented value of x used 
during the numerical integration 
in Region III. Equation (3-62). 

Downwind extent of Region I 
normalized by the turbine disk 
radius. Equation (3-39). 

Downwind extent of Region I from 
Abramovich solution normalized 
by the turbine disk radius. 

Equation (3-47). 

Downwind location at which wind 
speed profiles are to be calculated 
normalized by the rotor radius. 
Equation (3-83) . 

Lateral coordinate measured from 
the vertical center of the wake 
in physical units. 

Increment in the lateral coordinate 
used in generating the wake profile 
in physical units. Table 3-2. 

Increment in the lateral coordinate 
used in generating the wake profile 
normalized by the turbine rotor 
radius. Table 3-2. 

Vertical coordinate measured 
positive upward from the ground 
in physical units. 
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Symbol 

z 


z 


V 


Z 


0 


AZ 


Az 


a 


Y 

AKE 


ri 

P 

a 

^9 


Computer 

Symbol 

Z 

ZVAL 


ZO 

DZZ 


DZZ 


ALPHA 


N 


ST 


Definition 


Vertical coordinated normalized 
by the turbine rotor radius. 

Altitude relative to the turbine 
hub, normalized by the rotor 
radius. Equation (3-85). 

Altitude of the lowest wind speed 
profile in physical units. Table 3-2. 

Altitude increment between 
successive wind speed profiles 
in physical units. Table 3-2. 

Altitude increment between 
successive wind speed profiles 
normalized by the rotor radius. 

Table 3-2. 

Wake growth rate due to ambient 
turbulence. Equation (3-28). 

Exponent for the power law wind 
speed profile. Table 3-2. 

Change in kinetic energy of the 
air passing through the turbine. 
Equation (3-10) . 

Non-dimentional radial parameter 
used in Region I. Equation (3-33). 

Mass density of ambient air. 

Equation 3-1. 

Pollution dispersion coefficient 
from plume theory. Equation (3-19). 

Standard deviation of wind direction 
in the ambient air. Table 3-2. 
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Appendix B 
PROGRAM LISTING 


The following pages contain a listing of the FORTRAN computer 
program used to calculate turbine wake progerties and to generate 
plots of wake properties. Comments in the program listing 
describe the sequence of calculations. Equation numbers given 
in the listing give the equation numbers in the text of this 
report . 

The main program is listed first with suboutines following. 
Subroutine R3 calculates the wind speed deficit parameter, Au^; 
the wake growth rate due to mechanical turbulence, (dr 2 /dx)^; 
and the turbine power factor, P, in Region III. Subroutine™ 

CALCU is used to calculate the wind speed profiles. The 
remaining suboutines are plotting suboutines used to generate 
appropriate plots. Subroutine XVSR2 is used to plot wake 
boundary parameters as a function of the downwind coordinate, x. 
Subroutine XVSDUC is used to plot the wind speed deficit 
parameter, Au^, as a function of x. Subroutine XVSUB is used 
to plot the wind speed at the wake center as a function of x* 
Subroutine XVSPR is used to plot the turbine power factor, P, 
as a function of x. Subroutine YVSUBR is used to generate 
the plots for the lateral wind speed profiles. Subroutine 
UVSZ is used to generate the plots for the vertical wind 
speed profiles. 
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TITLE t (MKE 

AUTHORS i DR. U. EBERLE 
A. C. SCOPPA 

L^KHEED niSSILES AND SPACE CO. 

DATE I UINTER. ISIl 

. . . THIS PROCRAH CALCULATES UAKE PARAHETERS FOR THE RECOUERV 
. . . OF UAKES OF LARGE WIND TURBINES. THE PROGRAH PLOTS UAKE 
. . . RADIUS. UAKE UIDTH, WAKE ALTITUDE. THE UIND SPEED DEFICIT 
. . . AT THE CENTER OF THE UAKE. THE UIND SPEED AT THE CENTER 
. . . OF THE UAKE. AND THE TURBINE POUER FACTOR AS FUNCTIONS OF 
. . . THE DISTANCE DOUNUIND Of THE TURBINE. IN ADDITION. UIND SPEED 
. . . PROFILES OF THE UIND SPEED AS A FUNCTION OF THE LATERAL 
. . . COORDINATE ARE PLOTTED FOR EIGHT ALTITUDES FOR EACH DOUNUIND 
. . . LOCATION SELECTED BV THE USER. A VERTICAL UIND SPEED PROFIU 
. . . IS GENERATED FOR EACH DOUNUIND LOCATION SELECTED FOR PLOTS OF 
. . . OF LATERAL UIND SPEED PROFILES. 


. . . ST - ATfWSPHERIC TURBULENCE PARAHATER. IF INPUT AS A 

POSITIVE NUNBER. IT IS THE STANDARD DEVIATION OF THE 

UIND DIRECTION. IF INPUT AS A NEGATIVE NUNBER, IT IS 

THE PASOUILL STABILITV CLASS AS FOLLOUSi 

-I FOR PASOUILL STABILITV CLASS A 

-2 FOR PASOUILL STABILITV CLASS B 

-3 FOR PASOUILL STABILITV CLASS C 

'A FOR PASOUILL STABILITV CLASS D 

-5 FOR PASOUILL STABILITV CLASS E 

-6 FOR PASOUILL STABILITV CLASS F 

IF ST IS INPUT AS AN ANBIENT TURBULENCE OF ZERO 

IS USED, AND THE RESULTING PROFILES ARE THE 

ABRANOVICH SOLUTIONS. 

. . . CPN - RATIO OF FREE STREAN UIND SPEED TO INITIAL UIND SPEED 
IN THE UAKE 

. . . AH - HUB HEIGHT IN ROTOR RADII OR IN PHYSICAL UNITS AS 
SPECIFIED BV 10 

. . . VH« - ANBIENT UIND SPEED AT HUB ALTITUDE . THIS NAV BE CIVEH 
IN PHYSICAL UNITS IF OUTPUT IS DESIRED IN PHYSICAL 
UNITS . AN INPUT VALUE OF !.• WILL GIVE OUTPUT 
NORNALIZCD BY THE FREE STREAN UIND STREAN. 

. . . UEXP - COEFFICIENT OF THE UIND SPEED POUER LAU OF THE FORN 
*V>VB<(Z/Z*)SSUEXP* 

. . . RRR - ROTOR RADIUS. THIS NAV BE GIVEN IN PHYSICAL UNITS. 

AN INPUT VALUE OF I.S UILL GIVE ALL OUTPUT 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Of UNITS OF LENGTH IN ROTOR RADII. AN INPUT 

VALUE Of .5 WILL GIUE ALL OUTPUT OF UNITS 

OF LENGTH IN ROTOR DIANETERS. AN INPUT VALUE GREATER 

THAN 2.S UILL GIVE ALL OUTPUT OF UNITS OF LENGTH IN THE 

PHYSICAL UNITS USED FOR THE ROTOR RADIUS. 

. J - NUniER OF DOUN UIND LOCATIONS AT UHICH UIND 

SPEED PROFILES ARE TO IE CALCULATED 

. 10 - INPUT OPTION FOR PARANETERS UITH PHYSICAL DINENSIONS 

(AH.XNPT.ZO.DZZ.DYY) 

1 FOR INPUT IN ROTOR RADII 

2 FOR INPUT IN PHYSICAL UNITS 

. HP - SPECIFIES HOU VELOCITY IS TO BE N0RNALI2ED FOR 

PLOTS Of WIND SPEED PROFILES IN LATERAL DIRECTION. 

NP • 9 NORNALIZES UIND SPEED BY THE FREE STREAN UIND 
SPEED AT THAT ALTITUDE. NP • 1 NORNALIZES UIND SPEED 
BY THE FREE STREAM UIND SPEED AT THE HUB ALTITUDE. 

. XNPT - DOUNUIND LOCATIONS AT UHICH UIND SPEED PROFILES 

ARE TO BE CALCULATED (ROTOR RADII OR PHYSICAL UNITS 
AS SPECIFIED BY 10 

. Z9 - niNinun altitude at uhich uind speed profiles are to be 

- CALCULATED (ROTOR RADII OR PHYSICAL UNITS AS SPECIFIED 
BY 10) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DZZ - INCRERENT IN ALTITUDES AT WHICH UIND SPEED PROFILES ARE 
TO BE CALCULATED (ROTOR RADII OR PHYSICAL UNITS AS 
SPECIFIED BY 10) 

DVY - INCREMENT IN Y BY UHICH CALCULATIONS ARE TO BE HADE IN 

GENERATING THE UIND SPEED PROFILES (ROTOR RADII OR PHYSICAL 
UNITS AS SPECIFIED BY 10) 

CONTROL INPUT i 

RECORD 1 - (50F19.4) 

C6 INPUT UARIABLESD 
ST,CPn,AH,VH0,UEXP,RRR 

RECORD 2 - (5912) 

C3 UARIABLES3 
J.IO.NP 


RECORD 3 - (59F19.4) 

CJ VARIABLES3 
XNPT 

RECORD 4 - (50F19.4) 

C3 VARIABLES3 
Z9.D22.DYY 

OUTPUT FILES » 

ONE DATA SET PER PLOT, NANING CONVENTION 'SUBNANE.DAT* 
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RISCELLAHCOUS t 

THE PlOTTINa PACKAGE 'OISSPLA* RUST •£ LINKED WITH 
•UAKE* BEFORE EXECUTION. THIS PROCRAH 
CAN BE RUN ‘AS IS* USING A POP-IB CONPILER AND A 
TEKTRONIX 4«14 EXTENDED GRAPHICS TERfllNAL. OTHERWISE 
EXTENSIUE nOOIFICATIONS UILL BE NECESSARY. 

COnnON /SUBCOn/ CPn.AH.Rt.RZ.OUC.UB.UBG.ORDX.PR. 

X XN.R22.R21.XH 

COfINON /PSCALE/' RRR.UHB 

COrmON /CRPA/ICRT, IPAPER 

DOUBLE PRECISION INFIU.UFILE 

REAL XNPT(20),R2SAU(2e) 

DATA ICRT/1 /, FLAG/-99 E/2. 7182818/ 

DATA lUBAR/'UBAR'/ 

. . . DETERNINE DATA SET SPECS FOR CONTROL INPUT 
URITE(5.5) 

FORNATC' ENTER LOGICAL UNIT AND FILE NANE UHERE'. 

* ' CONTROL INPUT RESIDES i') 

READ(5,7)INU. INFILE 
FORNATCIE.AS) 

. . . OPEN DATA SET CONTAINING CONTROL INPUT 

OPEN (UNIT-INU,FILE-IHFILE,ACCESS-'SEQIH',DISPOSE-'SAUE', 

* DEUICE-'DSK' ) 

. . . READ IN CONTROL DATA 

READ( INU, 1« ) ST.CPN, AH.UHB.UEXP.RRR 
FORNATCSef ia.4) 

READ(INU.20) J.IO.NP 
F0RnAT(Sei2> 

READ(INU.IB) (XNPTCI ).I«1.J) 

READ(INU.IB) Ze.DZZ.OVV 

. . . CALCULATE PARANETERS IN ROTOR RADII IF THEY HAUE 
. . . BEEN INPUT IN PHYSICAL UNITS 

IF(.NOT.(IO .EG. an GO TO 26 
C EOUATION (3-12) 

AH«AH/RRR 
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OZZ*DZZ/RRR 

DW-DVV/RRR 
00 86 I>1,J 

XNPTC 1 )»XNPTm/RRR 
COMTIMUE 
CONTIHUe 


EOtMTION (3-13) 
EQUATION (3-14) 
EQUATION (3-15) 


EQUATION (3-83) 


• . CALCUtATE CONSTANTS 


. .SET VALUE FOR ALPHA 
CONTINUE 
IF(ST .CT. 8. ) 

t ALPHA-1. 37X9. 031<E*X(.eStST) 

IF(ABSCST) .LT. l.E-19) 

X ALPHA-9. 

IF(A«S(ST*1. ) .LT. l.E-19) 

X ALPHA-1. 97X. 812 

IF(ABS(ST+8. ) .LT. l.E-19) 

X ALPHA-1. 97X. 156 

IFCABS(ST+3. ) .LT. l.E-19) 

X ALPHA-l.97X.194 

1F(ABS(ST*4. ) .LT. l.E-19) 

X ALPHA- 1.97X. 969 

IF(ABS(ST+5. ) .LT. l.E-19) 

X ALPHA-1.97X.059 

IF(ABS(ST+S. ) .LT. l.E-19) 

X ALPHA- 1.97X. 934 

CONTINUE 

. . SET VALUE FOR RADIUS CHANGE 
DR-.98 


CALCULATE INITIAL UAKE RADIUS IN ROTOR RADII 


EQUATION (3-28) 


R9-SQRT((CPN'fl. )/2. ) 

CALCULATE UAKE RADIUS AT END OF REGION I 


EQUATION (3-18) 


R21-R9/SQRT( . 8144'. 144XCPN ) 


EQUATION (3-34) 


CALCULATE DOUNUIND E)(TENT OF REGION I BV ABRANOVICH 


EQUATION (3-36) 
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c 

c . 
c . 
c 

c 

c . 

c 

c 


c 

c 


40 

c 

c . . 
c . . 
c 


c 

54 

55 
C 

C . 

c . 
c . 
c . 
c . 
c . 
c . 
c 


c 

c 


XHf1-R0«(l.+CPH)/ 

* (.a7<(CPn-i)ssoRT(.2i44.i44«cprm 

. . CALCUUTE OOUNUINO EXTENT OF REGION I IN PRESENCE OF 
. . AMBIENT TURBULENCE 

XH«0.5XR21/SaRT((e.5SR21/XHn)>ta4ALPHAt<2> 

, . CALCULATE UAKE PARANETER AT THE END OF REGION II 

EQUATION (3-47) 

AN-SQRT( .214+. 144XCPN )* 
t (I.-SORT( .134*.124»CPN))/ 

( ((l.-SGRT(.214T.144tCPn))I 

t S0RT(.134+.124XCPN)) 

EQUATION (3-46) 

XN-ANXXH 

EQUATION (3-49) 

R22>Re«AN«(R21-R0) 

X*XN 
R2*R22 
CALL R3 
CONTINUE 

. 5AUE UALUES OF UAKE RADIUS AT INPUT VALUES OF 
. XNPT(I) FOR THE XNPT(I) UHICH LIES IN REGIONS I OR III 

DO 55 I-l.J 

IF(.NOT.(XNPT(I) .LE. XN)> GO TO 54 

EQUATION (3-84) 

R2SAU( I ) -R0* ( R22-R0 )tXNPT ( I )/XN 
CONTINUE 
CONTINUE 


. CREATE FIRST SET OF PLOT FILES FOR UAKE RADIUS, 

. UAKE UIDTH, UAKE ALTITUDE. UIND SPEED DEFICIT AT THE CENTER 
. OF THE UAKE, UIND SPEED AT THE CENTER OF THE UAKE, AND 
. THE TURBINE POUER FACTOR AS A FUNCTION OF THE DISTANCE 
. DOUN UIND OF THE TURBINE 


OPEN(UNIT-10. DEVICE* 'DSK‘,ACCESS*'SEQOUT', 

NODE- ' ASCII '.DISPOSE* 'SAVE'. FILE*'XUSR2. DAT' ) 
OPEN(UNIT*ll, DEVICE* 'DSK',ACCESS*'SEQOUT', 

nOOE*'ASCII'.OISPOSE-'SAUE',FILE*'XVSDUC.OAT' ) 
0P£N(UNIT*12, DEVICE* 'DSK',ACCESS*'SEQOUT', 

NODE* 'ASCI I '.DISPOSE* 'SAVE '.FILE* 'XVSUB.DAT' ) 
OPEN(UNIT* 13, DEVICE* 'DSK', ACCESS* 'SEQOUT', 

fWDE*'ASCII' .DISPOSE* ' SAVE '.FILE-' XVSUBG . DAT ' ) 
OPEN(UNIT*14. DEVICE* 'DSK '.ACCESS* 'SEQOUT'. 

NODE* 'ASCII '.DISPOSE* 'SAVE '.FILE* 'XVSPR.DAT' ) 


FIRSTX-0. 
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FIRSTVR* 

FVa-R«*flH 
FY3*R«ta 

URITEClt.l*) FIRSTX.FIRSTV.FVa.Fva 
FIRSTV-t. 

URITECll.l*) FIRSTX.FIRSTV 
FIRSTV^l./CPfl 

URlTE<ia.l6) FIRSTX.FIRSTV 
FIRSTV-l./CPf1 

UHITE(13.ie) FIRSTX.FIRSTV 

XOLD-X 
RaOLD'Ra 

. . . SEARCH THROUGH X. TERfllHATE UHEN R2 IS GREATER THAN 8 . 

... OR X IS GREATER THAN 52 

C 

6« CONTINUE 

IF(.N0T.(R2 ,LE. a. .AND. X .LE. 52.)) GO TO 70 
IF(.N0T.(X .EQ. 0.))GO TO 61 
R2-R0 
0UC«l. 

UB-l./CPN 

UBG«l./CPN 

61 CONTINUE 

IF<,NOT.(X .EQ. XN)) GO TO 62 
R2*R22 
DUC-1. 

UB«i./cpn 

UBG»l./CPn 

62 CONTINUE 
RADDAH>R2«AH 
RNUl.2*R2S2. 

URITE(ie.ie) X.R2.RADDAH.RmJl.2 
URITEUl.lO) X.DUC 
URITE(12.10) X.UB 
URITEIia.lO) X.UBG 
IFC.NOT.CX .NE. 0.)) GO TO 63 
URITE(14.10) X.PR 

63 CONTINUE 

DO 69 1*1. J 

IF(.NOT.(XNPT(I) .GT. XN ) ) GO TO 68 
XF(.NOT.<X .GT. XNPT(I) 
t .AND. 

t XOLO .LT. XNPT(l))) GO TO 66 

RaNEU<((R2-R2OLD)/<X-XOlD))B(XNPT(Z)-X0LD)4^R20LO 

R2SAU(1)>R2NEU 



ooooooooo o o o o o o o o o o o 
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66 CONTINUE 

IF(.NOT.(X .CO. XNPTd))) GO TO 67 
RSSAUd )-Ra 

67 CONTINUE 

68 CONTINUE 

69 CONTINUE 
XOLD*X 
RaOLO*R2 
R2P-R2+DR 
CALL R3 

DROXE*SORT( ( (DRDX-»ORDXP)/a. )»a'»ALPHAIxa> 

XP<X*DR/0R0XE 

Ra*Rap 

X»XP 

OROX'DROXP 
GO TO 60 

70 CONTINUE 

. . . CLOSE FILES 

CLOSE(UNIT«10> 

cLosE(UNir*in 

CLOSE(UNIT-ia) 

CL0SE(UNIT-13) 

CLOSE (UNIT- 14) 

. . . INITIALIZE PLOTTER FOR OISSPLA, 1300 BAUD RATE. EXTENDED 
. . . GRAPHICS TERNINAL 

CALL TKTRN(ia0.n 

. . . PRODUCE PLOTS 

CALL XUSR2 
CALL XUSDUC 
CALL XUSUB 
CALL XUSUBG 
CALL XUSPR 


GENERATE WIND SPEED PROFILE PLOTS AS A FUNCTION OF THE 
LATERAL COORDINATE. FOR EIGHT ALTITUDES FOR EACH 
DOWNWIND LOCATION SEUCTED BY THE USER. ALSO. A 
VERTICAL WIND SPEED PROFIU IS GENERATED FOR EACH 
DOWNWIND LOCATION SELECTED. 


IUNITZ-J+14+1 

OPEN<UNIT-IUNITZ. FILE- 'UVSZ'.DEVICE-'DSK'. NODE- 'ASCII', 
X DISPOSE- 'SAVE'. ACCESS- 'SEOOUT' ) 

DO 160 I-l.J 

WRITEdUNlTZ.10) FLAG 
IUNIT-1*14 

ENC0DE(6.75.UF1L£) lUBAR.lUNlT 
76 F0RflAT(A4.ia) 
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c 

c 

c 


9 « 


134 


140 


OPEN(UtUr>lUNIT.DEUtCE-'DSIC',ACCCSS-'SEOOUT', 

nODE>'ASCII‘.DlSPOSE*'SAUC',EILE>UmE) 

Z-Z0 

Y-0. 

STEP THROUGH EIGHT ALTITUDES 

00 140 JJ«1.8 

URITE(1UHIT.10) FLAG 

V>0. 

UBAR*0. 

CONTINUE 

IF(.NOT.(ABS(l.-UtAR) .GT. .00001)) GO TO 130 
VUAL-V 

EQUATION (3-85) 

2UAL-AH-Z 

CALL CALCU ( YU AL . 2UAL , XNPT. RESAU, 1 . UB ARA > 
YUAL-Y 

EQUATION (3-95) 

ZUAL*AH4Z 

CALL CALCU ( VUAL. ZUAL.XNPT.R2SAU, I .UBARB ) 

EQUATION (3-98) 

IF(NP .EQ. 0) 

UBAR-UH0<(UBARA>UBARB-1. ) 

EQUATIONS (3-98) AND (3-99) 

IF(NP .EQ. 1) 

UBAR«UH0t(UBARA>UBARB-l. )« 

(Z/AH)ttUEXP 
URITE(IUNIT.ie) V.U8AR 
V.V+DYY 
GO TO 90 
CONTINUE 
V*6. 

UBAR>1. 

URITE(IUNIT,10)Y.UBAR 

2*Z*0ZZ 

CONTINUE 


CALCULATE UERTICAL UINO SPEED PROFILES 


YUAL*0. 

Z>0. 

145 CONTINUE 

IF(.NOT.(Z .LE. 5. ))GO TO 150 

C EQUATION (3-85) 

ZUAL*AH-Z 

CALL C ALCU ( VUAL , ZUAL , XNPT. R2SAU, I, UBARA ) 
ZUAL'AH'I’Z 

C EQUATION (3-95) 

CALL CALCU( VUAL. ZUAL. XNPT, R2SAU.1, UBARB) 

C EQUATIONS (3-98) AND (3-99) 

UBAR«UH0J:(UBARA«UBARB-l.)t(Z/AH)MUEXP 
URITEdUNITZ.lO) UBAR.Z 
Z«ZtDVY 
QO TO 14S 
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15« 

C 

C 

c . 
c 

c 

16« 

C 

C 

C . 
C 

C 

c . 
c 

c 


CONTINUE 

CL0SE(UNIT«IUN1T> 

FOR EACH OF THE •XHFT*, PLOT LATERAL UXNO SPEED PROFIUS 
CALL WSUIR(IUNIT,UFIL£,XHPT.I) 

COHTINUE 

CLOSE(UHIT«IUNITZ> 

PLOT VERTICAL WIND SPEED PROFILES 
CALL UVSZCIUHITZ, 'UVSZ'.AH.J) 

TERfllNATE PLOTTINC SESSION, END PROGRAN 

CALL DONEPL 

END 
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c 

c . 
c . 
c 


c 

c 

c . 

c 

c 


c 

c . 

c 

c 

c 

c . 
c 

c 

c 

c . 
c 


le 

c 

c . . 
c . . 
c . . 
c . . 

c 

c 


SUBROUriNC R3 

. . THIS SUBROUTINE CALCULATES UAKE CROUTH RATE, UIND SPEED AT THE CENTER OF 
. . THE UAKE, AND POWER RATIO IN REGION III OP THE UAKE. 

CONNON /SUBCOn/ CPN, AH.R«,R2,0UC.UB,UB6,DRDX,PR, 
t XN,R22,R21,XH 


. . CALCULATE UIND SPEED DEFICIT AT THE CENTER OF THE UAKE 

EQUATION (3-53) 

DUC-3 . 73* ( - . 2S8»CPn/ ( 1 . -CPN )-SQRT ( ( . 258ICPN/ ( 1 . -CPN ) )««2 
t >.53G/(R2<RES(I.-CPN)/R9«<2))) 

. . CALCULATE UAKE GROWTH RATE DUE TO NECHANICAL TURBULENCE 

EQUATION (3-6B) 

DRDX-.27/(2.*CPN/((CPN-l. )<DUC)-1.) 

. . CALCULATE NON-OINENSIONAL UIND SPEED AT THE UAKE CENTER 

EQUATION (3-94) 


UB-l.-DUC»(l.-i./CPN) 

UBG>UB 


. IS THE CENTER OF THE UAKE IN GROUND EFFECT 7 

IF(.N0T.(R2 .QT. Z.tAH)) GO TO 10 

UF- 1 . -DUC« ( 1 . - 1 . /CPN J « ( 1 . - ( 2 . SAH/'R2 Xtl .5 )»2 
UBG«UB«UF 
CONTINUE 

. CALCULATE THE RATIO OF POWER GENERATED IN A GEORETRICALLV 
. IDENTICAL TURBINE CENTERED IN THE UAKE OF THE UPUINO TURBINE 
. TO THE POWER THAT TURBINE WOULD GENERATE IF IT WERE IN THE 
. STREAN UIND 

EQUATION (3-75) 

B*DUC»<l.-l./CPn) 

EQUATION (3-78) 

Fl-(1 .-3.SB>3.SBS«2-l.tBtt3)/2. 

F2*(6.TB-12.SB»2«6.XBtx3)/’3.5 

F3«(-3.<BT18.tBS*2-15.SBXt3)/S. 

F4>(-12.«B«*2T20.*Bt«3)/6.S 

FS«(3tB**2-15.SBXS3)/8. 

F6-6.<BX«3/9.S 
F7-Btt3/(-ll. ) 

RI«SORT( CCPN»1. )/(2.»CPN)) 

RIR*RI/R2 

PR*2.<(Fl*F2XRIRtXl .5«F3«RIRX<3 

TF4«RIR»4.5FFS<RIRtX6«F6tRIRtX7.54F7<RIR><9> 

RETURN 

END 



163 


c 

c . 
c . 
c . 
c . 
c 


c 

c 

c 

c . 
c 


lee 

c 

c . 
c 


lie 

c 

c . 
c 


SUBROUTINE CALCU(y.2,XNPT0n.R2SU0rt,l,UBA<n 

. . THIS SUBROUTINE CALCULATES UINO SPEED IN THE MAKE FOR 
. . CIUEN UALUES OF THE LATERAL AND VERTICAL COORDINATE 
. . DISTANCES . THE FORN OF THE UINO SPEED PROFILE OF 
. . REGION I . REGION II . OR REGION III AS APPROPRIATE 

CONNON /SUBCOn/ CPN. AH.Re.R2,DUC,UB.UBG,DRDX,PR, 

X XN.R22,Ret,XH 

REAL XNPTDn<2e},N,R2SUDN(2e) 


R>SGRT(V«t2^Z*t2) 


I EQUATION O-SS) 


. . detehuine if xnpt is in region i 


IF{.NOT.<XNPTDrUI) .LT. XH)) GO TO 100 
R1 -Re-ReFXNPTDm I )/XH 
R2-R2SUDN(n 
N*(R-R2)/(R1-R2) 

IF(N .LT. e.) UBAR'l. 

IF(N .GT. 1.) UBAR-l./CPN 

IF(0. .LT. N .AND. N .LT. 1.) 
t UBAR-l./CPNtd.-l./CPN) 

t «(1.-N»1.5)t»2 

CONTINUE 


• EQUATION 0-87) 

! EQUATION 0-33) 

• EQUATION (3-89) 
) EQUATION (3-88) 


• EQUATION (3-96) 


. . DETERHINE IF XNPT IS IN REGION III 

IF(.NOT.(XNPTDN(I) .GT. XH ) ) GO TO 110 
R2-R2SUDn(I) 

CALL R3 

IF(R .GT. R2) UBAR'l > EQUATION (3-91) 

IF(R .LE. R2) UBAR-l.-DUC»(l.-l./CPri) 

X X(l.-(R/R2)t*1.5)tS2 ! EQUATION (3-93) 

CONTINUE 


. . DETERNINE IF XNPT IS IN REGION II 


X 

X 


X 

X 

X 


IF(.NOT.(XH .LT. XNPTDncl) .AND. XNPTDfl(I) 
Rt-e. 

R2-R2SUDN(I) 

H»(R-R2)/’(R1-R2) 

IF(N .LT. 0. ) UBARl'l. 

IF(N .GT. 1.) UBARl-l./CPN 
IF(0. .LT. N .AND. N .LT. 1.) 

UBARl . 1 . /CPN+ ( 1 . -1 . /CPN ) 
X(l.-N*tl.S)St2 
DUC«1. 

IF(R .GT. R2) UBAR2*1. 

IF(R .LE. R2) 

UBAR2«l.-DUC<(l.-l./CPfl) 
t(l.-(R/R2)»1.5Xt2 
UB AR • ( ( XNPTDfl ( I )-XH )/ ( XN-XH ) )SUBAR2 
♦( (XN-XNPTDN(1))/(XN-XH))SUBAR1 


.LT. XN)) GO TO 120 

• EQUATION (3-33) 

! EQUATION 0-89) 

I EQUATION 0-88) 

< EQUATION 0-90) 

• EQUATION 0-94) 

! EQUATION 0-93) 

! EQUATION 0-95) 


C CALL lOUAIT 

CALL HDCOPV 
C 

C . . . TERNINATE PLOT 
C 

CALL ENOPL(0) 

C 

C . . . RETURN TO HAINLINE 
C 

RETURN 

END 
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c 

c 

c 

c 

c 

c . 
c 


c 


5 

c 


c 

c . 
c 


c 

c . 
c 


c 

c . 
c . 
c . 
c . 
c 

c . 
c 


SUBROUTINE XUSR2 
COBfIOH /PSCALE/RRR.UHB 
REAL VUALO) 


NLINES*3 


. . SET riAxinun values for x and v axis 

XHAX-S0. 

vnAX'S. 

IF(.NOT.(RRR .CT. 2.)) GO TO 5 
XflAX'SO.XRRR 
VnAX«8.tRRR 
CONTINUE 


XSCALE-XNAX/SO. 

VSCALE-VNAX/8. 

. . RE-INITIALIZE PLOTTER 

CALL BONPLU) 

CALL PHVSORd.,1* » 
CALL PAGE(l'4.,ll. ) 


. . PLOT AXIS 

CALL HEICHT(0.2» 

CALL CONPLX 

CALL nXlALF( 'STAND','#') 

CALL nx2ALFl'L/CSTD', '•') 

CALL HX3ALF('INSTR','fc') 

CALL TITLEt' ',1, 

t 'DISTANCE DOUNUIND FROtl TURBINE , X/RlLHl.«Dt>', 

» I«a,'UAICE BOUNDARY PARANETERB', 

S 100,10. ,8.) 

CALL GRAPH 1 0 . , Xfl AX/" 1 0 . , 0 . , YNAX^S . ) 

CALL HEADIN('(fl<). NORHALIZED UAKE BOUNDARIES 8',100,1.,1) 


! ! DRAU ONE LINE EACH FOR UAKE RADIUS, UAKE 
. . RADIUS PLUS ALTITUDE. AND TUO TINES UAKE RADIUS 


. . OPEN DATA SET FOR EACH LINE 
DO 40 I’l.NLINES 

OPEN (UNIT* 10, FILE* 'XUSR8.DAT ' , ACCESS* 'SEQIN ' ) 


C 

C 

C 


• • 0 


READ FIRST CARD 

READ(10,10,END*30> KUAL, (YVAL( J),J*1,NL1NE3) 
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le F0«niAT(5»Flt.4) 

XPOS« (1*. /XfMX )XXUALSXSCALE 
VPOS* « ». /-VriAX >»VUAL( I )*VSCALE 
CALt STRTPT(XPOS.VPOS) 

C 

C . . . PLOT ENTIRE DATA SET 
C 

29 CONTINUE 

XPOS* ( 19 . /XNAX )tXUAUXSCALE 
VP0S-(8./VnAX)*VUAL(I )*VSCALE 
CALL CONNPTIXPOS.VPOS) 

READ(19,19,END«39> XUAL, (VUAL(J),J>1,NLINE5I 
GO TO 29 

39 CONTINUE 

C 

C . . . CLOSE DATA SET 
C 

CL0SE(UNIT«19) 

C 

49 CONTINUE 

C 

C . . . HAKE A HARD COPY AND ERASE SCREEN 
C 

C CALL lOUAIT 

CALL HDCOPV 
C 

C . . . TERHINATE PLOT 
C 

CALL ENDPL(t) 

C 

C . . . RETURN TO NAIN ROUTINE 
C 

RETURN 

END 
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SUIROUTINE XUSOUC 
C 

COnnON /PSCALE/RRR.UH4 
C 

C ... SET RAXinun VALUES FOR X AND V AXIS 
C 

XHAX'SS. 

YfWX-l. 

C 

IF(.NOT.(RRR .CT. 2.)) 00 TO S 
XMAX-SO.SRRR 

5 CONTINUE 
C 

IF(.NOT.(UH« .GT. 1.)) GO TO fi 
VNAX'l.SUHG 

6 CONTINUE 
C 

XSCALE*xnAX/S0. 

VSCALE-VNAX/l. 

C 

C. . . . OPEN DATA SET 
C 

OPEN ( UN IT • 1 1 , F I LE • ' XUS DUC . 0 AT ' , ACCESS* ' SEOIN' ) 

C 

C . . . RE-INITIALI2E PLOTTER 
C 

CALL BGNPL(l) 

CALL PHVS0R<1.,1.» 

CALL PACE(H..ll. > 

C 

C . . . PLOT AXIS 
C 

CALL HEIGHT(0.2] 

CALL COnPLX 

CALL NX1ALF< 'STAND', '»' ) 

CALL nx2ALF('L/CSTD','S' ) 

CALL nX3ALF< 'NATHENATIC','*' ) 

CALL flX4ALF('INSTR'.'»' ) 

CALL ZtUSE('«(ULLHl.SOSLXHXt-ULLHl.«CSLXHXt)/’S',ieG) 

CALL Z2USE( 'tCUlLHl .XOlLXHXt-UtLHl .tOSLXHX»)S',ie«) 

CALL TITLE(' ',1, 

t 'DISTANCE DOUNUIND FRON TURBINE , X/RSLHI.«Dtt', 

* ie«, 'WIND SPEED DEFICIT tZUZES', 

S 100,10. ,5.) 

CALL GRAPH ( 0 . , XNAX/ 10 . , 0 . . VNAX/S . ) 

CALL HEADIN('(02«). NORNALI2EO UINO SPEED DEFICITS', 100,1., 2) 
CALL HEADINC'AT THE UAKE CENTERS ',100,1., 2) 

C 

C . . . READ FIRST CARD 
C 

READ(11,10,END>30) XUAL,VUAL 
10 F0RnAT(2F10.<4) 

XPOS- ( 10 ./XNAX }SXVAL>XSCALE 
VPOS- ( 5 ./VNAX )«VVALtVSCALE 
CALL STRTPTI XPOS, VPOS) • 



167 


c 

C . . . PLOT ENTIRE DATA SET 
C 

a« CONTINUE 

XP0S> (!• . /XHAX XXUAUXSCALE 
VPOS- < S . / VNAX ) tVUALXVSCALE 
CALL CONNPT(XPOS.VPOS) 
REA0(n,te,EN0«3«) XVAL.YVAL 
GO TO 2« 

39 CONTINUE 

C 

C . . . WAKE HARD COPY THEN ERASE SCREEN 
C 

C CALL lOUAlT 

CALL HDCOPY 
C 

C . . . CLOSE DATA SET 
C 

CLOSE(UNIT-Il) 

C 

C . . . TERNINATE PLOT 
C 

CALL ENDPL(0) 

C 

C . . . RETURN TO NAIN ROUTINE 
C 

RETURN 

END 
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SUIROUriNE XUSUB 

c 

COmON /PSCALE/RRR.UHB 
C 

c ... SET nAxinun ualues for x and v axis 
c 

XHAX'SB. 

VHAX-1. 

C 

IF(.NOT.(RRR .QT. 2.)) GO TO 5 
XnAX*S«.SRRR 

5 CONTINUE 
C 

IF(.NOT.(UH« .GT. 1.)) GO TO 6 
VHAX't.SUHe 

6 CONTINUE 
C 

XSCAIE'XNAX/Sa. 

VSCALE-'/NAX/l. 

C 

C . . . OPEN DATA SET 
C 

OPEN <UNIT*12.FILE*'XUSUB.DAT'.ACCESS»'SEQ1N') 

C 

C . . . RE-INITIALIZE PLOTTER 
C 

CALL BGNPL(I) 

CALL PHVSORa.,1.) 

CALL PACE(14.,11. ) 

C 

C . . . PLOT AXIS 
C 

CALL HEIGHT(«.a) 

CALL CONPLX 

CALL nXlALF< 'STAND', '«*) 

CALL nX2ALF('L/CSTD','l') 

CALL nX3ALF('INSTR','l' ) 

CALL nX-4ALF<'NATHEnATIC'.'X') 

CALL Z1USE( 'tU&LHl.«C&LXHX«/ULLHl.X0*t',ie9> 

CALL TITLE! ' ' 1 

S 'DISTANCE DOUNUIND FRON TURBINE , X^RkLHl .fDtt' , 
t lee.'UIND SPEED AT UAKE CENTER, tZlS',ie«, 1«. ,5. ) 
CALL CRAPH{«.,XHAX^l«.,*.,VnAX/S. J 

CALL HEADINC '(BS*). NORNALIZED UIND SPEEDS', 109,1., 2) 
CALL HEADINC'AT THE UAKE CENTERS ',199,1., 2) 

C 

C . . . READ FIRST CARD 
C 

READ(12,19,END>39) XUAL,YUAL 
19 FORNATiaFie.'O 

XPOS* (19. /XNAX )SXUALSXSCALE 
YPOS> ( 5 . /YHAX ISVUALXVSCALE 
CALL STRTPT(XPOS,VPOS) 

C 

C . . . PLOT ENTIRE DATA SET 



ooo ooo ooo o o n o 
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c 

aa CONTINUE 

XPOS • ( la . /XNAX XXUALtXSCALE 
VPOS • < S . /VNAX )«VUAt« VSCALE 
CALL CONNPTI XPOS, VPOS) 
READ(t2.ia,EN0«3a) XUAL,VUAL 
00 TO sa 
3a CONTINUE 

. . . NAKE A HARO COPY THEN ERASE SCREEN 

CALL lOUAIT 
CALL HDCOPV 

. . . CLOSE DATA SET 

CL0SE(UNIT-I2) 

. . . TERNINATE PLOT 

CALL ENDPLta) 

. . . RETURN TO NAIN ROUTINE 

RETURN 

END 
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c 

c 

c 

c . 
c 


c 


5 

C 


6 

C 


C 

C . 

c 

c 

c . 
c 


c 

c . 

c 


c 

c . 
c 

le 


SUBROUTINE XUSPR 
COfinON /PSCRCE/RR«,OH« 


. . SET HAXinun VALUES FOR X AND V AXIS 

XNAX>59. 

VHAX*!. 

1F(.N0T.(RRR .OT. 3)) GO TO S 
XHAX'SO.XRRR 
CONTINUE 

IF(.N0T.(UH9 .GT. 1 )) GO TO 6 
yflAX-l.»UH« 

CONTINUE 

XSCALE-XNAX/50. 

VSCALE-YNAX/1. 

. . OPEN DATA SET 

OPEN (UNIT-1‘4. FILE-'XUSPR. DAT'. ACCESS-'SEOIN' > 

. . RE-INITIALIZE PLOTTER 

CALL BGNPLCl) 

CALL PHVSOR(l..l. ) 

CALL PAGEUP. .11. ) 

. . PLOT AXIS 

CALL HEIGHT(0.3) 

CALL CONPLX 

CALL nXlALFC 'STAND'. '«' ) 

CALL nX2ALF( 'L/CSTD'. '•' ) 

CALL HX3ALF( 'NATHENATIC'. 'X' ) 

CALL nX4ALF( 'INSTR'. 'I' ) 

CALL TITLE ( ' ' 1 

t 'DISTANCE DOUNUIND FRON TURBINE . X/RUH1.ID»S', 
t 100. 'POWER FACTOR. P/PILHI .XO«'.100.10. .5. ) 

CALL GRAPH(O..xnAX/lO..0..yNAX/5.) 

CALL HEADINC '(B4»). TURBINE POUER FACTORS '.100.1.. 2) 
CALL HEADIN('FOR AN UNBOUNDED UAKES',1B0.1..2) 

. . READ FIRST CARD 

READCM. 10. END-30) XUAL.VUAL 
FORNAT(3F10.4) 

XPOS-(10./XnAX)tXUALSXSCALE 
VPOS • ( 5 . /VHAX ) t WALtVSC ALE 
CALL STRTPT(XPOS.VPOS) 


C 

C 


0 0 0 


PLOT ENTIRE DATA SET 



ooo f> r> o o f> o o o o o 
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c 

29 CONTINUE 

XP0S*(l«./XI1AX)<XUA(.<XSCAt.E 
VPOS* (S./VNAX )SVUALtVSCALE 
CALL CONNPT(XPOS.VPOS) 
READ(14*10.END>3«) XVAL.VUAL 
00 TO 20 

30 CONTINUE 

. . . riAKC HARD copy THEN £<’ ■$£ SCREEN 

CALL lOUAIT 
CALL HDCOPV 

. . . CLOSE DATA SET 

CL0SE(UNIT-14> 

. . . TERNINATE PLOT 

CALL ENOPL(0> 

. . . RETURN TO HAIN ROUTINE 

RETURN 

END 
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SUBROUTINE VUSUIR( lUNIT.UFILE.XNPTDH. I ) 

C 

COftnON ^PSCALE^RRR.UM 
C 

DOUBLE PRECISION UFILE 
C 

REAL XNPTDH(2«) 

C 

INTEGER UORK(13),ITirLE(t3),PLTNUN 
C 

LOGICAL NEULIN 
C 

DATA FLAO>’-99.>’ 

... SET PLOT HUflBER FOR TITLE ANNOTATION 
PLTNUfl-I+4 

... SET NAXinUN UALUES FOR X AND V AXIS 

XNAX*6. 

VNAX-1.2 

IF(.NOT,(RRR .GT. 2. )> CO TO 1 
XNAX'S.BRRR 
CONTINUE 

IF(.HOT.(UH0 .GT. 1.)) GO TO 2 
VNAX»1 .2*UH0 
CONTINUE 

XSCALE-XNAX/6. 

VSCALE*VnAX/'I.2 

. . . OPEN DATA SET 

OPENIUNIT* lUNIT.FILE-UFILE, ACCESS* 'SEQIN') 

. . . 'CREATE’ TITLE FOR PLOT 

ENCOOe<65,5,UORK) PLTNUN.XNPTDNC I ) 
F0RfMT('(',I2, LATERAL UINO SPEED PROFILE 
« 'AT X/RELH1.0DILXHXB *',F7.e, '»') 

DECODE (6S. 6, WORK) ITITLE 
FORNATdSAS) 

. . .RE-INITIALIZE PLOTTER 

CALL BCNPL(i) 

CALL PHVSORa..!.) 

CALL PAGE(9.,9. ) 

. . . PLOT AXIS 
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c 

c . 

c 

19 

C 

C . 
C . 

c . 
c 

15 

C 


29 

C 


39 

C 


35 

49 

C 

C 


C 

59 

C 

C . 
C 

C 

C . 
C 

129 

C 

C . 
C 


CALL HEIGHT ( 9 . 2 ) 

CALL COfIPLX 

CALL nXlALFC 'STAND', 't') 

CALL nX2ALF('L/CSTO 

CALL NX3ALF('nATHEnATlC'.'X' ) 

CALL HX4ALF('INSTR'.'l') 

CALL TITLEC' '.1, 'LATERAL COORDINATE , V/R9LH1 . 90 **', 199, 
'UIND SPEED, U/UiLHl.»0«',199,6.,6. ) 

CALL GRAPHt9.,XNAX/6.,*.,VnAX/6. ) 

CALL H£ADIN(ITITLE(1),35,1.,2) 

CALL HEADIN(ITITLE(8),199,1.,2) 

READ FIRST RECORD 

READ(IUNIT,19,EN0«59) VUAL,UBAR 
F0RflAT(2F19.4) 


PLOT DATA SET 


CONTINUE 

IF( .NOT. (VUAL .EO. FLAG)) GO TO 29 
NEULIN-.TRUE. 

CONTINUE 

IF(.NOT. NEULIN) CO TO 39 

READ< lUNIT, 19,EHD-S9 JYVAL.UBAR 
XPOS* VUAL* ( 6 . /XNAX )*XSCALE 
YP0S-UBAR*(6./VflAX)»VSCALE 
CALL STRTPT<XPOS,VPOS) 

NEUL1N-. FALSE. 

CONTINUE 

IF(.NOT.(.NOT. NEULIN)) GO TO 49 
XPOS • VUAL* < 6 . /XNAX )*XSCALE 
YP0S-UBAR»(6./VnAX)*VSCALE 
IF< .NOT. (XPOS .LT. 19.)) GO TO 35 
CALL CONNPT(XPOS,YPOS) 
CONTINUE 

CONTINUE 


REA0(IUNIT,19,EN0<59)YUAL.UBAR 
GO TO 15 

CONTINUE 

CLOSE FILES 

CLOSE(UNIT-IUHIT,DISPOSE-'SAUE' ) 
HAKE A HARD COPY AND ERASE SCREEN 
CONTINUE 

RETURN TO MAIN LINE 


RETURN 

END 
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t 

c 

c 

c 

c 

c 

c 

c . 
c 

c 

c . 
c 

c 


s 

c 


6 

c 


c 

c . 
c 


c 

c . 

c 

7 

c 

8 

c 

c . 

c 


c 

c . 
c 

t 


SUBROUTINE UUSZ( lUNlT.FIlE. AH, J) 

COnnON ✓PSCALE/RRR.UHB 
REAL XNPToncaa) 

INTEGER UORKdai.lTITLEdai.PLTNUH 
LOGICAL NEULIN 
DATA FLAOZ-SO./ 

. . SET PLOT NUMBER FOR TITU 
PLTNUM-4*J-H 

. . SET MAXIMUM VALUES FOR X AND V AXIS 

XMAX*1.2 

VMAX-S. 

IF(.N0T.(RRR .GT. 2)) CO TO 5 
VMAX-S. »RRR 
CONTINUE 

IFC.NOT.fUHO .GT. 1.)) GO TO 6 
XMAX-1.2XVHO 
CONTINUE 

XSCALE-XMAXZl.2 

VSCALE-VMAXzS. 

. . OPEN. DATA SET 

OPEN<UNIT-IUNIT, FILE-FILE, ACCESS-'SEQIN', 
t DEVICE- 'DSX' ) 

. . PRODUCE TITLE 

ENCODE(60,7,UORIC) PLTNUM 

FORr'rtTJ ' ( 12. ' >. VERTICAL UIHD SPEED PROFILES 
t 'AT UA<E CENTERI') 

OECOOE(6e,8,UORK)ITITLE 
FORMAT (12A5) 

. .RE-INITIALIZE PLOTTER 

CALL BGNPLd) 

CALL PHVSORd.,1. > 

CALL PAGE(8.,8. ) 

. . PLOT AXIS 

CALL HEICHT(«.3> 
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CALL COflPLX 

CALL (1X1ALP<'STAMD','«M 
CALL nX2ALF<'L/CSTD','«') 

CALL nX3ALF<'HATH£nATIC'.'«') 

CALL f1X4ALF< 'INSTR', 'k' ) 

CALL TITLE<' M.'UINO SPEED, U/UILHl 
t 'ALTITUDE, 2/RkLHl.«D»S',ie0,6.,5. ) 
CALL aRAPH(«.,XnAX/6.,0.,vnAX/$.l 
CALL HEADIN(ITITLE(1),35,1.,2) 

CALL HEA0IN(ITITLE(8),ie«,l.,2> 

C 

C . . . READ FIRST RECORD 
C 

READCIUNIT, ie,END*Se) UUAL,ZUAL 
le FORHAT<2Fl».<) 

C 

C . . . 

C . . . PLOT DATA SET 

C . . . 

C 

15 CONTINUE 

C 

IF(.NOT.(UUAL .EQ. FLAG)) GO TO 20 
NEULIN*.TRUE. 

20 CONTINUE 

C 

IF(.NOT. NEULIN) GO TO 30 

READ! lUNIT, 10, END-50 )UUAL,ZUAL 
XPOS-UUAL»(S./XMAX)»XSCALE 
VPOS - ZUAL* ( 5 . /VriAX )»VSCALE 
CALL STRTPT(XPOS.VPOS) 

NEULIN-. FALSE. 

30 CONTINUE 

C 

IF( .NOT. ( .NOT. NEULIN)) GO TO 40 
XPOS-UUAL»(6./XriAX)»XSCALE 
VPOS-ZUAL* ( 5 ./VflAX )»YSCALE 
IF( .NOT.CXPOS .LT. 10.)) GO TO 35 
CALL CONNPT(XPOS,VPOS) 

35 CONTINUE 

40 CONTINUE 

C 
C 

READdUNIT, 10, END-50 )UUAL, ZUAL 
GO TO 15 
C 

50 CONTINUE 

C 

C . . . ORAU A STRAIGHT LINE AT Z-AH 
C 

XPOS-0. 

VPOS- ( 5 . / VrWX ) t AHtVSCALE 
CALL STRTPT(XPOS,VPOS) 

XPOS-6. 

CALL CO«NPT<XPOS,YPOS) 



non non o o o o ooo non ooo 
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. . . DRMJ A DASH AT 2-AH^J 

V^-ts.^VfWJOttAH-MJIVSCAlE 
CALL STRTPT(XPOS.VPOS) 
XPOS'4.5 

CALL CONHPT(XPOS.VPOS) 

. . . DRAU A DASH AT Z*AH-l 


Woi • ( 5 . /VHAX ) • ( AH- 1 )»VSCALE 
CALL STRTPT(XPOS.VPOS) 
XPOS'9.5 

CALL CONNPTCXPOS.YPOS) 


. . . CLOSE FILES 

CLOSE(UMIT-IUNIT) 

. . HAKE A HARD COPY AND ERASE SCREEN 

CALL lOUAIT 
CALL HDCOPY 

. . ; TERHINATE PLOT 

CALL ENDPL(e) 

. . . RETURN TO HAINLINE 


RETURN 

END 
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Appendix C 

WAKE CALCULATIONS ON THE TI-59 CALCULATOR 

This appendix contains a description of a version of the 
wake model for the Texas Instruments TI-59 programmable 
calculator. A Texas Instruments PC-IOOC print cradle is used 
with the calculator. Mathematically, the model is identical 
to the model presented in this report and impletmented as a 
FORTRAN computer program. There are some limitations 
associated with the calculator version of the model. First, 
no plots are produced. All output is in numerical form and 
printed on the PC-IOOC printer. Second, all units of length 
are normalized by the rotor radius, and all units of wind 
speed are normalized by the free stream wind speed. There 
are no options for output in physical units. Third, atmospheric 
turbulence must be input as the standard deviation of wind 
direction. There is no option for input as a Pasquill class. 

There are two programs described in this appendix. The 
first program calculates values of downwind distance, x, 
wake radius, r 2 > normalized wind speed deficit at the wake 
center, Au^, normalized wind speed at the wake center, u^, and 
turbine power factor, P. The second program generates numbers 
for wind speed profiles for Region III. 

It is assumed that the user of these programs is familiar 
with the TI-59 programmable calculator and is reasonably pro- 
ficient in programming the calculator. The TI-59 calculator 
was chosen because it is widely used in scientific institutions. 
The wake program will not fit on the TI-58 calculator. The 
program for wind speed profiles will fit on the TI-58. Program 
listings are given on the following pages. 


Wake Program 

For the wake program, the calculator must be properly 
partitioned. Input of 3 0P17 will properly partition the 
calculator with 720 program steps and 30 memory locations. 

The inputs to the program are the initial wind speed 
ratio, m, and the standard deviation of ambient wind direction, 
a„ . To run the program, enter m, and press A. The input 
value of m will be printed. When 0.02 appears on the display, 
enter OQ,and press B. The value of Oq will be printed. A 
value of Oq = -1000 should be used to^ obtain the Abramovich 
solution (i.e., to make da/dx = 0 in equation (3-20)). Output 
consists of sets of five numbers: downwind location, x; wake 

radius, r 2 ; normalized wind speed deficit at the wake center. 
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Au ; normalized wind speed at the wake center, u ; and turbine 
power factor, P. The first set of numbers is for the initial 
wake (i.e., x = 0) . The second set of numbers is for the 
end of Region I (i.e., x = x„) . The third set of numbers is 
for the end of Region II (i.e., x = Xj^) . The following sets of 
numbers are printed during- the numerical integration in Region 
III. Since the turbine power factor is undefined in Regions I 
and II, it does not appear with the first three sets of 
numbers. Table C-1 shows sample output. The program terminates 
if x>60. The user may terminate the program at any time by 
pressing R/S. 

The parameters stored in the 30 memory locations are shown 
in Table C-2. This information is not necessary for running 
the program but is given for users who may desire to modify 
the program. A TI-59 listing of the program is given below. 

The description of the program is given adjacent to the program 
listing . 


Input m 

000 
00 i 
002 
003 

76 LBL 
1 1 fl 
42 STO 
02 02 

Print m 

LI U 4 

005 

93 RDV 
99 PRT 

Set value of x at 

006 

6 t* 

which integration 

IJ IJ 7 
003 

U 0 0 

of Region III ends 

•-»i_ I* * 4 1 

O *“» 

For numerical 

0 1 0 

—i 

m: ' « 

00 0 

integration in 

0 1 } 

02 2 

Region III, set 

0 1 2 

42 STD 

Ar = 0.02 

013 

1 9 1 9 

Input 0g 

0 1 
0 1 5 

9 1 F-! X S 
76 LBL 

0 1 6 

12 B 
99 PRT 

Print Og 

ij 1 • 



TiTTT 

ET •-! 
. » .*1 


0 1 9 

24 


020 

65 


0 2 1 

93 


022 

00 

Calculate a 

LI ii- 

03 

from equations 

024 

54 

(3-20) and (3-28) 

025 

li— iL. 

and store in 

026 


location 01 

027 

65 


023 

93 


029 

00 


030 

03 


031 01 1 

032 65 X 

033 01 1 

034 93 . 

035 03 9 

036 09 7 

037 95 = 

033 42 STD 

039 01 01 

TJWrr 93 RDV 

041 00 0 

042 99 PRT 


Print X = 0 
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043 

53 

1. 



044 

53 

if 



045 

43 

RCL 



046 

02 

02 

Calculate 

r^i from 
(3-18) 

047 

85 

+ 

equation 

048 

01 

1 

and store 

in 

049 

54 


location 

03 

050 

55 

H- 



051 

02 




052 

54 

*1 



053 





054 

95 

= 



055 

A •-» 

STD 


Print r 


0- 


Print Au, 


n 



U b U 

43 

RCL 

Print u 

061 

02 

02 

c 

062 

35 

1/K 


0 6 3 


PRT 


U64 

25 

U L K 


0 6 5 

53 

if 


0 6 6 

93 

• 


067 

02 



068 

01 

1 


069 

04 

4 


070 

85 

- 1 - 

Calculate 

071 

93 

• 

/O. 214+0. 144m 

l”! p" 

0 1 

1 


073 

04 

4 

and store in 

1974 

04 

4 

location 05 

1975 

65 



076 

43 

RCL 


077 

02 

02 


078 

54 

} 


079 

34 

TK 


080 

95 

= 


081 

42 

STD 


0 8 2 

05 

195 


_ _ _ 




U c* o 

43 

k L L 


0 y *4 

03 

03 

Calculate 

085 

55 

-r 

from equation 

0 8 b 

43 

RCL 

(3-34) and store 

U87 

05 

05 

in location 04 

U88 

95 



089 

42 

STD 


090 

04 

04 


• Jb 


57 


0 3 
99 

01 
QQ 


03 

PRT 

1 

PRT 



U9i 


RCL 


0 b >■ 

i9j]: 

0 3 


093 

65 



094 

5 y 

if 


095 

01 

1 


l9 b K 

85 

+ 


097 

43 

RCL 


0 9 8 

02 

02 


0 9 9 

54 

’j 


1 00 

55 



101 

53 

if 


102 

93 

• 

Calculate (x„)„ 
from equation 

103 

02 

im. 

104 

07 

r‘ 

(3-35) and store 

105 

65 


in location 06 

106 

53 

if 


lO";^ 

43 

RCL 


108 

02 

02 


109 

“7 C! 
1 -J 

- 


1 10 

01 

1 


1 1 1 

54 



112 

65 



1 1 

43 

RCL 


114 

05 

05 


115 

54 

*;i 


116 

95 

= 


117 

42 

STD 


118 

06 

06 
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119 

43 

RCL 


120 

04 

04 


121 

55 



122 

53 



123 

53 

< 


1 ii!4 

43 

RCL 


125 

04 

04 


126 

55 



127 

43 

RCL 


128 

06 

06 


129 

54 

) 


130 

•— * •«’ 

M 

t 1 

Calculate Xpj 

131 

0 • J 

+ 

from equation 
(3-39) and store 

132 

133 

53 

02 


in location 06 

134 

65 



135 

43 

RCL 


136 

01 

01 


137 

138 

139 

54 

) 


54 



140 

34 

fK 


141 

95 

= 


142 

42 

STO 


143 

06 

06 

Print Xpj 

■ 144. 


RDV 

145 

99 

PRT 

Print T 22 

1 46 
147 

43 

04 

RCL 

04 


148 

99 

PRT 

Print Au 

c 

149 

01 

1 

150 

99 

PRT 


151 

43 

RCL 

Print u^ = 1/m 

1 52 

153 

02 

35 

02 

1/^K 


154 

99 

PRT 


155 

25 

CLR 


Calculate 

/0.134+.124m 

and store in 
location 07 


Calculate n 
from equation 
(3-48) and store 
in location 08 


156 

53 

r' 

157 

93 

• 

158 

01 

1 

159 

03 

.J: 

1 60 

04 

4 

161 

y 5 

+ 

1 62 

99 

• 

163 

01 

1 

164 

02 

k. 

165 

04 

4 

1 bb 

65 


167 

H ‘“i 

RCL 

168 

02 

02 

169 

54 


170 

34 

fK 

171 

95 

= 

172 

42 

STD 

173 

07 

07 

T7T- 

43 

RCL 

175 

05 

05 

176 

65 

sy 

177 

53 

( 

178 

01 

1 

179 

75 

- 

180 

43 

RCL 

181 

07 

07 

182 

54 

J 

183 

55 


184 

5 • 

r’ 

185 

tr -i 

. ."i 

r; 

186 

01 

1 

187 

75 

- 

188 

43 

RCL 

189 

05 

05 

190 

54 

) 

191 

65 

X 

192 

43 

RCL 

193 

07 

07 

194 

54 

> 

195 

95 

= 

196 

42 

STO 

197 

08 

08 
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Calculate Xj^ 
from equation 
(3-47) and store 
in location 09 


Print X, 


Calcualte 22 
from equation 
(3-50) and store 
in location 10 


T78 

43 

RCL 

199 

06 

06 

200 

65 

V * 

201 

**1 

RCL 

202 

OS 

08 

203 

95 

= 

204 

42 

STD 

205 

09 

09 

206 

98 

fiDV 

207 

99 

PRT 

208 

43 

RCL 

209 

133 : 

03 


Print r^^ 

Store r 22 in 
locations 11 and 
14 as the initial 
radius for R. Ill 

Print Au„ 


Print u = 1/m 
c 


Calculate Au fronf;f'' 
Subroutine and 
store in location 
12 23> 

Calculate (dr/dx)j^ 4 i 
from Subroutine 24 

B’ and store in 24 l 

location 13 


85 + 

4 3 R U L 

05 08 

6 5 X 

5 3 ( 

43 RCL 
04 04 

75 - 

43 RCL 
03 03 

54 ) 

95 = 

42 STD 
1 0 1 0 
99 PRT 
42 STD 
11 11 

42 STD 

14 14 

01 1 
99 PRT 

43 RCL 
02 02 
35 1/K 
99 PRT 
25 CLR 
71 SBR 

16 Fi • 
42 STD 
12 12 
71 SBR 

17 B* 

42 STD 
13 13 


Label D 

24^. 

247 

Calculate r' and 


Store r’ in " 
location 14 for 
Subroutine A* 
Call Subroutine 

A* for Au^* 

Call Subroutine 
B* for Cdr/dx)„' 
Store (dr/dx)^*"^ 
in location 17 


Calculate 
(dr/dx) 
and store in 
location 18 


49 

43 

RCL 

50 

19 

19 

c < 
• J i 

95 

= 

,J ^ 

42 

STD 

53 

15 

15 

54 

66 

f H U 

55 

42 

STD 

56 

14 

14 


71 

SBR 

58 

1 s 

FI ■ 

"5F 

71 

SBR 

60 

1 "7 

4. { 

B ■ 

b 1 

42 

STD 

62 

1 r 

A 1 

17 

hi"; 

53 

if 

64 

53 

if 

65 

53 

if 

66 

43 

RCL 

67 

13 

13 

6 y 

85 

+ 

6 9 

43 

RCL 

70 

17 

17 

"7 

i i 

54 

’"l 


55 

02 

54 

-T- 

( 0 

7 4 


”? cr 

r -j 

.."1 


7 6 

85 

+ 

»'* 

43 

RCL 

1 

01 

01 

79 


t '1 

so 

54 

> 

8 1 

y 4 

rx 

0 2 

95 

= 

0 

42 

STD 

84 

18 

18 
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o 3 1 5 • C 

332 01 1 




285 

43 

RCL 





286 

09 

09 






85 

+ 



Calculate x' 

31. O 

43 

RCL 



from equation 

289 

19 

19 



(3-62) 

and store 

29fi 

=i5 


Calculate 

first 

in location 09 

99 1 

43 

RCL 

coefficient of 



29? 

IS 

18 

equation 

(3-78) 



Q 

95 

— 

and store 

in 



294 

42 

STO 

location 

21 



295 

09 

09 





290 


flliV 



Print 

X’ 

ill 9 7 

QO 

PRT 





ill 9 

43 

RCL 



Print 

r ' 

M 

15 

15 





300 

99 

PRT 



Print 

Au ' 

30 1 

43 

RCL 




c 

3U2 

1 b 

. 1 6 





303 

99 

PRT 





3 0 4 

oi 

1 




U3 3 
65 X 
336 43 RCL 


335 65 


339 03 3 

340 65 X 

341 43 RCL 

342 20 20 


344 

345 
34 


< j 

43 RCL 
20 20 


'“I H -T" H C* 

.!« ^ r 4 .j 
340 03 

•~i A c* .-I 

4 7 . J “T 

350 55 

351 02 

352 95 


I 


42 STO 


3UD I'D 



306 

43 

RCL 


307 

16 

1 6 

Calculate u ' 

308 

65 


and print 

3 0 9 

53 

1 


3 1 0 

01 


31 1 

75 

- 


312 

43 

RCL 


313 

02 

02 


3 1 4 

35 

1/K 


315 

54 

*) 


3 1 6 

95 

= 


317 

99 

43 

PRT 

RCL 


1 b 


3 1 9 

1 b 

16 


320 

65 



:-i 1 

C “I 
.1.1 

l'* 

Calculate B 

3 2 2 

01 

1 

from equation 

8 2 2 

75 

- 

(3-75) and store 

324 

43 

RCL 

in location 20 

325 

02 

02 


326 

35 

ixX 


•-I 

O ill i 

54 

y 


!-I 

97i 

= 


Pi 2 9 

42 

STO 


3 3 0 

20 

20 




354 

21 

1 



555 

53 

i;* 



3 5 6 

06 

6 



• J ! 

65 

X 



350 

43 

RCL 



359 

20 

20 



3 6 0 

-7 cr 
i J 

- 



3 6 1 

01 

1 



3 6 2 

02 

•-1 



■-» D 

65 

X 



3 b 4 

•■f 

RCL 

Calculate 

second 

3 6 5 

20 

20 

coefficient of 

3 6 6 

O -I* 

X iZ 

equation 

(3-78) 

3 6 7 

o 5 

+ 

and store 

in 

3 6 3 

06 

b 

location 

22 

3 6 9 

6.5 

X 



370 
3 i 1 

43 

'Z* i”l 

RCL 



^ !_l 

45 

fc-U 

YX 



•-» r* -z* 

03 

3 



■I* f H' 

54 

;> 



“7 cr 
1 ’J 

55 

-f- 



3 7 6 

03 

o 



O ( (' 

93 

■ 



O i C* 

05 

.5 



379 

95 

= 



380 

42 

STO 



381 

ill 3 

22 
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;1| C; 

03 

,"i 



384 

94 

+ ..••• - 



•Z * «“» cr 

.1 II . 1 

65 




386 

43 

RCL 



‘ i* 

20 

20 



►Z • C* 

y 5 

+ 



389 

01 

1 



390 

08 

C* 

1^1 



391 

65 




392 

43 

RCL 



393 

20 

20 



•I* y A- 

•“t *”i 
•2« •!« 

* 

1***1 ^ 



395 

”7 cr 

r -J 

- 

Calculate 

third 

Q h:. 

01 

1 

coefficient of 

■-» y 

05 

5 

equation 

C3-78) 

•I' o 

65 

RCL 

and store 

in 

: 9 9 

43 

location 

23 

4 0 0 

20 

20 



401 

45 

Y ■■■■• 



402 

0 3 

3 



403 

54 

J 



404 

55 




405 

05 

5 



406 

95 

= 



407 

42 

STD 



408 

23 

23 



4 u y 

53 




410 

01 

1 



41 1 

02 

1^ 



412 

94 

+ - 



413 

65 




414 

43' 

RCL 

Calculate 

fourth'^ ] z 

20 

20 

coefficient of 

i b 



equation 
and store 

(3-78) 

in 

4 1 i' 
418 

O ■ J 
02 

+ 

location 

24 

419 

00 

u 



420 

65 




421 

43 

RCL 



422 

20 

20 



423 

45 

YX 



424 

03 

•T' 

•i.' 



425 

54 

> 



426 

55 

4- 



427 

06 

6 



428 

93 




429 

05 

5 



430 

95 

= 



431 

42 

STD 



432 

24 

24 



433 

53 '1 


434 

03 3 


435 

65 X 


4 3 6 

43 RCL 


437 

.•i •** o 

20 20 


*T O 

439 

•Z* O i*'*i 

75 - 

Calculate fifth 

440 

0 1 1 

coefficient of 

441 

05 5 

equation (3-78) 

442 

65 X 

and store in 

443 

43 RCL 

location 25 

444 

20 20 


445 

45 yx 


446 

0 3 3 


447 

• J “i* 


448 

55 


449 

08 8 


450 

95 


451 

42 STD 


452 

25 25 


453 

0 6 6 


454 

65 X 


.-I c*c- 

4 . J . J 

43 RCL 


456 

20 20 

Calculate sixth 

457 

45 y:^< 

coefficient of 

*4 5 y 

03 3 

equation (3-78) 

459 

r cr • 

•J ” 

and store in 

460 

j~! 

location 26 

461 

93 


462 

05 5 


463 

95 = 


464 

42 STD 


465 

26 26 


4 6* 6* 

01 1 


467 

94 +x- 


468 

65 X 


469 

43 RCL 

Calculate 

470 

20 20 

seventh 

471 

45 

coefficient of 

472 

03 3 

equation (3-78) 

473 

55 

and store in 

474 

01 1 

location 27 

475 

01 1 


476 

95 = 


477 

42 STO 


478 

i- ( £l 1* 
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479 

53 

( 


4S0 

53 



481 

43 

RCL 


488 

02 

02 


4 8 3 

85 

+ 


484 

01 

1 

Calculate 

4:=:'=: 

5U 


from equation 


55 


(3-72) and store 

UR? 

CT 


in location 28 

483 

02 



489 

65 

V 


490 

43 

RCL 


491 

02 

02 


492 

54 

) 


493 

54 

^1 


494 

34 

fK 


495 

95 

= 


496 

42 

STO 


497 

iiL o 



49 y 

24 

CE 

Calculate r_/r- 

499 

55 

-r 

and store 

500 

43 

RCL 

in location 29 

501 

15 

15 


502 

95 

= 


503 

42 

STO 


504 

29 

29 


505 

02 

w 


506 

65 



507 

53 

1, 


508 

43 

RCL 


509 

21 

21 


510 

85 

+ 


511 

43 

RCL 


512 

I— 


Calculate F from 

5 1 3 

65 

X 

equation C3-78) 

514 

43 

RCL 


515 

29 

vQ 


516 

45 

$ $ 
1 


517 

01 

1 


518 

93 

• 


519 

05 

cr 


520 

35 

+ 



521 

43 

RCL 


522 


o •-* 


•J vl I* 

524 

b‘j 

43 

RCL 


525 

ka -* 

29 


526 

45 

Y 


527 

03 

.^1 


523 

8 5 

+ 


529 

43 

RCL 


530 

24 

24 


531 

65 

X 


532 

43 

RCL 


533 


29 


534 

45 

V 


535 

04 

4 


5 3 6 

93 

■ 


5 o 7 

05 

5 


538 

35 

+ 


539 

43 

R ij L 


540 

25 

25 


541 

65 

V* 

Calculate P 

542 

43 

RCL 

from equation 

543 

29 

29 

(3-78) 

544 

45 

YX 

(concluded) 

545 

06 

6 


546 

Pi 5 

+ 


547 

43 

RCL 


548 

2 6 

26 


549 

65 

X 


550 

4c> 

RCL 


551 

29 

k— m' 


552 

45 

yx 


553 

' 07 

“7 

1* 


554 

93 

■ 


555 

of: 

Fi 


556 

85 

+ 


557 

43 

RCL 


553 


•“*"7 
w 1 


559 

65 

s..* 


560 

43 

RCL 


561 

29 

29 


562 

45 

1 • 
1 


563 

09 

9 


564 

54 

;i 


565 

95 

= 


3 F 5 . 99 PRT 


Print F 
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^•6? 

43 

RCL 


56S 

15 

15 

Set r = r* 

569 

42 

STD 


570 

1 1 

1 1 


5 1 

43 

R L L 

Set Au^ = Au^' 

572 

16 

1 6 

c c 

573 

A •"! 
*-► ^ 

STD 


574 

12 

12 


575 

43 

RCL 

Set (dr/dx)^ 

576 

17 

17 

in 

£i77 

42 

3TD 

= (dr/dx)j^' 

573 

13 

13 


579 

43 

RCL 

I£ X < 60, 

530 

09 

09 

GO TO D 

C I-I ^ 

■JO i 

532 

0 o 

u. 

1 i 

IHV 

GE 


533 

14 

n 

Stop 

534 

91 

R^S 


535 

b 

LBL 


536 

16 

H ■ 


537 

0 3 



533 

93 

■ 

Subroutine A' to 

589 

07 


calculate Au^ 

590 

0 3 

.^1 

by equation 

591 

65 


(3-58) 

592 

5--1 

1. 


593 

9 o 

S 


594 

02 

i:! 


=i35 

OFi 



596 

03 

8 


597 

65 



598 

43 

RCL 


599 

02 

02 


600 

55 

•r 


Subroutine A' 
(concluded) 


60 1 

■ J •••’ 

( 

6 0 2 

43 

RCL 

603 

02 

02 

604 

75 

- 

605 

01 

1 

6 0 6 

54 

*) 

607 

75 

- 

60S 

C* --I 

•J O' 

( 

609 

24 

CE 

6 1 0 

O' -O 

•. .* “■ 

61 1 

75 

- 

612 

93 

■ 

6 1 3 

05 

ET 

614 

03 


615 

06 

6 

6 1 6 

65 


617 

A •”i 

•-to 

RCL 

6 1 3 

03 

03 

619 

O 0* 

X - 

620 

55 


621 

53 

(* 

622 

43 

RCL 

623 

14 

14 

624 

O' o 


625 

65 


626 

53 

( 

627 

43 

RCL 

623 

02 

02 

629 

"n c 
1 -J 

- 

630 

01 

1 

6 3 1 

54 

) 

632 

54 


6 3 3 

54 


634 

34 

r 

i i*'i 

6 3 5 

54 


6 3 6 

95 

= 

637 

42 

STD 

t' O 

1 6 

1 6 

639 

Qb 

•* 

RTH 
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640 

76 LBL 



641 

17 B' 



642 

93 . 



643 

02 2 



644 

07 7 



645 

55 



646 

53 



647 

02 2 



6 4 8 

65 X 

Subroutine B' to 

649 

43 RCL 

calculate 


6 5 0 

02 02 

(dr/dx) 

by 

f.Fl ! 

5Fi -r 

equation 

C3-60) 

652 

53 < 



653 

53 



654 

43 RCL 




02 02 



656 

75 - 



657 

01 1 



658 

54 > 



659 

65 X 



6* 6 Lf 

43 RCL 



661 

16 16 



662 

54 > 



6 6 3 

75 - 



6 6 4 

01 1 



665 

54 > 



6 6> 6 

95 = 



66*7* 

92 RTH 
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Table C-1 

OUTPUT OF WAKE PROGRAM FOR 
TI-59 CALCULATOR 


Description 

Output 

Sy mb o 1 

Input value o£ m 

■ 

m 

Input value of 

u 

10. 

^0 


0. 

x=0 

Initial wake 

1. 224744S71 


1 . 



0. 5 

u "" 
c 

Wake at end of 

8. 15901383 


1. 728597064 


Region I 

1 . 



0. 5 



12. 25580717 



1 . 9.3 1 590667 

r^- 

Wake at end of 

1. 

0 5 

2 2 

Region II 

Au 



u 



c 


12. 40933542 

X 

First point in 

ill K U IJ 1 5 *3 i J t' 7^ 

^2 

numerical integration 

. 9695134231 


in Region III 

. 5152432384 



. 2916912548 

P^ 


12. 5656277 

X 

Second point in 

2 . U 2 1 5 9 U b 6 7 

^2 

numerical integration 

. 9411535052 

#c 

in Region III 

. 5294207474 



. 3019249019 

p 


12. 72450863 


Third point in 

2. 041590667 


numerical integration 

. 91 4393 i i t. 3 
. 5423031 119 

in Region III 

^c 

.311 9 6 6 6 5 5 

F 
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Table C-2 

STORAGE LOCATIONS USED FOR 
WAKE PROGRAM FOR TI-59 CALCULATOR 


Storage 

location 

Symbol 

Definition 

01 

a 

Wake growth rate due to ambient 
turbulence 

02 

m 

Initial wind speed ratio, U^/Uq 

03 

""o 

Initial wake radius 

04 

^21 

Wake radius at the end of Region I 

05 

/O. 214+0. 144m 

06 

and 

Downwind extent of Region I from 


Abramovich solution and downwind 


extent of Region I 

07 


/O. 134+0. 124m 

08 

n 

X(,/Xh 

09 

X 

Downwind coordinate 

10 

^22 

^2 

Wake radius at end of Region II 

11 

Wake radius 

12 

Au 

c 

Wind speed deficit parameter 

13 

(dr/dx)^ 

Wake growth rate due to mechanical 
turbulence for wake radius, T 2 

14 

^2 

Wake radius used by Subroutine A' 

15 

^2 ' 

Incremented value of wake radius 

16 

Au^' 

Wind speed deficit parameter for 
wake radius, r' 

17 

Cdr/dx) 

Wake growth rate due to mechanical 


ill 

turbulence for wake radius, r2* 

18 

1 

Cdr/dx) g 

Effective growth rate of the wake 
radius 

1 

19 

Ar 

Increment in wake radius for integration 

20 

B 

Au^ (1-1/m) 

21 


First coefficient for equation (3-78) 

22 


Second coefficient for equation (3-78) 
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Table C-2 

STORAGE LOCATIONS USED FOR 
WAKE PROGRAM FOR TI-59 CALCULATOR (concluded) 


Storage 

location 


Definition 

23 


Third coefficient for equation (3-78) 

24 


Fourth coefficient for equation (3-78) 

25 


Fifth coefficient for equation (3-78) 

26 


Si:xth coefficient for equation (3-78)’ 

27 


Seventh coefficient for equation 
(3-78) 

28 


Radius of stream tube in the free 
stream flow 

29 


Ratio of r^ for the downwind 
turbine to the wake radius of 
the wake of the upwind turbine 
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Wind Speed Profiles 

The program for wind speed profiles generates the wind 
speed in the wake for the lateral wind speed profiles or for 

the vertical wind speed profile. The inputs to the program 

are the initial wind speed ratio, m; the height of the rotor 
hub Cin rotor radii), h^; the power law coefficient for the 
free stream wind speed profile, y; and the wake radius T 2 - 
The wake radius may be obtained for a given downwind location 
from the wake program. For horizontal wind speed profiles, 
the altitude at which the profile is to be calculated, z, is 
also input. The output of the program is the normalized wind 
speed in the wake at intervals of 0.1 rotor radii. The wind 

speed in the wake is normalized by the free stream wind speed 

at the hub altitude. An input value of y = 0 allows the 
wind speed to be normalized by the wind speed at the altitude 
of the profile. 

To operate the program: 

Enter m Press A 

Enter h Press R/S 

Enter y Press R/S 

Enter V 2 Press B 

For the lateral wind speed profiles: 

Enter z Press C 

For the vertical wind speed profile: 

Press E 

For lateral wind speed profiles at an additional altitude: 
Enter new z Press C 

To change X 2 ' 

Enter new T 2 Press B 

The program will terminate when the point on the profile 
is in the free stream wind. The last number printed is for 
the free stream wind. Table C-3 shows a sample output. Table 
C-4 shows the parameters stored in the memory locations. A TI-59 
listing of the program is given on the next page. The description 
of the program is given adjacent to the program listing. 



191 



UUI.I 

(' 6 

LBL 

Input m 

0 0 i 
002 

1 1 
42 

fl 

:3TD 


0 0 3 

01 

01 


004 

99 

PRT 

Print m 

0 0 5 

91 

R.--S 


on 4 

42 

STD 

Input 

007 

02 

02 


000 

Mb 

PRT 

Print h 

a 

009 

91 

R.-S 

Input Y 

010 
01 1 

42 

13 

STD 

13 

Print Y 

012 

9 9 

PRT 


0 1 3 

93 

a 

Set Ay=0.1 

014 

015 

01 

42 

1 

STO 


016 

03 

03 

Put 1 in 

U1 / 

01 

1 

t register 

0 1 3 
0 1 9 

9 1 

KIT 

R.--S 


020 

?^' 

LBL 

Input 

022 

12 

42 

B 

STO 


0 2 

04 

04 


024 

98 

flDV 

Print r- 

025 

98 

HDV 

L 

026 

99 

PRT 

Calculate rp 

02'^ 
023 
0 2 3 

53 

43 

( 

( 

RCL 

from equation 

030 

III 

01 

(3-18) and 

il 

o D 

+ 

store in 

032 

01 

1 

y 

location 05 

03:3 

54 


l”i -A 

55 



1 1 .’:I 

02 



C* 3 b 

54 

) 


!”| "y 

:34 

TK 


030 

95 

= 


0:39 

42 

STD 


040 

05 

05 


Calcuate Au 
from equati 
(3-58) 


g 


n 


041 

0:3 

‘Z* 

042 

9;3 

• 

|“l .4 



1 J . -| 

U i' 

i 

044 

03 

►-I 

045 

65 


U 4 o 

53 

( 

04"’ 

C" --1 

•-* O 


043 

9:3 

■ 

049 

02 

o 

050 

05 

5 

051 

08 

o 

ns'-' 

65 




053 

43 

RCL 

1-1 ET < 

O 4 

01 

IJ ,_i i-i- 

1 J X 

055 

c cr 
• J -J 


056 

cr 

.J .ji 

( 

057 

43 

RCL 

053 

01 

01 

059 

75 

- 

060 

01 

1 

0 6 1 

54 

j 

062 

54 



“? c* 


1,1 1« 

I' -J 


064 

cr -i 
•-* O 

l'' 

0 6 5 

24 

CE 

066 

O -Z* 

K^ 

06":^ 

75 

- 

063 

9:7 

B 

069 

05 

C 

.J 

070 

03 


071 

06 

6 

072 

65 


073 

43 

RCL 

074 

05 

05 

075 


K ^ 

076 

55 


077 

C •“» 

• J O 


073 

43 

R L- L 

079 

04 

04 

080 

:3 3 


031 

65 

*•••* 

082 

53 

i'.' 

083 

43 

RCL 

084 

01 

01 

085 

75 

- 

086 

01 

1 

087 

54 

J 

083 

54 

J 

089 

54 

J 

090 

34 

TK 

09 1 

54 

*) 

0 9 :3 

95 
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Store Au in 

033 

42 

STO 

location^06 

034 

06 

06 

Input z for lateral-' 3 j 

31 

Fl-'S 

wind speed profile 

036 

7 S 

LBL 


037 

13 

c 


0 3 3 

42 

STO 


033 

07 

07 


1 00 

38 

RDV 

Print z 

101 

3 8 

flDV 


102 

33 

PRT 


103 

00 

0 

Initialize y 

1 04 

42 

STO 


Jjl5 

03 

08 


106 76 LBL 


Calculate z from 
(3-85) 


equation 
store in 
09 


107 
03 

location j,-|g 

I i 0 

I I 1 


112 

113 

114 

115 

Calculate u from 
Subroutine A* and 117 
store in location 113 

10 113 

Calculate z^* from 120 
equation (3-96) and.21 
store in location l£2 
09 123 

124 

125 

1 ^ TZb 

Calclated u* from i27 
Subroutine A* and joc; 
store in location Ijlo.^ 

r~T~ IW 

Print y < i 


1 33 

Increase y by Ay 

135 
1 3 6 
137 
133 
133 


17 

B' 

43 

RCL 

02 

02 

75 

- 

43 

RCL 

07 

07 

35 

= 

42 

STD 

03 

09 

71 

SBR 

1 6 

FI* 

42 

STD 

10 

10 

43 

RCL 

02 

02 

85 

+ 

43 

RCL 

07 

07 

35 

= 

42 

STD 

03 

03 

71 

SBR 

1 6 

fi * 

42 

STD 

1 1 

1 1 

43 

RCL 

03 

08 

CiO 

yw 

RDV 

99 

PRT 

43 

RCL 

03 

03 

44 

SUM 

08 

08 


Calculate U 
from equation 
(3-98) and 
print 


Go to B’ if not 
in free stream 
flow 


Slop 

Vertical wind 
speed profile 


140 

43 

RCL 

141 

10 

10 

142 

65 


143 

43 

RCL 

144 

1 1 

1 1 

145 

65 

X 

146 

53 

if 

147 

43 

RCL 

143 

07 

07 

143 

55 

-r 

150 

43 

RCL 

151 

02 

02 

152 

54 

y 

153 

45 

V 

154 

43 

RCL 

155 

13 

13 

156 

35 

— 

157 

33 

PRT 


"T5S 43 RCL 
153 10 10 

160 65 X 

161 43 RCL 

162 11 11 

163 35 = 

164 22 INV 

165 67 EQ 

166 17 B* 

TSr 31 Rv’S 
~TSU 76 LBL 

163 15 E 

1 ,'i I on 0 

171 42 STD 


Initialize 7 


172 
1 73 
174 
-rnr 
176 


Calculate z 
from equation 
(3-85) and store 
in location 
09 


Calculate u from 
Subroutine A' 
and store in 
location 10 


08 
42 
07 
7 6 
10 


08 

STO 

07 

LBL 

E ■ 


177 

43 

RCL 

173 

02 

02 

173 

75 

- 

1 SO 

43 

RCL 

181 

07 

07 

182 

35 

— 

1 S3 

42 

STD 

184 

09 

09 

185 

71 

SBR 

1 36 

1 6 

fi * 

187 

42 

STD 

133 

10 

10 
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1 8 9 


RCL 





Calculate 

z 

* from^ 

02 

02 





equation 
and store 
location 

C3- 

in 

09 

96) 

191 

192 
1 9 3 

i“i cr 

C' *J 

43 

07 

+ 

RCL 

07 

Subroutine A' to 
calculate u 

o *■:> cj 
239 

7 6* 
16 
53 

LBL 
fi • 

r* 




194 

95 

— 


240 

43 

R C L 




195 

42 

STD 


241 

09 

09 




196 

09 

09 


242 

O 


Calculate 

11* 

from 

T97 

71 

SBR 

Calculate r/r- 

243 

85 

+ 

Subroutine A' and 
store in location 

198 

199 

16 

42 

fl ‘ 
STD 

from equation'^ 
(3-86) and store 

244 

245 

43 

08 

RCL 

08 

11 



200 

1 1 

1 1 

in location 12 

•“» . * 
il. C* 

O I' 

K ^ 



20 1 

43 

RCL 


*1* •** ”7 

54 

) 

Print z 



202 

07 

07 


'•“i .4 i“i 

ill ‘t C* 

34 

f7 




203 

98 

flnv 


249 

cr c- 
•J 

-r 




204 

99 

PRT 


250 

43 

RCL 




205 

43 

RCL 


251 

04 

04 




206 

10 

10 


252 

95 

— 




207 

65 



fi 3 

42 

STD 

Calculate 

U 

from 

208 

43 

RCL 


254 

12 

12 

equation 

(3- 

98) 

209 

1 1 

1 1 

T "P T* / y* 1 

255 


IHV 

and print 



2 1 0 
211 

65 

-1 

( 

X 3. X. j ^ ^ ^ J. y 

go to C' 

256 

257 

l‘* 1*’ 

IS 

GE 

c • 




212 

43 

RCL 

y =...I 

258 

01 

1 




213 

07 

07 

Return 

259 

92 

RTN 




214 

55 

H- 


260 

76 

LBL 




215 

43 

RCL 


261 

18 

c ■ 




2 1 6 

02 

02 


262 

01 

1 




217 

54 

*j 


263 

75 

- 




213 

45 

YX 


264 

43 

RCL 




219 

43 

RCL 


265 

06 

06 




220 

13 

13 

Calculate u 

266 

65 





i 

95 

= 

frm equation 

267 

53 

1 




ill il! 

99 

PRT 

(3-93) 

268 

01 




..*• , ~1 

43 

RCL 

269 

C 

1 

— 




224 

0 3 

03 


270 

43 

RCL 




“Z25 

44 

sun 


O'? 1 

1 X 

12 

12 




226 

07 

07 


•'*t 

45 

't 




.“i 1 

Im t 

43 

RCL 


.-I 

ii- 1 

oi 

1 




'O '1* C‘ 

w ll. 

10 

10 


274 

93 

■ 




229 

65 

X 


O ”7 cr 
w 1 • J 

05 

5 




230 

43 

RCL 


iC *• b 

54 

*) 

Go to E' 

if 

not 

231 

1 1 

11 


••"i 

C. f i ’ 

33 


in free stream 

•1* 'I* O 

95 

= 


•”i *”• 

b- 1 O 

95 


flow 



*1* *1* 

iL. 

INV 

Return 

2,-y 

92 

RTH 


ii. i-ii i n V 

234 67 EQ 

235 10 E‘ 


o b y 1 R y 


Stop 
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Table C-3 

OUTPUT OF WIND SPEED PROFILE PROGRAM FOR TI-59 


Description 



Cm cm Cn Cm C^< C'<! C 
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Table C-4 

STORAGE LOCATIONS USED FOR WIND SPEED 
PROFILE PROGRAM FOR TI-59 CALCULATOR 


Storage 

location 


Symbol 


Definition 


01 

02 

03 

04 

05 

06 

07 

08 

09 

10 
11 
12 


m 

Ay 

^2 

‘0 

Au^ 

z 

y 


z or z ’ 

V V 


u 


u’ 


r/r. 


Initial wind speed ratio 

Height of turbine rotor hub in 
rotor radii 

Interval between successive 
points in the profile 
calculation, in rotor radii 

Wake radius in rotor radii 

Initial wake radius in rotor 
radii 

Wind speed deficit parameter 

Altitude of wind speed profile 

Lateral coordinate 

Altitude of profile relative 
to hub of real or image 
turbine 

Wind speed in wake for wake of 
real turbine 

Wind speed in wake for wake of 
image turbine 

Radial coordinate in wake 
normalized by wake radius 


13 


Y 


Power law exponent for free 
stream wind speed profile 
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Appendix D 


ALTERNATIVE ALGORITHMS FOR CALCULATION OF 
THE DOWNWIND EXTENT OF REGION I 

In the development of the wake model, seven different 
approaches for the calculation of the downwind extent of 
Region I were considered. All of the approaches are based 
on the Abramovich approach for a jet in the absence of ambient 
turbulence. All of the approaches use the method of the 
square root of the sum of the squares for the adding of 
mechanical turbulence and ambient turbulence. For all of 
the approaches, the wake radius at the end of Region I is 
given by Abramovich equation (5.19) as 


21 


^21 ^ 

^0 ^d 


0 


/O. 214+0. 144m 


CD-I) 


This equation was given as equation (3-34) in the description 
of the analytic model. It is the result of a momentum balance 
between the initial wake and the end of Region I; therefore, 
it should be valid regardless of the presence or absence 
of ambient turbulence. 

* For mechanical turbulence, the downwind extent of Region 
I is given by Abramovich equation (5.20) as 


1+m (D-2) 

O.27(m-1)/0. 214 + 0. 144m 

This equation was given as equation (3-35) in the description 
of the analytic model. 


Definition of Approaches 

For the purposes of identification, let the seven 
approaches for calculating the downwind extent of Region I 
be defined as follows. 
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1. ri approach . - Add the effect of wake growth 
due to ambient turbulence to the erosion of 
the inner core as given by Abramovich, and 
define the end of Region I as the point at 
which the radius of the inner core becomes 
zero. 

2. V 2 approach . - Add the effect of wake growth 
due to ambient turbulence to the growth of 
the wake radius as given by Abramovich, and 
define the end of Region I as the point at 
which the radius of the wake becomes '^21’ 

3. b approach I . - Add the effect of wake growth- 
due to ambient turbulence to the boundary 
layer growth given by Abramovich, and define 
the end of Region I as the point at which 

the width of the boundary layer becomes ^^21’ 
The width of the boundary layer is 

4. b approach II . - Add twice the effect of wake 
growth due to ambient turbulence to the 
boundary layer growth given by Abramovich, 
and define the end of Region I as the point 
at which the width of the boundary layer 
becomes i’21* Twice the effect of wake 
growth due to ambient turbulence is added 
because ambient turbulence affects the 
boundary layer on both sides--at the potential 
core and at the wake boundary. 

5. Streamline ri approach . - This is the same 
as approach 1, except that erosion of the 
inner core is calculated relative to the 
streamline which passes through the initial 
wake radius, rather than relative to the 
line r = Vq as was done for the r^ approach. 

6. Streamline t2 approach . - This is the same 

as approach 2, except that growth of the wake 
radius is calculated relative to the stream- 
line which passes through the initial wake 
radius, rather than relative to the line 
r = Vq as was done for the T 2 approach, 

7. .5r7i approach . - This is the same as the 
first and second approaches, except that the 
erosion of the inner core and the growth of 
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the wake raditis are calculated relative to a 
line passing from the initial wake radius at 
X = 0 to r = .5r22 ^t the end of Region I, 

This line bisects the boundary layer. 
Mathematically, this approach is identical 
with approach 4. 

All of these approaches take an angle associated with 
the Abramovich solution for the initial region, augment 
the wake growth rate for a, the wake growth rate due to 
ambient turbulence Casing the square root of the sum of the 
squares for adding the turbulence components), and calculate 
the downwind e:xtent of Region I accordingly. The difference 
in the approaches is in the choice of angles from the 
Abramovich solution. The choice of angle for the various 
approaches is shown in Figure D-1. 


Description of Approaches 

By definition, the end of Region I is that point at 
which the potential core vanishes. Thus, the wind speed at 
the center of the wake at the end of Region I is Ug. A 
momentum balance between the initial wake and the end of 
Region I gives equation CD-I) . 


From the equation for boundary layer growth given by 
his equation C5.1), Abramovich gives the length of the 
initial region of the wake as equation (D-2). The m subscript 
denotes that the quantity is associated with mecahically- 
generated turbulence. Therefore, for mechanically-generated 
turbulence, the slope of the radius of the inner core is 
Ccf, Figure D-1) 


= ''^0 

dx U (X„)„ 


-rn/(x„) 


H^m 


(D-3) 


and the slope of the outer radius 


of the wake is 


^ ^2l'S ^ ^2l'^0 

dx /m CXpj)j^ 


(D-4) 


Equation CD-I) which gives the wake radius when the 
potential core has been completely eroded is based upon a 
momentum balance Ceased upon the fact that U = Ug at the 
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center of the walce at the end of Region I and the assumption 
of the wind speed profile given by equation (3-32)) and is 
therefore valid regardless of the presence or absence of 
ambient turbulence. Therefore, the presence of ambient 
turbulence only affects the distance, at which the 

end of Region I occurs. 

ri approach . - The first approach for calculating 
the value of aj. is the calculation of the effect of ambient 
turbulence on the erosion of the inner core and specifying 
the downwind extent of Region I as the point at which the 
radius of the inner core becomes zero. Under this definition 
the rate of erosion of the inner core is 





dx /m 


+a' 


(D-5) 


Since r^ decreases for rg at x = 0 to 0 at x = x^j, the downwind 
extent of Region I is given by 


X 


H 


(dr^^/dx) 


CD-6) 


where x^ is the downwind extent of Region I normalized by R^. 

r? approach . - The second approach for calculating x^^ 
consists of calculation of the effect of ambient turbulence 
on the growth of the outer radius of the wahe and specifying 
the downwind extent of Region I as the point at which the wake 
radius reaches the radius defined by a momentum deficit balance 
as given in equation (D-1). Under this definition. 


dr 


2 


dx 



dx /m 


+a‘ 




z 


(D-7) 


Since increases from rg at x = 0 to ^ the 

downwind extent of Region I is given by 


X 


H 



(dr2/dx) 


CD-8) 
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b approach I . - The third approach of calculating jCr, 
considers the growth of the boundary layer, applies 

the wake growth due to ambient turbulence to the boundary 
layer. For mechanical turbulence, the growth rate of the 
boundary layer is 


^ ^'21 j. 

iJm Cxj,)^ 

With ambient turbulence, the growth of the boundary layer 
is 

db 


dx 

Since the boundary layer grows from 0 at x = 0 to r 2 ^ at x = x„ 
the downwind extent of Region I is 


21 




+a‘ 


(D-10) 


X 


H 


^21 

db/dx 


CD-11) 


b approach 1 1 . - The fourth approach for calculating 
uses the same criterion as the third approach, but augments" 
the Abramovich solution for 2a, since the ambient turbulence 
affects both sides of the boundary layer ambient tur- 

bulence acts to decrease the radius of the inner core and 
to increase the outer wake radius). Therefore, for this 
approach 

db 


dx 

and 


21 




+ (2a) 


H^m 


CD-12) 




r 


21 


db/dx 


CD-13) 
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Streamline ri approach . - Figure D-2a shows Region I 
for the Abramovich solution. The three areas shown are the 
free stream flow, the potential core, and the boundary layer 
between the free stream flow and the potential core. Also 
shown is the streamline which passes through the initial 
wake boundary. It is assumed that the boundary layer develops 
on both sides of this streamline. In order to calculate the 
growth of the boundary layer in Region I, it is necessary 
to calculate the value of r , the radius of the streamline 
which passes through the initial wake radius at the end of 
Region I. The radius of the streamline, r , is determined 
by conservation of mass between the initial wake and the 
end of Region I. Therefore, 


1 

TTro^Ufl = 2-nfu(T)r dr 

*'0 

Since 11^ = 11 /m, the equation is written as 


(D-14) 


'0 


r 

0 ^ 


m 




= 2 I rU (r) dr 
0 


CD-15) 


The wind speed profile in Region I is given by equations 
(3-32) and (3-33). Multiplying equation (3-32) by (Uq-U^)/U^ 
and using the definition ofm given in equation (3-9) 
gives 


or 


1 U 


m U \m 
00 ' 


l)(l-n^-")^ 


rl m-1 1.52-1 

u = u_ I - + — (i-n ) 

mm 


(D-16) 


(D-17) 


From equation (3-33), recalling that r^ = 0 at the 
end of Region I Csince the end of Region I is defined as that 
point at which the potential core vanishes) and T 2 = ^ 2 ^ at 
the end of Region I, 


r = r22(l-D) 


(D-18) 
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and 


dr = do 


CD^19] 


Using equations (P’'17] , (_I>'18) , and (1^-19) in equation (D-IS) 
gives 



m 





Cl-'n)Ii+ C to'1) ' 




(D-20) 



Integrating gives 

2n^ 2n^ • ® n** n®-. ^ 

n + +_ — CD- 2 2) 

2 2.5 3.5 4 5 d 

’^s 

Since (rg/r 22 ) is a function of m only, is a function of 
m only. 

Equation (0-22) was solved numerically for four values 
of m as listed below. 

^ Ds 

1.5 0.390 

2.0 0.383 

2.5 0.378 

3.0 0.375 

value of rig = 0.38 is reasonable for the range of values of 
for wind turbines. From equation (D-18) 

r^ = 0.62r22 (D-23) 
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It is assumed that the wind speed profile given by 
equation C3*'32) is valid in the presence of ambient tur- 
bulence, With this assumption, the relationship of equation 
CD-23] is valid in the presence of ambient turbulence. Figure 
D-2b shows wake growth in the presence of ambient turbulence. 

Let be the distance from the streamline which passes 
through the initial wake boundary to the line between the 
boundary layer and the potential core. Then from Figure D-2a, 
the growth rate of due to mechanical turbulence (i.e., the 
Abramovich solution}'^is 

/dgA T 

(— ) = ^ (D-24) 

\dx /m C^H^m 

where is given by equation CD’'2) . 

The wake growth rate due to ambient turbulence is cx. 
Adding the ambient turbulence to the mechanical turbulence 
by the square root of the sum of the squares of the ambient 
turbulence and the mechanical turbulence gives 


d$^ 

dx 



(D-25) 


Since 6, = r at the end of Region I, the downwind extent 
of Region I is 


X 


H 


0 . 62r 22 





(D-26) 


Streamline T2 approach . - For the sixth approach for 
calculating the downwind extent of Region I, the portion 
of the boundary layer above the streamline which passes 
through the initial wake boundary is considered. Let 3o be 
the distance from the streamline which passes through the 
initial wake boundary to the line between the boundary 
layer and the free stream. Then from Figure D-2a, the growth 
in due to mechanical turbulence is 
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dg. 


^2l'^s 


CD-27) 


\dx /m 

Adding ambient turbulence in the manner described above 


gives 


dg. 

L 

dx 




Since 


^2 ^2l'^s 0-38r22 


the downwind extent o£ Region I is 


CD-28) 


CD-29) 


0.38r 




21 


' -I 


CD-30) 


0.5r7i approach . - For the seventh approach for calculat- 
ing the downwind extent of Region I, the boundary layer is 
assumed to develop about a line which bisects the boundary 
layer. The boundary layer is assumed to develop equally on 
both of its sides. In this case 

Vs = 0.5r2i CD-31) 

In this case, the line connecting the initial wake boundary 
and r at the end of Region I is not a streamline. The 
downwind extent of Region I is developed in the same manner 
as presented above for the previous two approaches. The 
results are identical for development based on g^ or on 
above because in this case g^ = g 2 * For both cases, 



CD-32) 
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It is noted that multiplication of the numerator and 
denominator of the right side of equation (J)-32] gives the 
identical result given by equations CD-12') and . 

Therefore, the b approach II and the 0.5r22 approach are 
mathematically identical. 


Comparison of Results 

Figure D-3 shows the donwwind eatent of Region I as a 
function of a, the wake growth rate due to ambient turbulence, 
for the first four approaches presented above. For Figure 
D-3, m = 3. Figure D^4 shows similar data for the last three 
approaches for calculating the downwind extent of Region I. 

For Figure D-4, m = 3. Figures D-5 and D-6 show information 
similar to that shown in Figures D-3 and D^4, respectively, 
but for m = 2. As mentioned previously, approaches 4 and 7 
are mathematically identical. 

Approaches 3 and 4 are identical, except that approach 

3 augments the Abramovich solution for a, whereas approach 

4 augments the Abramovich solution for 2a. From a physical 
point of 'view, the 2a approach is more reasonable than the 
a approach because ambient turbulence affects the growth 

of the boundary layer on both sides of the boundary layer. 

In Region I, mass continuity demands that the flow 
have an inward radial component. This is illustrated by 
the fact that the streamline which passes through the 
intital wake radius has an inward component. For approaches 
1 and 2, the wake growth rate is taken relative to the line 
r = r,,, which has no physical relationship to the flow 
direction or to the boundary layer. Relative to the flow 
and to the boundary layer, the line r = rp is a randomly 
drawn line. From a physical point of view, the measurement 
of the growth of the boundary layer from the streamline or 
from the line which bisects the boundary layer is more 
plausible . 

For these reasons, from a physical point of view, 
approaches 4, 5, 6, and 7 are the most plausible approaches 
for the calculation of the downwind extent of Region I. 

Figures D-4 and D-6 show that the results do not vary 
greatly between these approaches (recalling that approach 
4 is mathematically identical to approach 7 ) . For these 
reasons, approaches 4 and 7 were chosen as the approaches 
for calculating the downwind extent of Region I in the model. 



Radial coordinate, r Radial coordinate, r Radial coordinate 
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Downwind coordinate, x (x„l 

* '■ H-'m 


Figure D-1. - Angle used from Abramovich model to be 
augmented for ambient turbulence and 
used to calculate downwind extent of 
Region I 
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(a) Region I for Abramovich solution 
(no ambient turbulence) 



Figure D-2. Geometry of Region I of the wake. 
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b approach I 
approach 
b approach II 

T2 approach 


due to ambient 




Downwind extent of Region 




Q.l 0.2 0.3 0.4 

Wake growth rate due to ambient 
turbulence, a 


Figure D-5, - Comparison o£ first four approaches 
for calculating x„ for m = 2. 




0.1 0.2 0.3 

Wake growth rate due to ambient 
turbulence, a 


Streamline 

r^ approach 

0 . 5 approach 

— Streamline 
I T 2 approach 


Comparison o£ last three approaches 
for calculating x„ for m = 2. 
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